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This  report  is  concerned  with  a  proper  treatment  of 
disturbances  in  the  ionosphere.  It  develops  the  treat¬ 
ment  of  the  excitation  of  extra-low-frequency  (magneto¬ 
hydrodynamic-like)  wave-motions  in  an  infinite,  homogeneous 
plasma  imbedded  in  a  constant,  unidirectional  magnetic  field, 
for  which  the  electrical  conductivity  tensor  is  easily  es¬ 
tablished  and  methods  of  analysis  are  highly  developed. 

The  problem  treated  in  some  detail  is  the  radiation  by  charges 
traveling  along  the  magnetic  field  with  velocities  comparable 
to  the  characteristic  phase  velocity  (Alfven  velocity)  for  low 
frequency  waves.  The  problem  resembles,  to  some  extent,  the 
ordinary  Cerenkov  problem  in  which  charged  particles,  by  a 
collective  mechanism,  generate  visible  radiation  when  passing 
through  transparent  substances.  The  ELF  analogue  --  which  we 
call  the  magneto  gas  dynamic  analogue  of  Cerenkov  radiation  -- 
has  several  distinguishing  features.  For  one,  the  radiation 
is  most  important  in  a  narrow  cone  of  angles  centered  on  the 
magnetic  field  line  of  the  source.  In  addition,  the  mag¬ 
netic  field  of  the  wave  is  transverse  but  the  electric  field 
is  longitudinal.  The  disturbance  is  thus  not  circularly 
polarized  as  these  are  not  transverse  waves.  Also,  in  addi¬ 
tion  to  the  usual  Cerenkov  component,  for  which  the  particle 
velocity  must  be  greater  than  a  (the  Alfven  phase  velocity) , 
there  exists  a  weaker  component  originating  from  particles 
less  than  a. 

The  latter  portion  of  the  report  represents  a  preparatory 
step  in  the  development  of  a  Microscopic  Approach  to  plasma 
problems,  by  solving  the  first  two  equations  in  the  BBGKY 
hierarchy  in  two  noncorrelation  limits.  The  possibility 
of  deriving  the  dispersion  relation  for  a  sound  wave  pro¬ 
pagating  in  a  hard-sphere  gas  directly  from  the  linearized 
Boltzmann  equation  for  hard  spheres,  is  also  discussed. 
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ABSTRACT 


This  report  is  concerned  with  a  proper  treatment  of  disturbances 
in  the  ionosphere .  It  develops  the  treatment  of  the  excitation  of 
extra-low-frequency  (magnetohydrodynamic-like)  wave-motions  in  an 
infinite,  homogeneous  plasma  imbedded  in  a  constant,  unidirectional 
magnetic  field,  for  which  the  electrical  conductivity  tensor  is  easily 
established  and  methods  of  analysis  are  highly  developed. 

The  problem  treated  in  some  detail  is  the  radiation  by  charges 
traveling  along  the  magnetic  field  with  velocities  comparable  to  the 
characteristic  phase  velocity  (Alfven  velocity)  for  low  frequency 
waves.  The  problem  resembles,  to  some  extent,  the  ordinary  Cerenkov 
problem  in  which  charged  particles,  by  a  collective  mechanism,  gen¬ 
erate  visible  radiation  when  passing  through  transparent  substances. 
The  ELF  analogue  —  which,  we  call  the  magneto  gas  dynamic  analogue 
of  Cerenkov  radiation  —  has  several  distinguishing  features.  For 
one,  the  radiation  is  most  important  in  a  narrow  cone  of  angles  cen¬ 
tered  on  the  magnetic  field  line  of  the  source.  In  addition,  the 
magnetic  field  of  the  wave  is  transverse  but  the  electric  field  is 
longitudinal.  The  disturbance  is  thus  not  circularly  polarized  as 
these  are  not  transverse  waves.  Also,  in  addition  to  the  usual 
Cerenkov  component,  for  which  the  particle  velocity  must  be  greater 
than  a  (the  Alfven  phase  velocity) ,  there  exists  a  weaker  component 
originating  from  particles  of  velocity  less  than  a. 

The  latter  portion  of  the  report  represents  a  preparatory  step 
in  the  development  of  a  Microscopic  Approach  to  plasma  problems,  by 
solving  the  first  two  equations  in  the  BBGKY  hierarchy  in  two  non¬ 
correlation  limits.  The  possibility  of  deriving  the  dispersion  rela¬ 
tion  for  a  sound  wave  propagating  in  a  hard-sphere  gas  directly  from 
the  linearized  Boltzmann  equation  for  hard  spheres,  is  also  discussed. 
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CHAPTER  I 
INTRODUCTION 

The  work  reported  here  consists  of  two  distinct  efforts e  The 

first  effort,  will  generally  be  referred  to  by  the  heading  "Excitation 

Problem".  It  aims  at  explaining  the  extra -low- frequency  signals 

observed  to  occur  naturally ^  and  also  in  conjunction  with  high  al- 

(2) 

titude  nuclear  explosions  in  the  ionized  regions  of  the  atmos- 

(3) 

phere,  Only  the  simplest  model  ionosphere  is  used  in  this  treat¬ 
ment  —  a  uniform,  infinite,  homogeneous  two-fluid  plasma,  neutral 
overall,  embedded  in  a  constant  uniform  magnetic  field;  only  Coulomb 
and  electromagnetic  forces  are  considered;  only  the  linearized  aver¬ 
age  field  (Vlasov)  approximation  is  used  to  determine  the  dielectric 
and  conductive  behavior  of  the  plasma;  and  collisions  are  entirely 
ignored. 

The  second  effort  aims  at  improving  the  description  of  the  plas¬ 
ma.  This  is  done  by  treating  the  inter-particle  forces  in  a  higher 
approximation,  the  result  of  which  is  to  bring  in  collisions  and 
in  other  ways  modify  the  dielectric  and  conductive  properties  of  the 
plasma.  The  modes  of  oscillation  and  propagation  which  then  exist 
for  the  plasma  in  interaction  with  the  electromagnetic  field  are 
expected  to  provide  new  insights  into  plasma  behavior,  and  possi¬ 
bly  also  into  ionospheric  phenomena,  which  are  very  diverse  and  little 
understood.  Because  the  microscopic  treatment  of  the  plasma  is  the 
only  feasible  one  in  this  endeavor  we  shall  refer  to  the  second 
effort  by  the  heading  "Microscopic  Approach." 

Recent  experimental  and  theoretical  progress  in  geophysics 
has  greatly  clarified  the  gross  features  of  the  ionosphere  and 
the  exosphere.  This  progress  has  come  about  largely  through  direct 
exploration  of  the  earth’s  neighborhood  by  instrumented  satellites, 
and  high-altitude  nuclear  and  thermo-nuclear  explosion  experiments. 

The  satellite  experiments  measure  the  normal  features  of  the  several 
ionized  regions,  together  with  naturally  occurring  disturbances  in 
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these  regions.  The  explosion  experiments  measure  the  dynamical- 
response  properties  of  these  regions  ,  and  are  therefore  complementary 
to  the  others.  The  picture  of  the  earth's  ionized  atmosphere  which 
emerges  illustrates  the  competition  between  the  influence  of  the  earth's 
magnetic  field  on  the  charged  particles  entering  it  and  the  particles' 
influence  on  the  earth's  field.  Where  the  particle  densities  are 
relatively  small,  (such  as  in  the  Van  Allen  regions  and  in  the  arti- 
fically  created  electron  shells)  ,  the  charged  particles  follow  the 
characteristic  Stormer  orbits.  However,  where,  the  earth  magnetic 
field  is  sufficiently  weak  the  field  lines  get  dragged  along  (so  to 
speak)  by  the  particle  fluxes  (as  far  as  we  know  entirely  from  the 
sun) .  The  region  around  the  earth  which  (apart  from  disturbances) 
is  dominated  by  and  moves  with  the  earth's  magnetic  field  is  often 
referred  to  as  the  magnetosphere.  Because  the  axis  of  the  earth's 
dipole  field  lines  cross-wise  to  the  plasma  streaming  from  the  sun  — 
(referred  to  as  the  solar  wind)  ,  the  earth  is  effectively  shielded 
within  its  magnetosphere.  Charges  are  deflected  away  from  the  earth, 
or  trapped  in  the  field,  or  focussed  towards  the  polar  regions.  As 
a  consequence,  not  all  of  the  magnetic  disturbances  seen  outside 
the  magnetosphere  lead  to  magnetic  disturbances  at  the  earth's  sur¬ 
face. 

The  magnetosphere  appears  to  be  typically  some  8-10  earth-radii 
(Re)  distant  (from  the  center  of  the  earth)  on  the  sunward  side, 
and  perhaps  20-30  Re  distant  on  the  opposite  side  of  the  earth. 

Within  these  limits  we  may  properly  take  the  earth's  magnetic  field 
to  be  constant,  of  dipole  shape,  and  unaffected  by  the  ionized 
atmosphere  which  is  distributed  (somewhat  unevenly)  throughout  it. 

A  proper  treatment  of  disturbances  in  the  magnetosphere  —  which  is 
the  object  of  the  present  work  —  is  hardly  possible  by  analytical 
methods  if  one  wants  to  include  the  dipole  nature  of  the  ambient 
magnetic  field  as  well  as  the  variations  in  particle  densities 
and  particle  flux  densities.  This  is  particularly  true  when  one 
wants  to  study  oscillatory  and  wave-like  phenomena  at  such  low 
frequencies  that  the  relevant  wave-lengths  are  comparable  to  or 
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over-shadow  the  structure  and  spatial  variations  of  the  magnetic 
field  and  particle  distributions.  However,  so  much  remains  to  be 
learned  of  a  qualitative  and  semi-quantitative  nature  about  the  inter¬ 
action  of  particles  (natural  and  artificial)  with  the  magnetosphere, 
and  about  the  propagation  of  disturbances  generated  within  it  and  at 
its  surface,  as  to  justify  the  study  of  more  tractable  situations. 

In  particular,  in  the  present  study,  we  develop  the  treatment  of  the 
excitation  of  extra-low-frequency  (magnetohydrodynamic-like)  wave 
motions  in  an  infinite  homogeneous  plasma  imbedded  in  a  constant, 
unidirectional  magnetic  field  for  which  the  electrical  conductivity 
tensor  is  easily  established  and  methods  of  analysis  are  highly 
developed. 

The  problem  we  treat  in  some  detail  is  the  radiation  by  charges 
traveling  along  the  magnetic  field  with  velocities  comparable  to 
the  characteristic  phase  velocity  (Alfven  velocity)  for  low  fre¬ 
quency  waves.  The  problem  resembles,  to  some  extent,  the  ordinary 
(4) 

Cerenkov  problem,  in  which  charged  particles,  by  a  collective 
mechanism,  generate  visible  radiation  when  passing  through  trans¬ 
parent  substances.  The  ELF  analogue  —  which  we  call  the  magneto¬ 
gas-dynamic  analogue  of  Cerenkov  radiation  --  has  several  distin¬ 
guishing  features.  For  one,  the  radiation  is  most  important  in  a 
narrow  cone  of  angles  centered  on  the  magnetic  field  line  of  the 
source.  For  another,  the  magnetic  field  of  the  wave  is  trans¬ 
verse,  but  the  electric  field  is  longitudinal.  The  disturbance  is 
not  circularly  polarized,  (i.e.,  these  are  not  transverse  waves). 

In  addition,  besides  the  usual  Cerenkov  component  for  which  v,  the 
particle  velocity,  must  be  numerically  greater  than  a,  the  Alfven 
phase  velocity,  there  exists  a  weaker  component  originating  from 
particles  of  velocity  less  than  a. 

We  believe  that  some  of  the  naturally  occurring  low  frequency 
disturbances  noted  on  magnetogram  records  --  especially  those 
from  observatories  located  in  the  more  northerly,  magnetically 
active  areas  —  can  be  understood  in  terms  of  the  Cerenkov  mechanism 
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with  reasonable  estimates  of  the  particle  fluxes  even  though  at  low 
frequencies  the  effect  is  in  general  very  weak.  Possibly  of  even 
greater  promise  is  the  application  of  our  calculations  to  the  hydro- 
magnetic  signals  definitely  known  to  accompany  high  altitude  nuclear 
explosions.  Both  of  these  applications  will  be  presented  in  the  sequel 
to  this  report. 

Chapter  II  is  devoted  to  the  development  of  the  excitation  formu¬ 
lation,  a  development  which  proceeds  naturally  from  an  earlier  report 
by  Cantor,  Keilson  and  Schneider ^  that  considered  a  non-local 
(temperature-dependent)  model  of  the  ionosphere.  The  formulation 
given  there  is  carried  over  and  brought  to  a  useful  stage.  That 
report  and  this  are  both  motivated  by  a  desire  to  found  the  theory 
of  ionized  gases  on  a  secure  microscopic  (kinetic)  basis,  to  set 
aside  the  doubtful  macroscopic  equations  in  favor  of  less  doubtful 
'molecular'  equations  that  properly  include  the  interparticle  forces* 
For  the  problem  treated  --  the  zero-temperature  collisionless  plasma  -- 
there  is  no  practical  advantage  to  the  philosophy;  the  results  of  the 
two  approaches  are  identical  as  far  as  the  electromagnetic  response 
functions  are  concerned.  This  is  due  to  the  fact  that  in  the  earlier 
report  the  basic  kinetic  approximation  was  the  average-field  or 
Vlasov  approximation,  that  treats  each  particle  as  if  it  moved  in 
the  average  field  of  all  others.  The  consequence  for  the  zero 
temperature  plasma  (and  also  the  low-temperature  plasma)  is  that 
the  molecular  encounters ^  --  the  true  sign  of  'kinetic'  behavior  — 
were  neglected.  The  later  sections  of  this  report  explore  the  possi¬ 
bility  of  bringing  a  more  refined  description  into  the  theory  (as 
an  improvement  in  the  electromagnetic  response  function) . 

In  Chapter  III  the  linear  electric  field  response  to  an  impulse- 
current  source  is  computed  in  some  detail  for  the  frequency  range 
below  the  ion  cyclotron  frequency. ^  This  "radiation  Green's 
Function"  is  the  simplest  function  to  illustrate  the  directional 
properties  of  Alfven-like  waves  generated  in  a  compressible  con¬ 
ducting  fluid.  The  coupling  of  electromagnetic  and  kinetic  pro¬ 
cesses  is  by  way  of  the  Vlasov  average-field  approximation . 
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The  radiation  Green's  function  is  an  isotropic  matrix  of  functions 
exhibiting  five  distinct  kinds  of  response  to  a  point  excitation. 

First,  there  is  the  excited  field.  This  occurs  even  in  the  absence 
of  an  external  magnetic  field.  It  exhibits  isotropic  propagation 
for  those  compone^s  whose  frequencies  are  above  the  plasma  frequency, 
and  damping  for  frequencies  below  the  plasma  frequency.  Second, 
there  is  a  signal  which  is  propagated  at  the  Alfven  phase-velocity 
equally  in  all  directions.  This  signal  is  characteristic  of  the 
compressible  fluid.  Third,  there  is  a  signal  which  is  propagated  in¬ 
side  a  double-cone  centered  on  a  magnetic  field  line  with  apex  at  the 
source  point,  and  which  is  damped  outside  of  the  cone.  The  narrowness 
of  the  cone  varies  with  the  emitted  frequency,  becoming  most  narrow 
at  the  lowest  frequencies.  It  is  this  which  corresponds  to  the  Alfven 
mode  (propagation  along  the  magnetic  field)  observed  in  incompressible 
conducting  fluids.  The  modification  due  to  compressibility  is  such 
as  to  make  this  part  of  the  radiation  field  singular  on  the  surface 
of  the  cone.  The  fourth  and  fifth  kinds  of  behavior  are  caused  by 
the  interference  or  competition  of  the  Alfven-mode  with  each  of  the 
first  and  second  kinds  of  behavior.  When  the  interference  is  with 
the  (first)  damped  mode,  the  propagation  in  the  cone  is  lost.  When 
the  interference  is  with  the  (second)  isotropically  propagating  mode, 
the  cone  boundary  becomes  diffuse  and  the  radiation  is  enhanced, 
especially  along  the  magnetic  field  direction. 

The  Normal  Cerenkov  Effect  is  examined  in  Chapter  IV,  and  several 
of  the  mathematical  details  that  will  be  needed  to  calculate  the 
Magneto-Gas-Dynamic  Analogue  of  the  Cerenkov  Effect  are  presented. 

This  calculation  —  the  main  effort  of  the  report  —  is  given  in 
Chapter  V.  We  have  already  stated  some  of  the  ways  in  which  the 
Magneto-Effect  differs  from  the  Normal  Effect.  In  both  cases  the 
radiation  appears  only  at  a  specific  angle  with  respect  to  the  motion 
of  the  charge  (same  as  the  direction  of  the  ambient  magnetic  field) . 
This  angle  depends  only  on  the  velocity  (v)  of  the  charge  relative 
to  the  Alfven  velocity  (a)  .  But  whereas  normally  v  must  be  greater 
than  a,  now  it  may  also  be  less  than  a,  though  this  "anomalous" 
effect  is  negligibly  small  compared  to  the  main  effect.  (However, 
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it  has  important  theoretical  consequences,)  The  radiation  is,  in 
both  cases,  linearly  polarized,  but  unlike  the  normal-effect  the 
electric  field  of  the  radiation  is  only  along  the  ambient  magnetic 
field  (Bq)  while  the  magnetic  field  of  the  radiation  is  only  per¬ 
pendicular  to  it.  The  "radiation therefore,  is  in  a  "hydromagnetic" 
mode.  Furthermore,  the  rate  of  radiation,  as  a  function  of  angle,  is 
greatest  near  the  magnetic  field  direction,  when  v  is  very  close  to 
but  slightly  different  from  a.  However,  it  goes  to  zero  when  v  =  a. 
T^e  latter  behavior  is  also  true  of  the  normal  Cerenkov  radiation, 
but  in  the  magnetic  case  it  comes  about  through  a  cancellation  of 
two  terms  which,  when  v  is  not  too  close  to  a,  give  rise  separately 
to  the  "main”  and  "anomolous"  effects. 

Chapter  VI  is  given  over  to  a  discussion  of  the  temperature 
corrections  that  may  be  expected  to  play  a  part  in  ELF  propagation, 
and  they  are  found  not  to  be  important  in  the  examples  we  have  con¬ 
sidered,  though  this  is  by  no  means  always  the  case.  More  will  be 
said  on  this  point  in  the  sequel . 

The  final  two  Chapters  represent  a  preparatory  step  in  the 

development  of  a  Microscopic  Approach  to  plasma  problems.  Chapter  VII 

shows  how,  in  terms  of  a  sequence  of  distribution  function  f  ^  (rvt)  , 

(2) 

f  (rv,r,v,t) ,  etc.,  describing  the  average  properties  of  an  assem¬ 
bly  of  particles  in  successively  greater  detail,  one  can  formulate 

(7) 

the  microscopic  laws  of  the  system.  The  resulting  equations  are 
solved  in  two  non-correlation  limits  —  the  first  leads  to  the  Vlasov 
approximation?  the  second  leads  to  the  (generalized)  Boltzmann 
Equation,  which  includes  the  Vlasov  limit.  To  illustrate  its  relation 
to  the  more  customary  Boltzmann  Equation  we  apply  it  to  the  case  of 
hard-spheres,  actually  using  a  limiting  form  of  interaction  poten¬ 
tial  to  describe  collisions.  In  the  sequel  we  shall  apply  it  to  the 
two  fluid  plasma  with  special  attention  to  the  rederivation  of  the 
plasma  conductivity.  In  Chapter  VIII  we  discuss  the  possibility  of 
deriving  directly  the  dispersion  relation  for  a  sound  wave  in  the 
hard-sphere  gas  rising  the  linearized  Boltzmann  Equation  for  hard 
spheres,  and  find  that  the  difficulties  are  great,  even  in  the  one¬ 
dimensional  approximation  (which  tas  been  treated  with  some  success 
(8) 

by  other  methods) . 
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Summary  and  Outlook 

To  summarize:  an  improved  understanding  of  ionospheric  (and 
other)  plasma  phenomena  requires  (a)  detailed  calculations  of  the 
electromagnetic  response  to  realistic  excitations  and  (b)  develop¬ 
ment  of  methods  for  including  in  the  plasma  description  other  conse¬ 
quences  of  the  collective  behavior  than  just  the  average-field, 
space-charge  forces.  With  regard  to  (a)  we  may  say  that  although 
there  has  been  a  great  deal  of  attention  to  the  magneto-ionic 
modes,  the  more  difficult,  and  tedious,  excitation  problems  have 
been  generally  sidestepped.  In  the  present  work  we  tackle  the 
problem  of  ELF  wave-motions  generated  by  charges  moving  along  the 
magnetic  field,  and  in  the  sequel  we  shall  include  an  additional 
spiralling  motion  to  duplicate  the  actual  behavior  of  charges  injected 
into  the  earth's  field  without  preferential  direction. 

With  regard  to  (b)  we  can  say  that  our  aim  is  to  include,  with¬ 
in  the  framework  of  a  general  formalism,  the  modifications  of  the 
electrical  susceptibility  and  conductivity  tensors  which  are  a  con¬ 
sequence  of  the  ’’higher  correlations”  of  particle  behavior.  The 
most  important  of  such  modifications  will,  for  our  purposes,  most 
likely  be  those  due  to  internal  collisions.  For  the  ionosphere, 
another  kind  of  collision  is  important  —  external  collision,  i.e., 
of  electrons  and  ions  with  the  neutral  molecules  of  the  gas.  One 
of  our  motivations  in  treating  the  hard-sphere  gas  is  to  provide 

some  idea  for  the  analytical  treatment  of  external  collisions.  In 
(3) 

earlier  work  a  Fokker-Planck  model  was  used  for  this  purpose, 
but  even  though  it  had  several  nice  features  it  led  to  no  practical 
improvement  over  the  simpler  methods  used  in  the  Appleton-Hartree 
theories  of  ionospheric  conductivity  of  a  cold  plasma. 

Finally,  in  the  next  report  we  shall  present  the  application 
of  the  results  of  our  calculations  to  both  natural  ionospheric  ELF 
phenomena,  and  to  the  artifically  created  ELF  signals  of  the  various 
high-altitude  nuclear  explosion  experiments „ 
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CHAPTER  II 

FORMULATION  OF  THE  EXCITATION  PROBLEM 

In  line  with  our  program  to  deduce  the  macroscopic  plasma  behavior 
from  molecular  laws,  we  shall  write  Maxwell’s  equation  for  vacuum 
processes  and  deduce  the  dielectric  and  conductive  parameters  from 
the  net  response  of  the  charges  and  currents  in  the  vacuum  to  ap¬ 
plied  fields. 


V  X  E  = 
** 

_ 

at 

(2.1) 

V  •  E  = 

P/£0 

(2.2) 

V  x  B  = 

I2£  +  a  j 

2  at  +  ^ojI 

c 

(2.3) 

II 

t> 

0 

(2.4) 

Throughout  this  work  we  shall  employ  an  implicit  vector  nota¬ 
tion.  No  special  mark  will  designate  a  vector,  and  it  will  be  under¬ 
stood  that  the  electric  field  E,  the  magnetic  field  B,  the  electric 
current  density  j,  and  the  gradient  operation  V  are  vectors,  while 

the  electric  charge  density  p  is  a  scalar  and  the  time-derivative 

3 

is  a  scalar  operation.  Furthermore,  p  and  the  separate  compon¬ 
ents  of  E,  B  and  j  are  functions  of  the  vector  position  r  =  (x,y,z) 
and  the  time  t.  The  ’’dot”  denotes  scalar  product,  the  ’’cross”  vector 
product . 

-7  -9 

The  constants,  jjt  =  4tt  x  10“  henries/meter  and  e  ^  1/3677  x  10 

farads/meter,  are  related  by  the  velocity  of  light  (c  =  3  x  10  meters/'se 

according  to 

e^xo  =  1/c2  .  (2.5) 
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Consequently,  charge  conservation  follows  from  the  above: 


jj£+  v 

3t 


j  =  0 


(2.6) 


Finally,  the  Lorentz  force  (a  vecotr  F)  on  a  particle  of  charge  q, 
velocity  v  (a  vector) ,  at  rt  (space-time  coordinates)  takes  the  form: 


F  =  q(E(rt)  +  v  x  B(rt)) 


(2.7) 


[The  applied  fields  will  be  a  constant  magnetic  field  and  the 
field  of  a  prescribed  charge  and  current  source,  or  as  we  shall  call 
it  an  external  source 0]  That  part  of  p  and  j  which  is  not  prescribed 
will  necessarily  have  been  induced  in  the  plasma: 


P  = 


ext 

P 


+  P 


ind 


j 

a* 


J 


ext 


+ 


j 


(2.8) 

(2.9) 


The  induced  charge  and  current  densities  are  due  to  the  motion 
of  the  constituent  particles,  the  equations  of  motion  for  which  con¬ 
tain  the  total  fields  E  and  B. 

The  plasma  we  have  in  mind  is,  in  its  undisturbed  state,  uniform, 
homogeneous,  and  in  equilibrium.  The  presence  of  a  constant  mag¬ 
netic  field  disturbs  only  the  isotropy  of  the  plasma’s  response,  but 
not  its  homogeneity  or  equilibrium.  This  term  is  characterized  by  no 
net  charges  or  currents,  no  electric  fields,  and  no  magnetic  fields 
other  than  Bq.  There  will  be  fluctuations  of  these  quantities  be¬ 
cause  of  the  molecular  nature  of  ions  and  electrons,  but  we  shall 
ignore  these  fluctuations  relative  to  the  disturbances  which  we  do 
treat.  For  our  purposes  the  plasma  will  be  described  by  entirely 
continuous,  slowly  varying  functions  (slow  relative  to  the  fluctua¬ 
tions)  . 


10 


S-2023-1 


A  disturbance  of  the  plasma  will  be  characterized  by  a  non¬ 
vanishing  E  and  B ,  together  with  accompanying  charge  and  current 
densities.  We  should  suppose  that,  in  the  sense  of  Ohm's  law,  there 
is  a  linear  relation  between  jln<^  and  JS,  B,  We  shall  discover 

that  so  long  as  the  only  source  of  anisotropy  is  B  ,  only  E  enters 
these  relations.  They  are: 


oo 


00 


pind(rt)  =  -  V-e  f  (d3r ' )  f  dt'  X(r-r* ,t-t’)  «E(r' t') 

r  m  O  J  /+  J  m  m  1  m  ** 


-OO 


-OO 


(2.10) 


CD 


.ind, 

J  (rt) 

A*  ** 


-  J 


(d3r ' 


-OO 


oo 


J 


-oo 


dt' 


cr  (r-r ' 


,t-t')  •E(r't') 


(2-11) 


The  quantities  X  (susceptibility  tensor)  and  cr  (conductivity  ten¬ 
sor)  will  be  understood  to  be  tensors,  and  dyadic  notation  will  be 
employed.  If  we  use  i,j  ,]£,•••  =  1,2,3  to  distinguish  the  components 
of  the  vectors  and  tensors?  then  in  more  detail  we  can  write: 


y  V.  f  (d3r ' )  f  dt'  X. . (r-r ' ,t-t 1 ) E . (r ' t ' )  (2.12) 

Lj  1  J  J  1J  jm.  Am  7  J 

ij  r'  t' 


jjnd(rt) 

m  1  * 


j 


dt’  a,  . (r-r' E. (r 1 t ' ) 
1J  m  A*  J 


(2.13) 


These  are  non-local  relations:  p  and  j  at  one  place  and  time 
depend  on  E  for  all  other  places  and  times.  Actually,  microscopic 
causality  implies  that  only  fields  at  earlier  times  determine  p  and 
j.  We  may  express  this  by  the  demand 


3 

X(t-t')  =  0  t  <  t' 

^f(t-t')  =  0 


(2.14) 


t  <  t' 


(2.15) 
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and  continue  to  extend  all  the  integrals  from  -ooto  +oo.  Thus  the 
induced  charge  and  current  are  given  by  fourfold  convolutions  of  the 
electric  field  with  tensors  characterizing  the  plasma's  responsiveness. 
That  X  and  <3  are  functions  only  of  r-r  r  and  t-t 1  is  true  because  the 
equilibrium  state  is  independent  of  r  and  t.  It  is  especially  conven¬ 
ient,  in  such  cases,  to  use  Fourier  integral  transforms  since  the  trans¬ 
form  of  a  convolution  is  simply  the  product  of  the  transforms  of  the 
factors . 

For  notation  we  shall  use  the  same  symbol  for  the  transform  as  for 
the  original  function,  and  shall  use  a  standard  conventions 


f  (rt) 


/ 

-CD 


(d3k) 

(2TT)  3 


% 

V 


-OO 


ch)  ik*r-icot 

2?  S  *  - 


f(kjo) 

A* 


(2.16) 


Then 


pind(ko)  =  -ik 


’£  )C(kuo) 

O  4* 


■E  (kco) 

A* 


(2.17) 


jind(kco)  -  cf(ka))  cE(kco) 

Ji  *  **  **  * 


(2.18) 


As  we  may  say,  these  relations  are  ’’local  relations  in  Jc -co-space . 11 
In  ordinary  parlance,  oo  is  the  (angular)  wave  frequency  and  k  the 
wave-vector  of  the  kco  component  of  the  electric  field. 

From  conservation  of  charge  there  emerges  an  important  relation 
^  ^  ext  ext 

between  X  and  cf.  Since  p  ,  j  separately  obey  the  conservation 
equation,  p^nc^,  j^n<^  do  so  too.  Substituting  (2.10)  and  (2.11)  into 
(2.6)  ,  we  realize  that  there  is  sufficient  arbitraries  in  E(rt)  to 

’  m  » 

conclude  that 


V 

r 


•cr(r-r ' 


,t-t ' ) 


(2.19) 
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or 

k^(kw)  =  -ie  6ok*^(koo)  (2.20) 

*  0  4+  ** 

in  the  everpresent  dyadic  notation. 

From  Maxwell’s  Equation  we  deduce  the  following  second-order 
equation: 

2 

,  n2  .  1  9  x  _  In/  ext  .  indx  5  ,  .ext  f  .indx 

(-V  +  — 2  - »)  J  =  -  —  V(p  +P  >  ■  at(J  +  J  ) 

c*  gt^  eo  °  0r—  — 

(2.21) 


or, 


or, 


(k2  -  ^)E  =  -  J=(pext(to)  +  pind(kw)) 

Z  6  to* 

C  O 

+  ICOLL  (jext(km)  +  jind(kto) )  (2.22) 

O  **  to 


(k2  -  %:)  E  +  k(k*X*E)  -  iu  ootf-E  =  -i  (k  ® -  -  u  cojext 

Z  to*  tot  **  O  to*  tom  £_  o  ** 

C  O 


(2.23) 

making  use  of  (2.17)  and  (2.18) .  Because  of  the  conservation  laws, 

ext  ext 

j  contains  whatever  significant  information  is  in  p  .  (The  con- 

ext 

verse  is  not  true.)  An  equation  containing  j  only  is  remarkably 

«+  2  2  2 

symmetrical.  Multiply  (2.23)  by  k*  ,  and  factor  out  (k  -  co  /c  ) 
after  using  (2.6)  and  (2.20)  : 

k*E  +  k-X-E  =  -  i  p6Xt/e  •  (2.24) 

<*»  m  mm  m.  O 
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Hence  upon  eliminating  p  ,  we  get 


(k2  -  ^)E  -  k(k-E)  -  ip.  aicr.E 

2  *+  m  m  O  0k 

C 


ext 


(2.25) 


Let  us  introduce  the  tensor  X  by  the  relation 


a  =  -  it  wx 

O 


(2.26) 


Then 


[<*2 '  *  '  -i-  ‘  ?  *  •  i  - 


(2.27) 


and 


3  * 

k*X  =  k.X  . 

*m  ««■» 

* . „?i,  * 


(2*28) 


The  symbol  1  is  the  unit  dyadic?  X  is  a  dyadic;  1*E  =  E,Xk-E  =  £(*•£)  • 

The  quantity  in  brackets  in  (2*27)  is,  in  Cartesian  components, 

a  matrix,  whose  inverse  allows  us  to  solve  for  E,  The  result,  when 

7  .  ext 

transformed  back  to  jrt- space  is  the  field  excited  by  the  current  • 

The  formal  solution  is 


!  (koo)  =  ip-  oo 


(k2  -  t  -  kk  -  ^  X  (ka3) 

c  c 


jext(ku))  .  (2.29) 

0+ 


Let  us  evaluate  the  matrix  before  j 
Introduce  the  scalars 


.ext 


and 


P  _  2  2.2 

k  -  oo  /c 


o>2/c2 

k2  -  *2/c2 


(2.30) 


(2.31) 
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Then 


A  = 


ik2  -  4  *  -  wc  -  <4 1 

C  C 


* 

1 


3 

M- 


-  0kk  - 


a  X 


(2.32) 


A 

We  first  treat  1  -  aX  as  a  unit,  and  formally  expand  the  denominator 
in  powers  of  {5 : 


* 

A  = 


p 


% 

i  - 


1 

aX 


* 

1 


? 

5 — 1 

1  -  aX 


(2.33) 


The  first  term  is  irreducible?  the  second  term  is  the  outer  prod¬ 
uct  of  two  vectors?  the  third  term  is  like  the  second,  except  for  a 
scalar  factor?  each  higher  term  contains  another  identical  scalar 
factor.  ^  can  therefore  be  written  as  two  terms,  by  summing  the  geo¬ 
metrical  series  of  scalar  factors. 


(2.34) 


It  has  been  necessary  to  use  the  same  symbol  1  for  the  unit 
matrix  and  the  scalar  unity.  The  latter  appears  only  in  the  large 
parenthesis,  which  is  the  sum  of  the  series  of  scalars: 


1  +  (k 


jf 


•o- 


k)  +  ( 


)2+  ( 


)3  + 
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We  shall  introduce  the  symbol 


61 

*  *  = 


1  -  aX 


3 

((3  times  inverse  of  1 


(2.35) 


for  the  important  matrix  that  appears  four  times. 

±  ZZ  M  •  kk  •  M 

A  =  M  +  - **** —  .  (2.36) 

1  -  k»M°k 

m*  -44tc 


(The  denominator  is  a  scalar,  so  we  needn't  worry  about  its  position 
in  the  formula.  But  we  shall  see  to  it  that  vectors  and  matrices  are 
arranged  according  to  dyadic  notation.  It  will  then  always  be  clear 
how  to  convert  them  to  component  form.) 


To  proceed  further  will  require  a  knowledge  of  X,  i.e.,  essentially 
the  tensor  conductivity  of  the  plasma  in  its  ambient  magnetic  field. 

So  long^as  the  plasma  is  at  zero  temperature  and  contains  no  drifting 
parts,  X  is  identical  to  X.  That  is,  the  conductivity  and  suscepti¬ 
bility  tensors  contain  identical  information.  Referring  back  to  (2.26), 
this  means  that 


cf  =  -  ie  ooX 


(2.37) 


which  is  a  stronger  condition  than  charge  conservation  (2.20). 

When  we  ask  for  the  low  frequency  waves  and  oscillations  excited 
by  external  currents,  X  simplifies  immensely.  Let  us  write  (2.29)  in 
the  form 


E(koo)  =  in,  ooA(koo)  •  j 

*«*  ”  O  m  ** 


(2.38) 


ext 

If  j  contains  only  low  frequencies,  E  will,  in  fact,  contain  the 

^*1  “J* 

same  low  frequencies.  What  this  means  is  that  only  A ' s  correspond¬ 
ing  to  these  low  frequencies  are  relevant  for  the  solution  JE.  Hence, 
if  we  anticipate  the  later  use  of  either  low-frequency  excitations 
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or  the  low-frequency  coniponents  of  JE  and  J3,  it  will  suffice  to  sim¬ 
plify  A  according  to  its  frequency  dependence.  The  frequency  range  we 
have  in  mind  is  below  the  ion  cyclotron  frequencies,  but  above  colli¬ 
sion  frequencies.  This  is  the  Alfven,  or  magneto-hydrodynamic,  range, 

and  is  distinguished  by  one-dimensional  wave  propagation  along  the 

it 

magnetic  field  lines.  X  must  necessarily  be  highly  anisotropic  to 
achieve  such  a  state  of  affairs.  But  the  simplifying  feature  in  this 
extra  low-frequency  region  is  that  x  is  diagonal.  Hence ,  ^ 
is  diagonal,  and  M  is  diagonals  simply  the  matrix  of  reciprocals. 

In  a  Cartesian  coordinate  system,  let  $  be  the  magnetic  field 
direction.  Then  the  non-vanishing  components  will  be ^ ^  9 ^ 


X  =  X 
XX 


yy 


■i 


»2. 

_Ei 

i 


(2.39) 


Xzz  =  " 


Y  X 

4  LI 


JBk 


oo 


(2.40) 


The  sum  is  over  the  various  components  of  the  plasma;  for  example, 

electrons  and  H+  ions,  oo^  is  the  plasma  frequency,  one  for  each 

component;  Q  is  the  gyro-frequency,  one  of  each  component.  We  shall 

let  the  total  plasma  frequency  be  oo  ; 

P 


i 


(2.41) 


It  is  approximately  the  same  as  the  electron  plasma  frequency.  Now 


Mxx  ~  Myy  1  -  aX. 


(2.42) 


xx 


17 


S-2023-1 


M 


zz  1  -  aX, 


(2.43) 


zz 


The  Alven  velocity  appears  in  M^;  it  is  designated  by  the  symbol  a 


2 

a  = 


1  + 


/  2\ 
rnf CO  \ 

V[  -J2  | 

iWii 


(2.44) 


Hence , 


-  M 

T  -  xx 


k2  _  a>2/a2 


(2.45) 


M  =  M  =  — r- 
z  zz  ,  2 


k^  +  o)2/c2  -  co2/ c2 
P 


(2.46) 


In  this  special  case,  M*k  =  k*M,  a  vector  with  the  components 


k*M  =  (k  Mm,  k  M_ ,  k  M  ) 
sm,  x  T’  y  T’  z  2 


(2.47) 


Further , 


k«M*k  =  (k2  +  k2) M_  +  k2M 
**,  x  y  T  z  z 


(2.48) 


We  can  let  kT  be  that  part  of  the  vector  k  which  is  transverse  to 
the  magnetic  field  lines: 


kT  =  (kx’ky’0) 


k2  -  k2  +  k2 
kT  x  +  Ky 


(2.49) 
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Then 


D  = 


1  -  k-^.k  =  1 


co‘  ..2 

(— E.  _ 

'  2  2' 
c  c 


k 


2 

T 


-  oo2/a2 


k2  +  oo2/c2  -  oo2/c2 
P 


MT^z 


(2.50) 
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CHAPTER  III 

THE  RADIATION  GREEN'S  FUNCTIONS 

We  now  turn  to  the  problem  of  determining  the  radiation  fields 
caused  by  a  single  charged  particle  passing  through  the  ionosphere* 

The  eye,  or  any  other  receiving  instrument,  is  sensitive  to  the 
frequency  components  of  the  resultant  radiation*  The  frequency  com¬ 
ponents  of  the  electric  field  at  a  point  r,  due  to  an  excitation  at 
a  point  rT ,  are  written  as  a  convolution  integral 


E(r,co)  =  i  co  [  ?(r-r',  co)  •  j ext  (r co)  d3r '  (3.1) 

mm7  °  J  m  4m  mm 


where  A  is  the  matrix  element  defined  in  Chapter  II,  which  appears 
explicitly  in  Equation  (2*36)  . 

Our  immediate  task  is  tj£  determine  A(r-r'  ,00)  .  The  elements  of 
A  in  wavenumber  Oc)  space,  A(k,oo)  can  be  obtained  directly,  as  was 
shown  in  Chapter  II,  Equation  (2.30) . 

* 

Then,  the  elements  of  the  matrix  A  may  be  written 


Ai*  -  Mi*  +  5  <Hij  kj  V 


(3.2) 


where 


D  a 


co2 

_£ 

2 

c 


2  k2 

,2  co  n  ^  2  T. 
X  -  -2  (i  +  c  -3) 

a  00 

P 


m^2 


(3.3) 


is  the  limit  of  Equation  (2.50)  when  co  <<  co  .  Throughout  this  paper, 

2  2  2  2  P  2  2 

we  have  neglected  co  /c  but  not  co  /a  with  respect  to  co^/c 

Explicitly 
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A 

_ 

M_ 

+  h2  mJ 

(3.4) 

XX 

T 

D  x  T 

A 

M_ 

(3.5) 

yy 

T 

D  T  y 

A 

M 

(3.6) 

zz 

z 

D  z  z 

A 

A 

1  2 

=  tM:  k  k 

(3.7) 

xy 

yx  u  t  x  y 

A 

A 

=  h  M_  M  k  k 

(3.8) 

xz 

ZX  D  T  Z  X  z 

A 

A 

=  i  M_  M  k  k 

(3.9) 

yz 

zy  d  t  z  y  z 

It 

is  readily 

seen  from  Equation 

(3.2)  and  from  the  nature  of 

that 

the  tensor 

a  is  symmetric. 

Substituting 

in  (3.4)  -  (3.8)  the 

expressions  for  MT,  Mz  and 

M, 


we  obtain 


Axx  -,.2  ...2  2 


1 -  +  c 


kt-  4*  , 

-  oo  /a  a), 


^yy  .2  2  ,  2 

k  -  oo  /a  oo, 


Z2  “  ,  2  2,2 

k  +  oo^/c 


2 

p<-4{ 

a  ' 

2 

c 

f  k2> 

,  ,  2  T 

1  +  c  ~2 

V  % 

k2 

V 

“p  k2  -  <4\ 

(..  2  kT 

a 

V  “p 

OJ 

CM 

C 

z 

2  o 

p  X2  -  ^1 
z  2 

a 

\  p. 

1  - 


.2,2  2,2. 

to  /a  +  a>^/c  ^ 


1  - 


k2  +  co2/c2 


fco2/a2  +  co  2/ c  2 ) 
W 


(3.10) 


1  + 


,2  +  co2/c2 

\r 


fco2/a2  +  co2/c2) 

H. 


(3.11) 


k2  -  co2/a2 


(3.12) 
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(oa2/c 


2  ,  2/2. 
+  co  /c  ) 
P _ _ 

^  2/2 
+  CO  /c 

p 


(3.13) 


(3.14) 


(3.15) 


^f(r-r',co)  is  then  obtained  from  A(k,co)  by  means  of  the  Fourier  integral 
theorem 


if(R,co)  =  r5A(k,co)  e  **' '*‘d3k/(27r)  3 


(3.16) 


where  JR  =  jr-rJ  . 

Note  that  in  the  expressions  for  ,  there  appear  five  different 
denominators ,  four  of  which  have  singularities  which  might  contribute 
to  radiation. 

We  shall  be  interested  in  integrating  these  five  expressions  over 
all  k-space.  We  need  concern  ourselves  only  with  the  denominators 
because  the  kx,  ky  and  kz  appearing  in  the  numerators  may  be  taken 
outside  of  the  integrals  as  appropriate  partial  derivatives  since 

0 

kx  transforms  to  -i  ^  ,  etc. 
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Let  us  obtain  the  Fourier  transforms  of  these  five  denominators 
which  we  shall  call  G1  to  inclusively.  It  should  be  noted  that  in 
later  calculations,  only  those  terms  that  fall  off  as  1/r  shall  be  retain¬ 
ed,  as  the  main  interest  is  in  the  radiation  parts  of  the  fields. 


Thus,  the  transformed  matrix  elements  take  the  form 

,  2 


2  A  2 
c  6 


XX 


Ayy  =  Gi 


Azz  ~  G2  "  .2 


0)2 

p 

ax2 

2 

c 

92 

2 

00 

P 

Q> 

ro 

2 

c 

92 

Is  9 ' 


r  6x2  5 


r1’ 


9y 

,  2 


2  5 


to  9z 

P 


2°3  -  (1  +  r’^G4 


(3.17) 


(3.18) 


(3.19) 


where 


xy  yx 


A  =  A 
xz  zx 


A  =  A 
yz  zy 


c2 

92 

'to2 

p 

9x9y 

2 

c 

92 

2 

9x9z 

c2 

92 

"to2 

p 

9y9z 

y'  9x9y 


(3.20) 


(3.21) 


(3.22) 


31  =  J 


d3k  eik>R 
(2n)  3  k2  -  to2/a2 


(3.23) 


»  ,3,  ik*R 

r  _  r  d  k  _ e _ 

^2  “  J  ,0  .  3  .2  ^  2,2 

(2ir)  k  +  to  /c 


(3.24) 
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where 


-J 


ik»R 


d3k _ 

12,7)3  v2  a! 

\  -  2  -  r 


-J 


d3k 

(2tr) ' 


ik»R 


e 


<*2  -  a!)  <*2  -  4 

a  a 


k2 

Y 


J 


ik-R 


d3k _ 

( 2tt)  3  _  os2  „  2 

(k2  +  -§)  (k2  -  Sj 
c  a 


£ 


co  2 

T  -  2  2 

CO  c 


(3.25) 


(3.26) 


(3.27) 


(3.28) 


We  must  specify  the  manner  of  evaluating  the  singularities  in  the 
inverse  transforms  and  G^  in  order  to  have  them  correspond 

to  real  physical  problems.  The  method  is  simple.  An  ambiguity  in  the 
inverse  transform  originates  from  the  use  of  a  Fourier  integral  not 
absolutely  convergent.  One  must  supplement  the  original  integrand  with 
a  damping  factor.  For  example, 

oo 

G (go)  =  J  dT  elaVr  G (T)  e~e  e  >  0  (3.29) 

-OO 

where  T  =  t-t  ’  is  the  usual  choice.  In  addition,  we  know  that  G(T) 
relates  the  field  at  time  t  to  a  source  at  time  t’  .  Causality  im¬ 
plies  that  G(t-t')  =  0  if  t  <  t’  (T  <  0)  .  Hence, 
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oo 

G  (oo)  =  f 

V 

O 


dT  e 


iooT  _-£  T 


G(T) 


dT  ei(oH-ie)T  G(T) 


(3.30) 


Thus,  for  causal  functions,  the  necessary  modification  is  to  replace 

G  (oo) — ►  G  (a>+  i  £ )  .  (3.31) 

Then  the  integral  (3.30)  is  absolutely  convergent,  and  in  the  final 
calculation,  one  can  let  £  — +  0  after  it  has  been  used  to  define 
the  singularities . 

Evaluation  of  G^ 


(3.32) 
(4) 

function  propagator 

(3.33) 


Gi(JS-“>  -  J 


e^'B-  la3*) 

k2  -  co2/a2  (2t r)  3 


This  is  the  well-known  radiation  Green's 
i  (oo/a)  R 


G1  - 


4fR 


where  R  is  the  magnitude  of  R  . 

Since  the  solution  of  G^  is  known,  it  is  advantageous  to  put 
the  integral  in  Equation  (3.32)  into  another  form  which  will  be 
encountered  repeatedly  in  future  calculations. 

The  denominator  of  the  integrand  in  (3.32)  can  be  "exponentiated." 
If  as  discussed  above,  the  00  in  the  denominator  is  replaced  by  (oo+ie)  , 


i  _  f  f  ids  exp  f-is (k2-  {-  7e  ))  exp (-2  — £  s)  (3.34) 

1  J  (2ir)  3  t  l  a2  J  a2 
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2ecos 
2 

The  e  a  serves  as  a  damping  factor  and  makes  the  integral 
(3,34)  absolutely  convergent. 

Letting  e — *0,  all  terms  involving  e  can  be  disregarded  and 


oo  2 

3^  =  J  - ^  J  ids  exp  -i  (k2  -  — ^)  s 

( 2rr)  a 


(3.35) 


The  k-integral  is  evaluated  by  completing  the  square  in  the 
exponential 


1  =  ^  UV  exp(is  ^2)exP(i  J  exp[-is(k  +  fi>  ] 


,3, 
d  k 


.-1/2  .2  2 
X~2/2  J  ds  s"  exp(i  ^  s)exp(i  j^) 

QU  q  a 


(3.36) 


However,  since  is  known,  the  solution  for  this  particular 
type  of  integral  in  (3.36)  can  be  obtained. 

In  general  for  A  >  0,  B  >  0. 


oo 


J  ds  s~3/2  eisA2  e1  4s  =  ^-g-A- 


B 

2 


(3.37) 


This  result,  Equation  (3.37) ,  will  be  used  in  subsequent  calculations. 


Evaluation  of  G, 


The  second  Green’s  function  G2  is  recognized  as  the  Fourier 
transform  of  the  Debye-Yukawa  potential 
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ou 

exp  (-  -p  R) 

G2  -  - 4^if -  <3*38> 

The  denominator  is  n exponentiated1'  to  put  the  integral  into  a 
form  to  be  compared  with  later  results.  In  this  case,  as  contrasted 
with  ,  the  denominator  is  not  singular  and  thus  a  real  exponential 
can  be  used. 


“  -  (d3k)  .v  _  .2.  u)2 

G2  =  J  dt  J  ^  e  ^  e"k  eXP^  “2 

2  J  J  (2rr)  3  c2 


(3.39) 


Again,  the  k-integrai  can  be  evaluated  by  completing  the  square 
in  the  exponential.  Hence 


GO 


G2  =  8TT3/2 


1  f  ,  jp  .»  ,  R  .  -3/2  .. 

3/2  J  exP(_  2  t)exP("  4t}t  dt 


(3.40) 


Since  is  known,  a  solution  for  the  integral  in  (3.40)  is 
easily  obtained.  For  C  >  0,  D  >  0 


00 


2  -  21 

Ax. 


-C‘t  -  4t  . -3/2  _  _£H~~CD 


dt  = 


D 

2 


(3.41) 


This  result,  Equation  (3.41)  ,  will  also  be  used  later, 


Evaluation  of  G_ 


The  determination  of  the  inverse  transform 


-1 


5 


Un) 3  ^  k| 

kz  '  a2  *  Y 


(3.42) 
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proceeds  along  lines  similar  to  those  of  and  G2*  It  shall  be 
shown  shortly  that  two  cases  must  be  considered,  one  of  which  is 
analogous  to  the  radiation  Green’s  function  ,  the  other  related 
to  the  Debye-Yukawa  potential,  G2 . 

Let  k T  denote  that  portion  of  the  propagation  vector,  k,  which 
is  transverse  to  the  magnetic  field  direction.  Similarly,  let  R^, 
denote  that  part  of  the  vector  R  which  is  perpendicular  to  the 
magnetic  field  direction. 

’’Exponentiating 11  the  denominator  in  (3.43)  and  arranging  terms, 
one  obtains 


ids  f 
(27 r)3  J 


dk 


z 


lk  Z  ik2s 
z  z 

e  e 


.2. 

d  kT  e 


ikTRT  , .  J"T 


exp(i—  s)exp(i^  s) 


Y 


(3.43) 


Completing  the  square  in  both  the  kT  and  kz  integrals  leads  to 


ids 
(2tt)  3 


i  77  \ 1//2  EX  i  ■ 

(u)  ltlexP(l 


CO 


s)  exp  J-(Z 


[t5< 


(3.44) 


2  2  —  2 
Let  IT  -  R^y  =  R  . 

If  R  is  positive,  the  integrand  in  Equation  (3.44)  is  exactly 

2  2  2 

the  same  as  the  integrand  in  Equation  (3.37)  ,  with  A  =  oo  /a  and 


Since  the  solution  of  the  integral  in  Equation  (3.37)  is  known, 
one  can  immediately  write  down  the  solution  of  (3.44) 


g3  =  =£-e 


i  «  R 
a 


Y  exp(i 


rtY) 


4ttR 


4tt^|z2- 


rt7 


(3.46) 
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This  is  analogous  to  the  radiation  Green’s  function  ,  except  that 

the  isotropic  radius  vector  R  has  been  replaced  by  the  quantity  R. 
-2 

If  R  is  negative  —  consider  the  s-integral  as  a  contour  integral 
along  the  positive  real  axis  in  the  complex  plane. 

If  one  integrates  along  the  positive  real  axis  from  0  to  Q 
(see  Figure  1) ,  and  then  along  a  circle  of  radius  rQ  in  the  first 
quadrant  of  the  complex  plane  to  P  and  n  takes  the  limit  as  ro*— *00, 
no  additional  contribution  is  made  to  the  original  s-integral. 


IM 


Figure  1 


For  along 

by  r  where 
o 

The  term  e 

as  r  00  0 
o 


the  curved  path  OP,  s  is  replaced  in  the  integrand 

0  <  8  <  \  . 

2  2  Z 

i  (w  /a  ) s  nQw  contains  a  damping  factor  for  all  0 


The  term 


iR2/4s 


=  exp 


iR2 

4r 


-10 


(3,47) 


goes  to  unity  for  all  values  of  0 ,  as  rQ— »oo.  Hence,  the  extra 
integral  gives  no  contribution.  Note  that  there  are  no  singularities 
within  the  quadrant  to  prevent  use  from  deforming  the  path  of  inte¬ 
gration  to  the  positive  imaginary  axis. 

—2  2  2 

Hence,  if  R  <0,  i.e.,  R^y  >  Z  ,  there  is  no  contribution  to 
the  s-integral  along  OP. 
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If,  instead  of  integrating  along  OQ  and  then  QP  the  integration 
is  performed  directly  along  the  imaginary  axis  from  O  to  P,  s  is 
replaced  by  it  in  the  integrand  in  Equation  (3.44) , 


oo 


G3  = 


^-372  J  exp(-  %  t)exp(+  f^)t-3/2 


(3.48) 


—2 

But,  for  R  <0,  this  integral  is  exactly  the  same  as  the  integral 
expression  of  the  Debye-Yukawa  potential  as  given  by  Equation  (3.41) 
with  C2  =  oo2/a2  and  D2  =  -R2  .  Therefore 


(3.49) 


—2  —2 

What  is  the  physical  significance  of  the  conditions  R  >  0  or  R  <0? 

—2  2  2 

R  =  0,  i.e.,  Z  =  Rty  defines  a  cone  in  a  three-dimensional  space < 
If  for  simplicity  it  is  assumed  that  the  source  is  located  at  the  origii 
(r^  =  0;  z'  =  0)  ,  then  this  cone  is  further  defined  by 


z 


2 


(3.50) 


2  2  2 

where  r  =  z  +  rT  is  the  distance  to  the  field  point.  Inside  the 
2  A  2 

cone,  where  z  >  r_v,  the  fields  as  given  by  are  radiation-like 

T  ?  ?  -1/2  ** 

expressions  which  fall  off  as  (z^-r^y)  , similar  to  the  1/r  dependence 

of  a  typical  radiation  field.  On  the  surface  of  the  cone  the  field 

2  2 

is  singular  and  infinite.  Outside  of  the  cone,  where  rTy  >  z  ,  the 
fields,  as  given  by  ,  are  damped  exponentially  and  do  not  corres¬ 
pond  to  radiation.  Therefore,  if  there  is  any  radiation  associated 

with  the  Green’s  function  G-,  it  is  contained  within  the  cone  defined 
2  2 

by  z  =  ripY . 
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Evaluation  of 

The  denominator  in  the  integral  expression  for  is  a  product 
of  the  terms  appearing  in  the  expressions  for  G-^  and  G^  respectively. 


•% 

v 


(d3k) 


(2tt)3  -2  a/* 

K  '  2  k^ 

a  z 


2/  2 
a)  /a 


Y 


(3.51) 


Each  term  in  the  denominator  can  be  "exponentiated"  and  when  like 
terms  are  grouped 


(27 r)  3 


oo  oo 

J  at  J  as 

o  o 


exp  ( ikzZ)  exP  ( -ik 


2 

z 


(s+t) 


J  d2kT  exp(ikTRT)  exp  |-ik2  (s  -  ~)J  exp  i  ^(s+t) 


(3.52) 


Again,  the  kz  and  kT  integrals  can  be  evaluated  by  completing 
the  squares  in  the  exponents 


-w0  oo  oo 

o  o 


ds 


[iR2 

r  .  -2  1 

r  2  i 

exp 

T 

exp 

l  z 

exp 

i  ~2  ( s+ 1) 

4 ( s-— ) 

4  (s+t) 

L_ _ xU 

(s  -  p  (s  +  t) 


1/2 


(3.53) 


Now  consider  a  double  change  of  variable. 

Let 

s+t  =  u  (3 ,54) 
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The  Jacobian  of  this  tranformation  is 


J(s,t,U,v)  =  ^  + 


(3.56) 


In  determining  the  new  limits  of  integration  ,  it  is  noted  that 
the  infinite  quarter  plane,  bounded  by  the  positive  s  and  t  axes,  must 
be  covered  by  integrating  along  a  pair  of  oblique  axes,  u  and  v.  This 
is  accomplished  by  integrating  over  v  from  the  line  s  =  0,  which 
corresponds  to  v  =  -u/y ,  to  the  line  t  =  0  which  corresponds  to  v  -  u, 
and  then  along  the  u-axis  from  zero  to  infinity 


G4  = 


-  J du  u_1/2  exp(i  %  u)exp(i  J  dv 

8rr  a  y 

o  -u/y 


exp(i  j~) 


1+y 


(3.57) 


Consider  the  special  case  of  +0.  This  corresponds  to  the 
source  being  located  at  the  origin  and  an  observer  looking  along 
the  magnetic  field  direction.  In  the  v-integral ,  for  finite  the 
singularity  at  v  =  0  is  damped  out  symmetrically  for  |v|  <  R^/4  by 
the  presence  of  the  exponential.  This  defines  the  Cauchy  principle 
value  of  the  integral.  In  the  limiting  case  of  — ►  0 ,  in  order  for 
the  integral  to  converge,  one  must  bear  in  mind  that  we  are  concerned 
with  its  principal  value 


*  J  f  -in  M  J  “  i"Y  • 

-u/y  -u/y 


(3.58) 


Then 


g4(rt  =  0)  = 


,  i-3/2  ? 

1+Y  Y  „_3/2  J  du 


,-1/2 


(3.59) 
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Making  the  substitution  u  =  1/x,  the  remaining  integral  takes 
the  form 


f  j  -3/2 
dx  x 

t 


o 


2  2 
exp  ( i  —  x)  exp  ( i  — 


2  ' 
a  x 


This  integral  has  been  encountered  previously  in  the  calculation 
of  G i  and  and  its  solution  is  given  by  Equation  (3.37) . 

For  the  present  case 


A 


2 


Therefore,  from  Equation  (3.37) 


(3.60) 


oo 

r 


j 


dx  x 


-3/2 


o 


2  2 
exp  ( i  x)  exp  ( i  ^T") 
ax 


exp(i  -  Z) 
a 


(3.61) 


Then 

.  exp(i  j  Z) 

01  -  eTI0?  (3'62> 

a 

This  expression  differs  from  the  previous  radiation-like  Green's 

functions  that  have  been  calculated.  G^ (R^  =  0)  represents  a  plane 

wave  propagating  along  the  magnetic  field  axis.  There  is  no  decrease 

in  intensity  with  distance  away  from  the  origin.  This  result  is  not 

completely  unexpected.  represents  an  isotropic  mode  of  radiation 

and  is  analogous  to  the  isotropic  mode  derived  in  Chapter  I  for  the 

—2 

case  of  zero  sound  velocity  (Equation  (1.13))  .  G^  9  for  R  >0, 
represents  radiation  filling  a  cone  of  small  angle  about  the  mag¬ 
netic  field  line  and  is  related  to  the  one-dimensional  Alfven 
mode  of  Chapter  I  (1.14)  which  is  an  idealization  of  the  actual  physi¬ 
cal  situation.  G^  is  a  measure  of  the  radiation  resulting  from  the 
interaction  of  these  two  aforementioned  modes. 
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For  the  case  of  Rj  /  0,  an  interesting  integral  representation 
of  G^  can  be  obtained. 

2 

If  the  partial  derivatives  of  is  taken  with  respect  to  R^, 
Equation  (3.56)  becomes 

_2 


oo 


3R^  4  32f3//2  q 


172  J  du  U 
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(3.63) 


The  integral  over  v  can  now  be  performed  since  the  differential 
of  the  exponent  appears  in  the  integrand. 
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(3.64) 


If  again  the  substition  u  =  1/x  is  made,  both  expressions  in 
the  integrand  can  be  put  into  a  form  similar  to  that  of  the  integral 
in  Equation  (3.37)  ,  our  known  identity. 
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(3.65) 
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- 2 “co  1+y  [exP(i  a  ^Z2  +  rJ)  -  exp(i  ^  ^Z2  -  R^y)l  (3.66) 

8ttRt  —  L  J 


The  Green’s  function  is  then  obtained  by  integrating  over  all 


values  of  R^. 
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(3.67) 


where  the  term  outside  of  the  integral  is  the  value  of  the  Green's 
function  for  RT  =  0 ,  which  has  been  previously  calculated  (Equation  (3,62) 

2  2 

It  is  seen  from  this  integral  expression  for  G^ ,  that  for  RTy  <  Z  , 

the  Green’s  function  is  radiation-like ,  and  that  it  is  partially  the 

2  2 

sum  of  contributions  over  cones  described  by  Z  =  R^y  (for  fixed  Z)  . 

There  is  a  contribution  at  =  0  which  represents  a  plane  wave  solu¬ 
tion  on  the  magnetic  field  axis. 


Evaluation  of  G^ 

The  denominator  in  the  integral  expression  for  G^  is  a  product  of 

the  terms  appearing  in  the  expressions  for  G2  and  G^ •  Physically,  this 

Green's  function  represents  the  interaction  of  the  directional  mode 

2  2 

represented  by  inside  of  the  cone  Z  =  RRy  with  the  attenuated  mode 
represented  by  the  Debye-Yukawa  potential.  One  expects  the  product 
of  the  two  to  be  a  directional  damped  oscillation 
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Each  term  in  the  denominator  can  be  "exponentiated. M  An  im¬ 
aginary  exponential  must  be  used  for  the  term  containing  the  singularity. 
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(3.70) 


Previously,  we  have  obtained  exactly  for  the  case  of  R^  =  0. 

In  the  present  example,  we  shall  consider  the  case  when  Z  =  0  and  show 
that  Gj-  (Z  =  0)  can  be  determined  in  an  analogous  manner. 


Both  s  and  t  are  real  variables.  If  the  upper  limit  of  the  t- 
integration  is  extended  to  cover  the  circle  in  the  first  quadrant 

iQ 

by  lim  r  e  ,  it  is  easily  seen  that  there  is  no  contribution  to 
rD— ►co  o  ’ 

the  t-integral  along  the  extra  portion  of  the  path.  For,  if  t— ►t 

2  2  t 

+  ix ,  t  >  0 ,  the  term  exp  i  oo  /c  t  contributes  e  P'  which 

P  2  t 

approaches  zero  as  x— ►oo,and  the  term  exp (R,p/4  ( s-i— )  )  is  bounded 
in  this  limit.  Then,  as  was  shown  previously  in  the  calculation  of 
G^ ,  the  contour  can  be  deformed  and  the  integration  performed  along 
the  positive  imaginary  axis. 


It  should  be  noted  that  the  contour  cannot  be  deformed  in  this 
manner  when  Z  ^  0. 

iZ2 

For  then,  the  term  e  7-7 — — becomes  an  isolated  essential 
singularity  at  a  point  on  the  path  of  integration  —  namely  where 
t  =  is,  and  gives  rise  to  an  infinite  contribution  to  the  integral. 
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For  the  Z  =  0  case,  this  term  vanishes  and  the  above  difficulty 
is  circumvented. 

Then  t  can  be  replaced  in  the  integrand  of  Equation  (3.70)  by  iT< 

exp 
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(3.71) 


One  now  proceeds  along  lines  similar  to  those  used  in  the  evalu¬ 
ation  of  .  Introduce  the  double  change  of  variable 


S  — T  =  U 


s+—  =  v 

Y 


(3.72) 

(3.73) 


The  Jacobian  of  the  transformation  is  again 
The  new  limits  of  integration  are: 
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The  u-integral  has  an  integrable  singularity  and  can  be  per¬ 
formed  directly 
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Gr  = 


47 T 


3/ 


JL 


1+i 


i  »Jr 


OO  2  2 

J  dv  V"1/2  exp(-  v)exp(-  ^) 


(3.75) 


Jf  one  now  makes  the  substitution  v  =  ^  ,  the  resulting  inte¬ 
grand  is  exactly  the  same  as  that  in  Equation  (3.41)  ,  the  integral 
expression  for  the  Debye-Yukawa  potential .  Therefore 
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and  from  Equation  (3.41) 
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Thus,  G^(Z=0)  represents  an  ’’attenuated  plane  wave"  traveling  in  the 

direction  transverse  to  the  magnetic  field  lines.  This  result  is  not 

unexpected,  since  it  represents  the  "interference"  between  the  damped 

Debye-Yukawa  potential  and  the  exponentially  attenuated  oscillation 

2  2 

given  by  outside  of  the  cone  Z  =  R^y  . 

As  discussed  previously,  the  deformation  of  the  contour  is  not 
valid  for  Z  /  0,  and  one  cannot  proceed  in  the  manner  described  above. 
Therefore,  in  order  to  get  an  expression  for  G5(Z  /  0)  ,  consider  the 
following  scheme 
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Now  decompose  the  product  in  the  denominator  into  a  sum  of  terms  9  by 
the  method  of  partial  fractions 
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and  making  use  of  the  Fourier  transform  relationship  between  k  and 
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from  the  integral  expressions  for  the  transforms  of  G2  and  G3  . 
Therefore , 


GO 

exp(-  -£■  R) 
47TR 


exp  ( i  —  R) 
4frR 


(3.83) 


a 

The  next  task  is  to  find  the  Green’s  function  for  the  operator  - — 2  * 
g (Z,Zf )  ,  where 


g(Z,Z')  =  6  (Z-Z ' )  .  (3.84) 

dZ 


Then 


G5  =  -  J  g(Z,Z')  f (X,Y,Z ' )  dZ' 


(3.85) 


From  the  physics  of  the  problem,  we  know  that 

g(Z,Z')  =  g(Z-Z')  (3.86) 

2  2 

g(Z-Z’)  must  satisfy  5  /9Z  g(Z-Z’)  =  0  everywhere  except  at  Z  =  Z’ 

and  the  change  in  the  first  derivative  must  be  +1  as  one  passes 
through  the  point  Z  =  Z’  from  above.  This  Green’s  function  is  easily 
obtained.  Consider  the  function 


1 
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|z  -  Z'| 


(3.87) 


Its  slope  changes  by  +1  as  the  point  Z  =  Z'  is  traversed 
and  its  second  derivative  is  zero  everywhere  for  Z  ^  Z'  . 
this  function  satisfies  92/9Z2  |z  -  Z'|  =  b( Z  -  Z').  We 
that 


from  above, 
Therefore, 
shall  see 


g(Z  -  Z')  =  J  |z  -  Z’ 


(3.88) 
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without  any  further  homogeneous  contribution  gives  the  previously 
calculated  Z  =  0  result.  Then 

Z 

G5(Z)  =  "  2  J  (Z  "  Z,)  [g2{Z,)  +  y  G3 (Z)]  dZ' 
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-  2  tJt  J  <z'  -  a  [<VZ'>  +  r  G3<z,)]dz'  • 

z 

(3.89) 

G^,  it  must  be  remembered,  has  two  functional  forms,  one  denot¬ 
ing  propagation  and  the  other  damping.  The  appropriate  form  of 
must  be  used  in  the  integrals  in  Equation  (3.89) ,  depending  on 
whether  the  value  of  |z|  in  the  limits  is  greater  than  or  less  than 
RT  Therefore,  Equation  (3.89)  can  be  explicitly  written 

exp  (i  ^  ^Z2  -  R^y)  exp  (-  Z2  +  R2) 
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for  | Z  |  <  Rt  NP  ( 3  *  91) 


In  writing  down  Equations  (3,90)  and  (3.91)  we  have  made  use  of  the 
boundary  condition  that  all  the  exponentials  in  the  integrand  approach 
zero  at  Z1  =  t  oo  .  This  corresponds  to  the  conditions  of  the  actual 
physical  problem  where  even  the  radiation  modes  are  not  infinite  but 
go  to  zero  at  a  sufficiently  large  distance  from  the  source.  It  is 
seen  from  Equation  (3.90)  that  the  first  term  in  that  expression  does 
represent  radiation.  Thus,  the  interaction  of  the  radiation  term 


-  iy 


•  00  I 

e1 

I  z  -rtY 

4  TT 


with  the  damped  mode  described  by  the  Yukawa  potential ,  leads  to  a 

2  2 

radiation-like  term  inside  the  cone  Z  =  Rty  ,  in  G5  . 

The  reader  can  easily  verify  that  Equation  (3.91)  reduces  to  the 
result  of  our  previous  calculation  for  Z  =  0,  [c.f.  Eq.  (3.83)J  , 
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The  radiation  problem  has  been  so  set  up,  that  now  with  the 
elements  of  the  matrix  A  determined,  one  must  merely  evaluate  a 
single  current  source  j (r,t)  .  In  Chapter  VI,  the  radiation  fields 
arising  from  several  different  types  of  current  excitations  will  be 
calculated  and  discussed  regarding  their  physical  significance. 
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CHAPTER  IV 

RADIATION  FROM  UNIFORMLY  MOVING  CHARGES 

Among  the  processes  by  which  charged  particles  are  slowed  down  in 
passing  through  solid  and  liquid  material  media,  Cerenkov  radiation  is 
one  of  the  most  efficient.  Such  a  state  of  affairs  might  be  true  also 
for  charged  particles  passing  through  ionized  magnetized  'fluids  1  ,  or 
plasmas,  for  which  the  conditions  resemble,  in  some  respects,  the 
conditions  needed  for  Cerenkov  radiation.  It  will  be  the  purpose  of 
the  next  two  chapters  to  reveal  the  analogy  in  more  detail ,  and  to 
estimate  the  quantitative  importance  of  the  results.  This  work  bears 
on  the  dynamics  of  the  ionized  zones  surrounding  the  earth,  but  is  too 
restricted  to  be  applied  directly  to  that  problem.  The  charged  particles 
which  originate  on  the  sun,  and  reach  the  earth's  magnetosphere,  arrive 
in  bursts  of  plasma  of  densities  perhaps  comparable  to  that  of  the 
ambient  plasma,  rather  than  as  independent  particles.  What  such 
streams  do  upon  reaching  the  vicinity  of  the  earth  must  be  investigated 
by  other  means. 

A  moving  charged  particle  does  not  radiate  in  traversing  empty 
space  unless  it  is  accelerated.  The  latter  possibility  is  realized  if 
the  particle  has  a  component  of  velocity  across  the  lines  of  the  mag¬ 
netic  field.  But  such  "cyclotron  radiation"  is  usually  important  only 
for  extreme  relativistic  motions,  because  the  accelerations  involved 
are  not  very  great.  The  most  efficient  cases  of  the  radiation  result¬ 
ing  from  single  accelerated  particles  occur  when  fast  particles  are 
brought  to  rest  by  dense  solids  ("bremsstrahling" ) :  the  fast  particle 
is  deflected  strongly  by  the  Coulomb  fields  of  the  heavy  nuclei,  and 
radiates  considerably  in  virtue  of  this  acceleration.  But  for  gases, 
and  plasmas,  this  effect  is  also  to  be  considered  an  unimportant 
mechanism  for  the  "stopping"  of  the  particle. 

One  must  turn  to  "collective  behavior"  to  find  an  important 
energy-loss  mechanism.  In  crossing  the  plasma,  the  charged  particle 
accelerates  the  ions  in  its  neighborhood.  While  no  one  of  the 
individual  ions  is  accelerated  greatly,  there  are,  however,  many 
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which  are  set  into  motion ,  and  with  definite  phase  relations,  so  that 
the  actual  fields  can  be  quite  considerable.  This  radiation  of  the 
induced  currents  set  up  in  a  solid  transparent  substance  was  first 
observed  (accidentally) by  Cerenkov  in  1934,  Its  existence  had  not 
been  suspected  even  though  closely  related  calculations  had  previously 
been  made  of  the  field  around  a  charge  moving  faster  than  the  velocity 
of  light.1 The  field  lies  entirely  in  a  cone  behind  the  particle,  and 
is  singular  on  the  surface  of  the  cone ,  The  Poynting  vector  at  the 
surface  is  normal  to  the  cone  and  outwards,  but  the  rate  of  radiation 
could  not  be  calculated.  In  any  case,  the  theory  of  relativity  (that 
followed  the  calculation)  eliminated  the  possibility  of  charges  moving 
faster  than  light,  Franck  and  Tamm  in  1937  explained  Cerenkov ’s 
discovery  by  treating  a  charge  moving  faster  than  the  velocity  of  light 
in  the  medium.  All  of  the  details  of  the  medium  were  compressed  into 
one  function  —  the  frequency  dependent  phase  velocity  of  light  (or 
the  index  of  refraction).  Their  treatment  of  the  problem  is,  there¬ 
fore,  often  referred  to  as  the  ’’macroscopic  approach,”  since  the  only 
information  needed  about  the  material  is  its  index  of  refraction  for 
light  of  different  frequencies,  a  bulk  property  of  the  medium.  It  is 
a  method  of  extreme  generality,  however,  being  restricted  primarily  to 
a  linear  approximation  to  the  response  of  the  medium  to  excitations. 

But  within  the  domain  of  linearity,  one  can  treat  nearly  all  important 
problems  of  the  interaction  of  charged  particles  and  radiation  with 
neutral  and  ionized  substances 

(12) 

As  is  well  known  from  the  work  of  Alfven  ,  an  incompressible 
conducting  fluid  in  a  DC  magnetic  field  has  a  one-dimensional,  low- 
frequency  mode  of  propagation  of  electromagnetic  energy  —  a  trans¬ 
verse  field  propagating  along  the  lines  of  force  of  the  magnetic 
field  at  a  velocity  a  which  can  be  considerably  less  than  the  velocity 
of  light  in  vacuum  (c) .  This  Alfven  (hydromagnetic)  mode  is  not 
suitable  for  Cerenkov  radiation,  which  is  a  charateristic  of  an 
’’isotropic”  mode  —  that  is,  one  which  propagates  equally  well  in 
all  directions.  The  reason  for  this  is  that  the  instantaneous 
Cerenkov  radiation  lies  very  nearly  in  the  surface  of  a  cone  of 
finite  angle.  At  the  very  best,  the  direction  of  the  magnetic  field 
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lines  (and  therefore  of  the  propagation  vector  of  the  Alfven  wave) could 
be  parallel  to  only  one  generating  line  of  this  cone?  tT*re  can  be  no 
finite  density  of  radiation  in  one  line. 

For  a  compressible  fluid,  however,  there  exists  in  addition  to  the 
Alfven  mode  an  n isotropic”  mode  whose  phase  velocity  is  the  same  as  the 
Alfven  velocity.  Since  it  depends  on  the  presence  of  compressibility, 
it  is  often  referred  to  as  a  "hydromagnetic  sound  wave.”  It  is  this 
mode  that  will  be  investigated  for  Cerenkov-like  radiation.  It  is 
well  defined  so  long  as  the  actual  sound  velocity  in  the  plasma,  s,  is 
much  less  than  a?  a  larger  value  of  s  corresponds  to  a  more  rigid, 
incompressible  medium.  For  the  zero  temperature  gas  we  have  in  mind, 
s  is  near  zero.  The  difference  between  the  ordinary  Cerenkov  effect 
and  the  hydromagnetic  example  lies  in  the  interference  of  the  one¬ 
dimensional  mode  with  the  " isotropic”  mode. 

ORDINARY  CERENKOV  RADIATION 

We  have  discovered  that  the  tensor  X(koo)  characterizes  the  electro¬ 
dynamical  properties  of  a  " slightly  disturbed"  medium.  For  most  solids 
and  liquids  (in  the  frequency  regions  for  which  they  are  transparent) 
X(koo)  can  be  adequately  represented  by  a  positive,  diagonal,  isotropic 
tensor  which  depends  only  on  oo: 


X(koo)  - 


XQ  (oo)l 


(4.1) 


2 

Consequently ,  the  propagation  denominator  k 
2 

k2  -  ^(1+  X(oo))  . 
c 

The  quantity 


a  (to)  =i 


1  +  x  (co) 
o 


<  C 


becomes 

c 


(4.2) 


is  the  "phase  velocity"  for  light  at  frequency  oo.  In  terms  of  the 
notation  of  Chapter  II,  we  can  show  the  following  in  this  case: 
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We  have  used  the  conservation  law; 
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(4.6) 
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(4.8) 


The  radiation  of  the  external  sources  may  be  found  from  the  Poynting 
vector: 


S (rt)  =  —  E (rt )  x  B(rt)  (4.9) 

|XO 

Let  the  currents  be  localized  in  space  and  act  only  for  a  finite  time. 
Consider  a  large  sphere  drawn  with  the  radiating  currents  near  its 
center,  which  may  also  be  taken  to  be  the  origin  of  a  system  of  spheri¬ 
cal  coordinates:  0,  An  element  of  surface  on  the  sphere  will  be 

o 

dA  -  r  sin  0  d  0  d  Its  normal  will  be  along  the  radius  and  out¬ 

wards:  n  =  y-y  .  The  energy  flux  at  time  t  through  dA  per  unit  time 
will  be 
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The  total  energy  passing  through  the  fixed  element  (dA)  will  be 

oo 

dAn  .  j  dt  S (rt)  (4.11) 

LaJ 


Measuring  instruments  (e.g.,  the  eye)  generally  respond  to  the  frequency 
components  of  the  electric  and  magnetic  fields;  their  contribution  to 
(4.11)  may  be  found  by  introducing  the  transforms  for  E(rt)  and  3(rt) 
making  use  of  the  Fourier  integral  theorem: 

oo 

J  dt  e^^  e-0^  -  2tt  b  (go  +  rof)  (4.12) 

_  oo 

and  noting  that  for  the  transform  of  a  real  function,  f (t) ,  we  must 
have  a  condition  on  the  complex  conjugates: 

f  *  (oo)  =  f  (-00)  (4.13) 

We  learn,  then,  that  (4.11)  becomes 

oo 

dAn  •  j  2m  ~~^o  Rej^E(roo)  x  B*  (rto)J  (4.14) 

o  ^  ° 

j  ^ 

When  we  realize  that  vjp  =  dv  ,  where  v  =  (£>/2v  is  the  frequency  of  the 
observed  light,  we  can  identify 

dl(v)  =  ■  dA  Re  ft  .  Jjs(rto)  x  B*(rv)j  (4.15) 

with  the  frequency  distribution  of  the  energy  flux  of  the  electro¬ 
magnetic  fields  through  the  surface  element  on  the  sphere.  If  dl 
approaches  a  finite  non-zero  limit  as  the  radius  of  the  sphere  is 
increased,  the  flux  may  be  called  radiation,  since  it  is  lost  energy. 
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A  particle  moving  uniformly  will  radiate  continuously.  To  use  it 
in  the  above  scheme  requires  that  we  formally  treat  it  as  being  on  for 
a  time  T.  We  take  the  path  as  linear,  with  velocity  vector  v  ,  passing 
through  the  origin  at  time  t  =  0: 


pext(rt)  =  q&  ^ (r-vt) 

jext(rt)  -  qv  b  ^ (r-vt)  =  v  pext(rt) 


(4.16) 


We  shall  return  to  (4.6)  and  (4.7)  and  rewrite  them  in  the  rco 
representation,  by  means  of  a  convolution  in  which  the  singular  de¬ 
nominator  is  isolated. 


Using  Equation  (4.16),  we  obtain 

E(rco)  =  ipo  (crs  +  a^i^)J  d^r  'G^  (r-r  ')  pext  (r  rco) 


(4.17) 


B(rto) 


-  y-r 


V  X  V 


% 


d^r  ?G^  (r-r  f )  pext  (r  Too) 


(4.18) 


where 


G1(r-r  f) 


i 


e 

4t r 


a 


r-r 


r-r  »’[ 


(4.19) 


As  is  usual  in  radiation  problems,  we  take  the  asymptotic  form  of 
.  The  source  point  rf  is  supposed  to  be  near  the  origin.  For 
large  r, 


;-r  ' 


A 

n.  r 


? 


Therefore , 


oo 


,3llr-r')~  4F 


A 

n.r 


i 


(4.20) 


(4.21) 
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In  the  evaluation  of  the  fields  E  and  B  o^  Equations  (4.17,  18)  the 

i  —  r 

gradients  are  allowed  to  act  only  on  e  ,  since  otherwise  the 

powers  of  r  produced  would  not  be  appropriate  to  a  finite  contribution 
to  the  radiation.  Hence,  in  Equations  (4.17,  18)  we  may  replace 


V  — »  i  ®  V|r| 

a 


.  0)  A 

=  1  £  n 


(4.22) 


To  evaluate  the  r  '  integral  we  write 


oo 


ext 


(r'oo)  -  q  |  dt  vt)e^a:>t 


-  oo 


so  that 


.  0)  A 

Jd3r'  e  a  ’  pext(r'oo)  =  q  J  dt  e 


00  .  .  ,  CO  A  . 

it  (oo  -  —  n.  v) 


_  00 


(4.23) 


(4.24) 


=  27rq5  (oo  -  -  fi.v) 


Combining  these  results  leads  to  the  radiation  fields 


E(roo)  =  iu  qoa(fi  -  ~) 


v,  e 


.  oo 

1  a  r 


2tt5 


hi- 


(4.25) 


and 


.  oo 

1  a  r 


B(eoj)  =  i|ioq“  v  x  n  e4Tr  - —  2rr&|co(l-  n.^)]  (4.26) 


which  contribute  to  the  power  flux 
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Re  —  E  ( rco) 


x  B*(ru>) 


=  “  ft(v2_  a2)  b  co(l-  2tt6  (0) 

37?  |  r  |  2a  a 


(4.27) 


The  infinite  quantity  2irb  (oo— oa )  =  2tt?)  (0)  .  arose  from  the  square  of  the 
6-  function  in  E,  which  in  turn  arose  from  the  time-integral  in  Equation 
(4.23).  The  use  of  a  source  which  is  on  for  all  times  must  lead  to  an 
infinite  amount  of  power  radiated.  Only  the  rate  of  radiation  is 
physically  significant  in  such  a  case.  Formally,  we  can  resolve  the 
difficulty  if  we  let  the  source  be  on  only  in  the  time  interval  -T/2 
to  T/2,  where  T  is  chosen  arbitrarily  large.  Then  the  integral  in 
Equation  (4.24)  is  still  2nb  (0)  to  sufficient  accuracy,  but  the  value 
of  2nb  (0)  at  its  singularity  is  only  the  overall  time  T.  With  this 
understanding 

2Tf6  (0)  =  T  (4.28) 


The  rate  of  radiation  becomes 


dR  ~  2tt  00  4^a  (v2-  a2)?)(l  -  cos  0)  sin  OdOdpf 

(4.29) 

in  which  we  have  chosen  the  vector  v  along  the  polar  axis  of  spherical 
coordinates,  and  have  used 

5  (ax)  -  y~-y  b  (x)  (4.30) 


Using  it  again, 

jx  2  2 

dR  =  q  dv  vVj^  6  (cos  0  -  y^p)  d(cos  0)  dd 

(4.31) 

So,  unless  v  is  greater  than  a  (the  phase  velocity  in  the  medium  at 
frequency  v  =  there  is  no  radiation  at  all.  The  total  radiation 

per  unit  time  is 
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2  2 

R  =  n^oq2  J  dv  vV  TvVVj  (4.32) 

the  integral  being  over  those  frequencies  for  which  v  >  a  (v).  All  of 

the  radiation  (at  frequency  v)  lies  in  the  surface  of  a  cone  of  angle 

0  =  cos"*1  (y—y)  around  the  direction  of  the  particle,  A  transparent 

substance  is  capable  of  supporting  visible  Cerenkov  light  by  virtue 

of  its  being  transparent.  The  efficiency  of  the  process  is  a  consequence 

2 

of  the  high  values  of  v  in  the  optical  region.  As  we  shall  see,  the 
corresponding  process  in  the  plasma  behaves  differently  at  low  frequence 


53 


S-2023-1 


CHAPTER  V 

THE  MAGNETO  GAS  DYNAMIC  ANALOGUE  OP  CERENKOV  RADIATION 

Let  us  now  proceed  to  determine  the  total  radiation  produced  by  a 
single  point  charge  moving  with  uniform  velocity  v  through  a  partially 
ionized , inf inite ,  anistropic  plasma,  in  the  direction  of  an  external 
magnetic  field.  The  plasma  under  consideration  is  again  the  upper 
region  of  the  earth’s  ionosphere.  The  point-charge  excitation  might 
be  one  of  the  countless  hydrogen  ions  which  are  emitted  near  the  sun? 
it  might  traverse  inter-planetary  space  in  corpuscular  streams,  and 
arrive  at  the  terrestrial  exosphere  with  speeds  of  100  km/sec  to  1,000 
km/sec.  in  this  region.  Therefore,  one  might  expect  the  uniformly 
moving  charge  to  radiate  in  a  manner  analogous  to  the  Cerenkov  effect. 

In  the  normal  Cerenkov  case,  the  electric  field  produced  by  the 
particle  lies  within  a  cone  given  by  cos  0  =  a/v  around  the  direction 
of  the  particle,  and  all  of  the  radiation  lies  in  the  surface  of  this 
cone.  As  was  shown  previously  in  the  calculation  of  the  Green’s 
function  ,  the  plasma,  for  a  certain  mode  of  oscillation,  behaves 
like  a  conical  wave-guide  with  propagation  allowed  within  a  cone 
described  by 

,  ^  1  _  a)  c  y/  , 

tan  Q  =  -  =  —  -  <<  1. 

i -  co  a 

V Y  p 

This  cone,  which  we  shall  call  the  Aliven  cone,  arises  because  of  the 
basic  anisotropy  of  the  medium  and  is  a  pure  magnetohydrodynamic  effect. 
In  general,  for  v  >  a,  the  Cerenkov  cone  lies  outside  of  the  Alfven  cone 
when  the  two  vertices  coincide.  However,  when  the  particle  velocity 
becomes  very  near  the  phase  velocity,  the  two  cones  will  coincide  and 
one  might  expect  to  find  another  mode  of  Cerenkov  radiation  which 
propagates  in  the  conical  waveguide ,  above  and  beyond  the  analogy  to 
the  normal  Cerenkov  effect. 
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CALCULATION  OF  THE  ELECTRIC  FIELD 

Both  collisional  and  thermal  damping  are  negligible  for  ultra-low 
frequency  magnetohydrodynamic  waves  propagating  along  the  lines  of  the 
earth ’s  magnetic  field,  for  the  collision  frequencies  and  temperatures 
encountered  in  the  exosphere*  The  latter  case  will  be  discussed  in 
detail  in  the  next  chapter.  Therefore,  we  can  use  the  collisionless, 
zero-temperature  Green’s  functions,  introduced  in  Chapter  III,  to 
determine  the  radiation  produced  by  a  uniformly  moving  point  charge 
as  described  above. 

The  physically  observable  quantity,  the  frequency  components  of 
the  electric  field,  can  be  written  as  a  convolution  of  the  external 
currents  and  the  various  Green’s  functions 


E(ru>)  =  X(i qC0  J  d3r’  A(r-r’,oo)  .  jext(r’,oo)  (5.1) 

ext 

j  is  the  current  produced  by  the  uniformly  moving  point  charge. 
Such  a  point  charge,  moving  in  the  z-direction  with  velocity  v  is 
described  by 

jext(rt)  =  gvz5^(r-vt)  =  qvz6 (x)  b  (y)  6  (z-vt) 

(5.2) 

ext 

Then  j  (rtu)  is  obtained  by  use  of  the  Fourier  integral  theorem 


j  (roo)  =  J  eia>t  j  (rt) dt  =  qfc  b  (x)  b  (y)  exp  (i  ^  z) 

(5.3) 


Thus  from  Equation  (5.1) 


E^rto)  =  ip,Qcoq  Jdz  •  A^(x  ,y  ,z-z  1  ,co)  exp(i  ^  z  ') 
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where 


and 


xz 


00 


2  axaz  3 


yz 


00 


dydz  G3 


zz 


G2  - 


00 


9z 


2  G3 


(1  +  ~) 

Y 


^G4 


Y 


60  2  2 

_e_  §_ 

,  2  2 

00  c 


(5.5) 


(5.6) 


(5.7) 


(5.8) 


G^,  ,  and  G^  are  the  zero-temperature  collisionless  Green’s 

functions  discussed  at  length  in  Chapter  III. 

G2(roo)  is  the  exponentially  damped  Debye-Yukawa  potential, 


exp  ( -  -^  I  r  | ) 
4tt  frl 


which  can  be  neglected  insofar  as  radiation  effects  are  concerned. 
Normally,  this  term  would  be  important  for  the  case  of  the  static 
charge,  but  at  zero-temperature,  the  static  electric  field  is  also 
zero. 


Making  use  of  the  result 

i  - 

J  dz  *  f  (z-z  ’)e  v 


z  ’ 


i  —  z  ’ 

-  J  dz  '  e  V  f  (z-z  ') 

i  -  z  ' 

i  J  dz  'f  (z-z  ')e  v 


(5.9) 
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[with  the  integrated  terms  vanishing  because  we  demand  that  f(z-z')  -  0 
at  z  1  -  l.00]  ,  the  electric  field  can  be  rewritten 


E(r,co)  =  ip,  coq  (-  i  V  j  dz  '  G3(xyz-z')e  v 

co 


.  CO 

1  —  z  ' 


(5.10) 


+  iu  coq  -  *  1  ^  z  J  dz '  G.(xyz-z')e  v 


i  z. 


2 

~  + 


E^£ 


(5.11) 


where 


(5.12) 


(4) 

is  that  part  of  the  field  due  to  ,  and  E  is  that  part  of  the  field 

due  to  G. .  Furthermore,  since  iooB  =  Vxfi  from  Maxwell's  equations,  and 

4  (4  j 

the  curl  of  a  gradient  is  zero,  B  is  determined  entirely  from  E  . 

(4) 

Our  first  step  is  therefore  to  calculate  E  ,  The  Green's  function 
G^  can  be  written  as  a  convolution  of  and  G^ . 


G4  (xyz-z  ') 


dr  '  ' 


G1  (r-r  '  0  G3  (r  ’  f-r  ')dr  '  ' 


(5.13) 


since  in  wave-number-frequency  space,  is  just  the  product  of  the 
other  two,  by  definition. 


We  shall  use  to  best  advantage,  an  integral  representation  for 
G^  that  was  developed  earlier,  rather  than  its  explicit  form  which 
denotes  propagation  inside  and  damping  outside  of  a  cone  of  angle 


tan 


-1  1 


a/T~ 
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G  (x  '  'y"  ,z  '  '-z  ') 


i-1/2 
3/2 


oo 


I 


87 T  '  OS 


ds 

-J/2  exp  1 


2  (z  '  '-z  ')  2— r  '  '2y 

— —  =5+  - — - 

2  ^  4s 

a 


then 


(5.14) 


.1/2  , 

, ,  4)  1  a)3  Y  ±  1 

'  .,3/2  2  Y 

2  77  '  v 


oo 


o  s 


ds 

3/2 


dz  '  dr  1  ' 


( z  '  '  ~  z  ' ) 

2  “j1  4s 

.a 


x  exp  (i  —  zJG^Cr-r'') 

Performing  the  z  *  integration  first 

oo  r 


.  00 

exp  l  —  z  f  exp 


— “(z  1  r — z  ») 
4s 


=  exp ( i  ^  z ' ' ) 


_  OO 


(5.15) 


exp(- 


s)  (47Tis) 


(5.16) 


is  the  isotropic  radiation  propagator 


G1(r-r  '  ') 


exp.  (i  ^  | r-r ' ' | ) 
477  |  r-r  1 1 


(5.17) 


where  r  is  the  point  at  which  the  field  is  measured  and  r '  '  is  the 
source  point,  relative  to  a  suitable  chosen  origin. 

As  we  are  interested  in  the  fields  at  large  r  and  in  the 
radiation  through  a  sphere  at  infinity,  let  us  take  the  asymptotic 
form  of  G^  for  the  case  |r|  >>  | r '  ' | .  Then,  Equation  (5.17)  becomes 
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exp(i  ^  r)  exp(-i  ft  .  r  '  ') 

Gx  (r-r  ')  (large  r)  =  - 4rrTrl -  (5.18) 

where  ft  is  a  unit  vector  in  the  direction  of  ir.  Then 

_(4)  ^o  Y  +  1  0)3  exP(i  5  |r|  a 

E  _  ___  v-  z 


^7  Jdz '  'dx  '  'dy  '  '  exp  i 

(5.19) 


foo 2  co  .  rT  '  'Y 

~~ 2  -  ^)S  -  “47“ 
a  v 


exp  l 


2  z’ ■  -  S  ft 
v  a 


Next  perform  the  z ' f  integration  which  leads  to  a  6-function 


* 

dz  f  f  exp. 


z'  1  - 


CO 

a 


n  z  f  1 
z 


27T6  (a  . ;  ftz) 


(5.20) 


where  nz  is  the  z-component  of  the  unit  vector  n  by  z.  The  x'  f  and  y !  1 
integrals  are  similar 


I- 


dx  r  1  exp 


T  ./Y  .  ,2  oo  a  ,  _  ,4ttsx1/2  ,  .  a/  a  2  s, 

|-^4ix"  ‘I  VC")J  =  expd  nx  -) 

A  2  ^ 

The  y  '  '  integral  yields  a  similar  result  with  ft  replaced  by  ft  . 


(5.21) 

2 

"X  --r  ~y 

Substituting  the  results  of  these  three  integrations  into  Equation  (5.19), 
we  obtain 
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f  .  CO  I  I  .  00  f 

+  expd  -Irl)  A  f 

Y - [71 -  z  Jds  exp  1S< 


2  2  2 

foo  oo  .  o>_ 

'  2  ~  2}  2 

a  v  a 


r ^  2,  a  2« 

n*  +ny  1 
Y 

✓ 


»  *C0  CO  A  < 

6  (-  -  —  n  )  = 

v  a  z 


Y^L  » 


"'i  |  «  «  ,60  60  _ . 

-i-d)  “  COS  0) 

FI  Za7ST7 

2^  V  '  2 

a  1  v 


(5.22) 


where  we  have  replaced  nz  and  +  n  ^  by  cos  0  and  sin^O  respectively. 

Equation  (5.22)  gives  the  electric  field  due  to  the  term  containing  G . . 

1  * 

It  lies  entirely  in  the  z-direction  and  falls  of  asjpp  like  a  typical 
radiation  field.  It  should  be  noted  that  the  only  approximation  used 
in  obtaining  this  result  was  to  replace  the  isotropic  propagator  by 
its  asymptotic  value  for  large  r  since  we  are  concerned  with  field 
measurements  taken  at  points  a  great  distance  from  the  origin  (or  source) 
Otherwise,  the  calculation  is  exact. 

Next,  we  must  calculate  the  scalar  potential  0  which  represents 
that  part  of  the  electric  field  ;due  to  the  Green's  function  G^ . 


0  =  dz  *  G^  (xyz-z  ')  exp(i  ~  z  ') 


(5.23) 


Again,  using  the  integral  expression  for  we  can  perform  the  zT 
integral  to  give 


=  •(4jriJ/2/-  exp(i  ^  z)  Jds 


T.  ,6o2  co2  1  ,  .  rt  Y 

exp  s  exp(-i  -j—) 

-  a  v  J 


(5.24) 


Leaving  ^  in  this  integral  representation,  the  electric  field 
can  now  be  written  as 
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.1/2 

i  '  I 


E(ro3)  = 


^  CO2 


47 T 


03 


2  v 


S-7 


exp(i  ~  z) 


00  r  ' 

fds  .  m 

J—  expih 

o  La 


-  ffij)  sj  exp  ( -  — *— ) 


si 


irT  y 


J-iQq<  Y  +  1  a  exP(i  flrl)  6(v  "  a  cos  °) 

+  —T-  y”  z  — na - - - - a 


,03 


00 


2v 


2  2 

oo"  n  sin  Q.  Q) 

2U  y  2 

a  T  v 


(5.25) 


CALCUALTION  OF  THE  MAGNETIC  FIELD 

The  magnetic  field  is  obtained  from  the  electric  field  by  use  of 
the  Maxwell  equation 


-  Z  x  E 

ioo 


(5.26) 


(4) 

As  mentioned  earlier y  only  E  contributes  to  B  because  the  curl  of  a 
gradient  is  zero. 


Therefore , 


iooB  =  Vxz\}r  (5,27) 

where  f  is  the  magnitude  of  E^.  But  using  the  well  known  vector 
identity 


Vxz^J/  =  -  S  x  Vj 


(5.28) 


and  keeping  only  the  leading  or  1/ [ r  t  term  after  differentiating 

(which  brings  down  a  factor  of  — ) ,  we  obtain 

a  ' 

B  =  -  z  x  (5.29) 

a 

where  n  is  a  unit  vector  in  the  r-direction.  Explicitly 
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qw3  y  +  1  exp(i  ^|r [ ) 


2v2a 


Y 


j.  /CO  00 

8  (-  -  -  cos  0) 

V  Q. 

2  2  2 

co  ,,  sin  __  co 

2U  y  }  2 

a  '  v 


(5.30) 


where  x  and  are  unit  vectors  in  the  respective  directions. 

CALCULATION  OF  THE  TOTAL  POWER 

The  Poynting  vector  S(rto)  which  denotes  the  total  energy  radiated 
per  unit  area  is  constructed  from  the  electric  and  magnetic  fields 
according  to  the  prescription 


S(nco)  =  E(roo)  x  B*(roo) 


(5.31) 


The  magnetic  field  is  entirely  due  to  the  Green's  function  . 

The  electric  field  consists  of  terms  involving  G^  and  G^ ,  namely 

and  respectively.  Therefore,  the  Poynting  vector  will  contain 

terms  involving  the  products  G^G^  and  G^G^.  We  shall  show  that  the 

total  power  due  to  the  terms  involving  G^G.  vanishes  over  a  surface 

(4)  6  4 

at  infinity  and  that  only  the  E  part  of  the  electric  field  con¬ 
tributes  to  the  radiation.  The  Green 's  function  G^  represents  pro¬ 
pagation  inside  and  damping  outside  of  the  Alfven  cone  of  small 
angle  0  -  tan~^  j  .  We  expect  that  only  the  z-component  of  the 
part  of  the  Poynt?ng  vector  involving  G^  will  contribute  to  the  total 
power,  since  the  transverse  components  lying  wholly  outside  the  conical 
waveguide  will  be  completely  attenuated  at  large  enough  distances. 
Therefore,  we  will  construct  S  ,  the  z-component  of  the  Poynting 
vector,  which  is  independent  of  terms  involving  the  product  G^G^ , 
integrate  this  expression  over  the  infinite  plane  parallel  to  the 
xy  plane,  at  z=zq  and  then  show  that  this  latter  quantity  vanishes 

as  z  — 
o 


i 
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(5.32) 


dRoo  =  —  Jnz#  S  dx  dy  is  the  total  energy  radiated  per  unit  frequency 
through°an  infinite  plane  parallel  to  the  xy  plane  at  z=z  . 
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exp- (- 


(5.33) 


The  6-function  imposes  a  condition  on  0.  Let  us  make  the  double 
change  of  variable  from  x  ,y  — >  r  $  and  then  another  change  from  r  — >  rT 
where  of  course 

x  =  r  cos  ^  sin  0  (5,34) 


64 


S-2023-1 


y  =  r  sin  </>  sin  0 


(5.35) 


-4 


2  2 
r  ~-M rT2  +  zQ2 


and  dxdy  — >  rdrd0  — > 


(5.36) 

(5.37) 


Furthermore 


A  4.  9 

x  dx  5y 


T  9r„ 


(5.38) 


There  is  no  azimuthal  angle  dependence  and  the  integration  of  0 
yields  2 tt. 

Manipulating  the  6-function  which  now  involves  rT  (since  cos 


0  =  i — ^ - 2}  according  to  the  formula 


rT  +  Zo 


b  (f (x) )  = 


6  (x-x  ) 

“ H 

9  x 


where  f (x  )  =  0 
o 


(5.39) 


leads  to 


co 
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oo  v  o  r 
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oo 


(rt2+zo2^ 


-i  r^/z.2+rm2)  (z  2+r  2) 
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(5.40) 
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where 


and 


v 


2 


2 


v 


(5.41) 


(5.42) 


Therefore 


dR 


00 


12 
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_O0 

00 


3  2  2 

_  2_  2_ 
2  ya  3 
'  v 


°°  i  z  2  2 

exp(i  ~  zQ)  zq2  J  exp(iv2s)  exp(-  -1) 


o  s 


(5.43) 


We  are  primarily  interested  in  the  dependence  of  the  total  power  on 
zq  ,  for  we  wish  to  show  that  the  above  expression  will  vanish  as  we 
remove  the  plane  to  infinity. 

The  s-integral  is  of  the  form 

°°  2  .2 
J  ^  exp(iv2s)  exp.  (-  J  dt  exp(-i62) exp(-^jf— ) 

os  ° 

(5.44) 


where  the  change  of  variable  s  -  T  has  been  made  and 


2  .  V  ,v2 

.  -  —  (“J  -  1) 

a 


Now  effect  another  change  of  variable 


6  t  -  x 


(5.45) 


(5.46) 
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and  the  integral  becomes 

oo  _ 

,  r  .  2.  2 

"2  J  dx  exp  (-  ix)  exp.  (“j — )  (5.47) 

o 

1  12 

The  — r  in  front  of  the  integral  contains  a  — ~  which  cancels  the  z 
1 1  £  o 

o  zo 

in  Equation  (5.43) 

Thus ,  the  entire  zq  dependence  of  the  total  power  (except  for  the 
harmonic  term)  is  contained  in  the  integral  in  (5.47) 

If  the  path  of  integration  is  extended  into  the  complex  plane  to 
include  the  infinite  quarter-circle  in  the  lower-right  hand  plane, 
there  is  no  added  contribution  to  the  integral  along  this  path  since 
if  x  is  replaced  by  x-iy,  y  positive,  exp(-ix)  — >  exp(-ix)  exp(-y) 
which  approaches  zero  as  y  — >°°,  and  the  other  term  in  the  integrand 
is  bounded.  This  argument  is  analogous  to  that  given  for  extending 
the  contour  when  evaluating  the  integral  expression  for  in  Chapterll. 
Then,  since  there  are  no  poles  contained  within  the  path  of  integration, 
we  can  deform  the  contour  and  integrate  along  the  negative  imaginary 
axis.  This  is  equivalent  to  replacing  x  by  -iy  and  integrating  from 
o  to  co.  Thus,  the  integral  becomes 
00 

2  2 

-i  J  dy  exp(-  y)  exp(-  -)  (5,48) 

o  y 

from  which  it  is  readily  seen  that  the  expression  tends  toward  zero 
2  2  2  2 

as  v  b  gets  very  large.  For  if  v  b  -  0,  the  integral  has  the  value 
2  2 

1  x,  and  as  v  6  starts  to  grow,  it  diminishes  the  y  =  0  contribution 

2  2 

from  exp(-y)  .  But  b  is  proportional  to  zq  ,  and  thus  as  zq  — ? 
both  the  integral  in  Equation  (5.48)  and  the  expression  for  dR^  the 
remainder  of  which  is  independent  of  zq ,  tend  to  zero. 

Therefore ,  that  part  of  the  Poynting  vector  which  contains  the 
product  ,  when  integrated  over  an  infinite  plane,  parallel  to 

the  x-y  plane,  at  zQ ,  gives  no  contribution  to  the  total  power  when 
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the  plane  is  taken  at  very  large  distances  from  the  source.  Thus, 
the  only  part  of  the  Poynting  vector  contributing  to  the  total  radiated 
power  is  that  which  contains  the  product  G4G4. 

The  remaining  portion  of  the  Poynting  vector,  which  we  shall  denote 

(4) 

as  S  ,  results  in  a  uniform  flux  through  a  sphere  even  when  the  radius 
of  the  sphere  is  infinite.  This  flux  may  be  considered  radiation,  since 
it  is  energy  lost  to  the  system.  As  we  shall  see,  the  total  power 
radiated  is  a  rather  small  quantity,  but  it  seems  to  be  the  dominant 
process  for  removing  energy  from  the  incoming  particle,  since  the  power 
lost  by  synchrocyclotron  radiation  at  nonrelativistic  velocities  is 
negligible,  and  the  Bremsstrahlung  effect  is  absent. 

(4)  (4) 

S  ,  which  depends  on  E  and  B,  contains  a  double  6 -function. 
Then,  as  was  discussed  in  Chapter  V,  the  integral  over  the  double 
6-function  yields  a  time  T,  which  is  the  value  of  27f6  (0)  ,  accounting 
for  the  fact  that  the  current  source  is  not  on  for  all  time  but  over 
an  arbitrarily  large  time  interval  T. 


S(4)  =  E(4)  x  B*: 
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(5.49) 


and  integrating  over  a  sphere  of  radius  r  gives  the  total  energy 
radiated 


dl  -  — 
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There  is  no  azimuthal  dependence  and  the  ^-integration  yields  2tt. 

T  is  the  actual  value  of  b  (0)  and  represents  the  time  interval  over 
which  the  current  acts. 

(4) 

Since  S  is  a  function  of  r  and  oo  ,  dl^  represents  the  total  energy 
radiated  rather  than  the  total  power.  To  obtain  the  latter,  we  must 
divide  the  above  expression  by  the  time  interval  T.  Thus, 


ooa4  1 

2  3  2  2 
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b  (cos  0  -  — ) 
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d(cos  0 


(5.51) 


and 
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u  4 
ro  c _ 

2  3 
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oo^doo 


b  (cos  0  - 


d(cos  0) 


(5.52) 


Again,  as  in  the  normal  Cerenkov  effect,  there  is  no  radiation 
unless  v  is  greater  than  a(oo).  Then,  all  of  the  radiation  lies  in 
the  surface  of  a  cone  of  angle  0  =  cos~^  ~  around  the  direction  of 
the  particle. 


The  total  radiation  per  unit  time  is 


oo 


=  !w.ai  r 

2  v3  »  4  A 


1  - 


a 


(5.53) 


Here ,  the  total  radiated  power  varies  as  the  sixth  power  of  the 

frequency  as  opposed  to  the  normal  Cerenkov  case  when  the  frequency 

2  2  2 
dependence  went  as  (x>  .  Furthermore ,  the  factor  of  v  -  a  ,  which 

appeared  in  the  numerator  in  the  normal  Cerenkov  case  now  is  found 

in  the  denominator.  This  suggests  that  the  radiation  becomes 

arbitrarily  large  when  v  — >a.  But,  when  v-a ,  cos  0=1  and  the 

radiation  would  be  directed  along  the  magnetic  field  lines.  From 
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the  symmetry  of  the  problem,  one  expects  the  radiation  to  fall  to  zero 
along  the  direction  of  the  incoming  particle.  We  must  look  more  closely 
at  the  fields  to  determine  the  correct  expression  for  small  angles,  0, 
in  the  limit  of  v  — >a. 

In  the  calculation  of  E  the  asymptotic,  large  |r|,  far-field  ex¬ 
pression  for  was  used.  This  asymptotic  form  is  valid  only  if  the 

source  term  mutiplying  G-,  in  the  integrand  is  localized.  In  the  normal 

ext 

Cerenkov  effect,  the  source  term  is  merely  the  current,  j  =  qv5(z-vt) 
which  is  indeed  localized.  However,  in  the  present  case,  the  source 
term  is  the  product  G^j  which  is  not  localized  since  G^  is  singular 
on  the  surface  of  the  Alfven  cone.  Thus,  the  asymptotic  approximation 
for  G^  should  be  valid  at  large  distances  provided  that  we  remain  also 
outside  of  the  cone.  These  expressions  are  incorrect  for  very  small 
angles  lying  outside  the  Alfven  cone  and  an  alternative  method  must 
be  used. 

(4) 

Returning  again  to  the  expression  for  E  ,  Equation  (5.10)  , 
let  use  for  G 4  the  integral  expression  developed  in  Chapter  III, 

.-3/2  ^  ,  2  2 

G4(xyz)  =  Jo  ^172  exp(i  ^2 u+i  fc) 
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(5.54) 
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zxVE 


(5,56) 


It  will  be  more  convenient  here  to  work  with  B  rather  than  E. 


6E(4)  6E(4 

and  B  will  involve  — r -  and  — t — 

y  oy  dx 


and  the  taking  of  these  derivatives 


will  enable  us  to  do  the  C-integral  and  also  the  u-integral( 


i  JL  f(4)  _  2ix  _9 p(4) 
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(5.57) 


and  a  similar  expression  for  B  with  y  replaced  by  -x.  The  only  r_ 

(4)  y  1 

dependence  in  E  is  contained  in  the  ^-integral. 
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(5.58) 


Denoting  the  u- integral  by  I(u) ,  it  becomes 
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(5.59) 


By  making  the  substitution  u  “  1/x,  each  of  these  integrals  can  be 
put  in  the  form 
oo 

P  dx 

J  "  3/2  exP^  Xx)  exP(i  — )  whose  solution  is 

o  v  '  ^  4  X 
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Therefore 


Ku)  =# 


exp ( i  ^  y (z-z  ')  2+r  2) 


-  exp(i  ^  ^(z-z  ')  2-rT2y) 


(5.60) 


Then,  from  Equations  (5 . 60)  ,  (5 . 57)  and  (5.55) 
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(5.61) 


We  now  assume  that  the  external  current  is  localized  within  the 


2  2  2 

Alfven  cone,  i.e,,vt  <z  -  r, 


Then 
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Performing  the  t-integral ,  we  obtain 
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2  2  1/2 

where  again  |r|  -  (z  +r_  ) 


9  9  V2 

and  I r I  -  (z  -rT  y) 


The  leading  term  remaining  in  each  expression  after  taking  two 

derivatives  is  just  the  original  expression  multiplied  by  the  factor 

2  2  2  2 
GO  Z  GO  Z  — 

-  “2  — 2  f°r  the  term  involving  r  and  -  —  for  the  term  involving  r. 

a  r  r 

We  are  interested  in  the  radiation  part  of  the  field  and  hence  keep  only 
those  terms  that  exhibit  a  1/r  dependence. 
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(5.65) 


In  order  to  see  the  angular  dependence  of  the  field,  let  us  convert 
to  spherical  coordinates,  r,  0, 

The  6 -function  appearing  in  the  first  term  can  be  written 


b  (w-00  — 


ITT 


)  =  5  (cos  0 


a»  _a_ 
v;  vco 


(5.66) 


and  the  other  6 -function  can  be  cast  into  the  form 


6  (co-co  J  -ft) 
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(5.67) 
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making  use  of  the  relation 

b  (x-x„) 


Mf  (X)) 


C4 


where  f (x  )  r  0. 
o 


Then,  making  use  of  Equations  (5.66)  and  (5.67)  in  (5.65) 
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(5.68) 


The  two  terms  in  this  expression  are  mutually  exclusive.  If  one 
contributes  to  the  magnetic  field,  the  other  does  not.  The  first  term 
contributes  only  for  v/a  >  1,  since  the  argument  of  the  6-function 
vanishes  when  cos  0  a/v.  However,  the  second  term  contributes  when 
v/a  <  1  as  is  easily  seen  from  examining  its  6-function.  Thus,  the 
first  expression,  resulting  from  the  exp(i  co/a|r[)  term  in  ,  is 
analogous  to  the  normal  Cerenkov  effect.  The  magnetic  field  is  con¬ 
tained  in  a  cone  of  angle  0  ~  tan-1  a2~~  ,  and  occurs  only  when 

c  a 

the  particle  velocity  is  greater  than  the  phase  velocity n  The  second 
term  leads  to  a  strikingly  new  effect  which  we  shall  henceforth  call 
the  anomalous  Cerenkov  Effect.  Here  the  radiation  field  is  contained 
within  a  cone  of  angle 


0  f  “  tan  ( 


i  2,2 

1-v/a 


1/2 


and  occurs  only  if  the  particle  velocity  is  less  than  the  phase 
velocity.  The  highly  anisotropic  plasma  is  capable  of  sustaining 
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a  radiation  field  regardless  of  the  velocity  of  the  moving  charged 
particle,  whereas  normal  Cerenkov  radiation  is  possible  in  an  isotropic 
medium  only  when  the  moving  charge  traverses  the  medium  at  speeds 
greater  than  the  phase  velocity  of  light  in  the  medium.  Moreover,  note 
that  the  cone  accompanying  the  anomalous  Cerenkov  radiation  always  lies 
within  the  Alfven  cone  and  coincides  with  it  m  the  limit  of  v/a-yO. 
Mathematically,  this  is  a  consequence  of  the  fact  that  this  new  effect 
arose  from  the  conical  waveguide  term  exp(i  ^  \r\  )  appearing  in  G^.  Pre¬ 
viously,  this  term  did  not  appear  because  of  our  inability  to  treat  the 

Alfven  cone  properly.  Its  physical  manifestation  is  of  unknown  origin „ 

-1/2 

The  overall  factor  of  y  '  multiplying  this  term  partially  explains 
the  fact  that  such  radiation  has  never  been  detected.  For  go  =  lOsec*"^, 
the  amplitude  of  the  magnetic  field  associated  with  the  anomalous  effect 
is  10,000  times  smaller  than  the  field  giving  rise  to  the  normal  Cerenkov 
radiat ion 

However ,  it  is  evident  that  our  treatment  of  the  problem  is  not 
adequate  for  very  small  angles.  Note  that  as  ©—*-(),  each  term  becomes 
arbitrarily  large  as  v  approaches  a  form  above  and  below,  respectively. 
This  results  because  each  term  is  proportional  to  esc  0  which  in  each 
case  is  dependent  on  the  quantity  - - -  .  As  will  be  shown  shortly. 

*  (4) 

the  expression  for  the  electric  field  E  will  also  have  two  terms 
which  are  both  independent  of  sm  0  and  therefore  finite  as  0 — V0o 
However,  the  product  of  E  and  B,  in  the  Poynting  vector,  will  result 
m  a  finite,  nonzero  total  power  for  Q  =  0.  But,  we  expect  zero 
power  as  we  look  along  the  path  of  the  moving  charge  because  there 
is  no  transverse  direction.  Moreover,  the  Green’s  Function  G^  is 
known  exactly  for  the  case  of  0  =  0  (rT  =  0) .  Thus,  we  can  calculate 
the  exact  expressions  for  the  electric  and  magnetic  fields,  and  check 
that  the  power  really  is  zero. 

exp  (i  (z-z  ’) ) 

G4(rT=°> - -  1^7  r  (5-69) 
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( 5  0  70) 


which  is  not  a  radiation  field  since  it  does  not  exhibit  the  classical 

1/r  behavior,  but  rather  more  like  a  plane  wave  solution.  The  magnetic 

field  for  small  angles  can  be  obtained  by  returning  to  Equation  (5.61) 

and  expanding  exp(i  —  |rl)  -  exp(i  — [ r [ )  in  a  Taylor  series  about  small  r^, 

a  a  l 

Then 
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which  goes  to  zero  as  y  — >  0.  Likewise,  which  is  proportional  to 

x  also  vanishes  in  this  limit.  The  magnetic  field  thus  vanishes  for 

0  -  0  because  when  treated  correctly,  the  two  terms  contribute  together 

and  oppose  each  other.  In  Equation  (5.65)  the  two  terms  cannot  cancel 

each  other,  even  in  the  limit  of  0  — >  0,  because  they  do  not  exist 

simultaneously,  i.e.,  the  existence  of  one  precludes  the  existence  of 

the  other.  Yet,  when  the  problem  is  treated  properly  for  small  rT , 

we  see  that  the  two  terms  contribute  together  and  the  net  result  is 

zero.  This  is  to  be  expected  since  when  0 ,  r  r  r  ~  z  and  the 

two  terms  are  identical.  The  inconsistency  present  here  is  a  consequence 

of  our  inability  to  treat  the  isotropic  propagator  exp(i  — [ r [ )  when 

a 

[r|  is  im ide  the  Alfven  cone  and  0  is  near  zero.  The  mechanism  that 

prevents  the  field  calculated  in  Equation  (5.67)  from  growing  arbi- 

—1  /  2  2 

trarily  large  comes  into  play  when  the  0  given  by  0  -  tan  Va  -v 

a 

lies  within  the  Alfven  cone.  For  it  must  be  noted  that  even  though 
the  medium  behaves  as  a  conical  waveguide,  there  is  nothing  in  the 
normal  Cerenkov  term  to  denote  our  passing  from  the  damping  region 
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outside  of  this  cone  to  the  region  of  propagation  inside.  We  would 
expect,  if  the  calculation  were  performed  correctly,  that  the  normal 
Cerenkov  term  would  exhibit  a  change  in  character,  as  0  is  decreased 
to  lie  within  0  .  Bearing  in  mind  the  fact  that  the  magnetic  field 

a 

vanishes  for  Q  =  0  in  view  of  the  exact  calculation  and  the  heuristic 

arguments  given  above,  we  can  use  Equation  (5.67)  to  describe  the 

field  at  all  points  off  the  magnetic  field  axis.  Since  for  ELF  waves 

-4 

Qa  is  of  the  order  of  10  radians,  the  normal  Cerenkov  term  is  correct 
down  to  the  smallest  angles.  In  addition,  since  the  radiation  accom¬ 
panying  the  anomalous  effect  is  contained  within  the  Alfven  cone,  this 
radiation  would  appear  to  be  concentrated  in  infinitesimally  thin  tubes 
centered  about  the  magnetic  field  lines. 


(4) 

In  order  to  determine  the  electric  field  E  ,  we  use  still 
another  integral  representation  of  G . .  In  Chap.  Ill  we  developed 
an  expression  for  involving  an  integral  over  rT 
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(4) 

where  E  '  is  the  value  of  the  electric  field  for  ©  or  r_  =  0 
o  T 

m  C  4 )  __  ^  o  co  a  t/i  a  \  f  •  ^  r  n  j  t 

*0  “  "4 - 2  6(1  "  v}  exP(l  (5.74) 

V 

Again  we  assume  that  the  source  is  localized  within  the  Alfven 
cone  ,  i.e.  , 


24.2  ..2  2 
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and  expand  the  exponentials  as  before.  This  leads  to 
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In  order  to  perform  the  remaining  integrals,  we  write  the 

2  2 

5-functions  in  the  form  b  (u  -uq  ) 
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From  the  previous  calculation  of  the  magnetic  field,  we  know 

that  the  first  6 -function  contributes  only  when  v  is  greater  than 

2 

a,  and  the  latter  only  when  v  <  a.  Furthermore,  the  value  of  r,p 
in  the  limits  of  the  integration  must  be  large  enough  in  each  case 
so  that  the  point  where  the  6-function  is  nonzero  lies  within  the 
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limits  of  the  integral.  Then 
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where 
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and  zero  otherwise,  etc.  In  taking  two  z-derivatives  of  the  quantity 
in  the  brackets,  we  keep  only  those  terms  that  have  the  1/r  dependence 
required  for  a  radiation  field.  The  terms  of  interest  are  the  cross 
terms  obtained  by  differentiating  once  both  the  exponentials  and  the 
H-f unctions.  The  derivative  of  an  H  function  gives  a  6-function  of 
similar  argument,  i.e. , 
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where  sgn  z  denotes  the  sign  of  z. 


Furthermore , 
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Note  that  the  6-functions  appearing  here ,  when  written  in  spherical 
coordinates  are  the  same  6-functions  contained  in  the  expression  for 
the  magnetic  field. 


Thus,  taking  the  appropriate  derivatives  and  rearranging  the 
resultant  6-functions  and  converting  entirely  to  spherical  coordinates 
leads  to 
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Again  both  the  analogue  to  the  normal  Cerenkov  effect  and  the 
anomalous  Cerenkov  term  are  present  in  this  expression.  However, 
unlike  the  result  derived  earlier  for  the  magnetic  field,  the  elec¬ 
tric  field  does  not  become  arbitrarily  large  as  v  — >  a  from  above 

and/or  below  but  merely  approaches  a  finite  value.  If  Equation  (5.85) 

(4) 

were  entirely  correct,  the  field  would  approach  the  value  E  (0) 
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but,  as  can  be  easily  seen,  the  terms  involving  the  6-functions  do 

(4) 

not  vanish  as  v  ->  a.  As  a  consequence  to  the  fact  that  E  '  does 
not  reduce  to  its  true  value  as  rT  =  0,  the  total  power  will  not 
vanish  for  this  same  case,  but  instead  approaches  some  constant  value. 
Actually,  the  electric  field  at  rT  =  0  is  not  a  radiation  field  but 
a  plane  wave  solution  of  the  form  exp(i  ^Iz|)  and  hence  should  not 
contribute  to  the  total  power.  The  6-functions  terms  in  Equation  (5.85) 
keep  their  1/r  dependence  even  in  the  limit  of  r^  — >  0.  The  reason 

for  this  discrepancy,  as  stated  above,  is  not  known,  and  is  contained 
in  the  mathematics.  If  the  electric  field  depended  on  sin  Q,  the  term 

/\J. 1  -  ^~2  would  appear  in  the  numerator ,  and  the  r^,=  0  value  (v  =  a) 

of  the  field  would  then  reduce  to  its  true  plane  wave  value. 

In  computing  the  total  power,  we  shall  treat  the  two  effects  sep¬ 
arately. 


1.  ANALOGUE  TO  NORMAL  CERENKOV  EFFECT  -  v  ^  a. 


The  fields  contributing  to  this  effect  are  obtained  from  equations 
(5.85)  and  (5.68) 
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Again,  one  of  the  6-functions  is  identified  with  T,  the  overall 
time  that  the  source  is  on. 
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2TT  6  (0)  =  COT 


(5.89) 


The  total  energy  radiated  per  unit  frequency  is  given  by 
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and  the  rate  of  radiation  is  just  the  above  expression  divided  by  T. 

The  total  power  radiated  by  the  moving  point  charge  is  then 
obtained  by  integrating  over  all  solid  angles  and  over  the  frequency 
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The  integral  being  over  those  frequencies  in  the  ELF  range  for  which 
v  >  a  (v)  . 


Both  the  total  power  and  the  power  per  unit  solid  angle  are 

2  2 

independent  of  the  factor  v  -a  and  furthermore  are  constant  for  a 
given  v. 


The  normal  Cerenkov  radiation ,  as  given  by  Equation  (4-32)  is 
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which  for  the  case  of  v  >>  a,  reduces  to  the  expression  obtained 
for  the  present  case  (except  for  a  factor  of  four.) 

This  result  is  not  too  surprising.  Our  expressions  for  the 
fields  are  most  correct  in  the  region  far  outside  of  the  Alfven 
cone,  where  the  particle  velocity  is  much  greater  than  the  phase 
velocity.  Obviously,  the  above  expression  for  the  total  power, 
Equation  (5.91)  ,  does  not  go  to  zero  as  v  — >  a  ,  but  approaches  a 
finite  constant  value.  One  can  conjecture  that  a  factor  of 
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2  2 

v  -a  in  the  numerator  would  suffice  to  yield  an  expression  compatible 

with  the  Q  =  0,  v  — >  a,  case,  but  there  would  still  be  no  reference  made 

to  the  existence  of  the  Alfven  cone.  In  fact,  because  of  the  nature  of 

the  Alfven  cone,  one  should  expect  the  radiation  to  be  enhanced  as  the 

Cerenkov  cone  first  coincides  with  the  former.  The  radiation  should 

then  reach  a  peak  for  some  value  of  <  ©^ ,  and  then  decrease  to  zero 

c  a 

at  0  -  0.  The  mathematical  mechanism  that  does  this  is  absent  from 
the  expression  for  the  magnetic  field.  When  the  ratio  of  a/v  is  such 
that  the  two  cones  coincide,  the  magnetic  field,  as  given  by  Equation 
(5.65)  ±sJ\fx'  times  as  great  as  its  value  for  large  v  >>  1).  This 
represents  a  large  buildup  of  the  field  strength.  However,  as  we 
approach  ©c  =  0 ,  this  expression  becomes  infinite  instead  of  cutting 
off  at  some  value  of  ©c  at  inside  the  Alfven  cone.  A  summary  of  these 
results  will  be  presented  at  the  conclusion  of  this  chapter. 


2.  THE  ANOMALOUS  CERENKOV  EFFECT 


Looking  now  at  the  fields  produced  by  the  anomalous  effect,  the 
terms  of  interest  are 


exp ( i  j I r 1 ) 
2  »rrl  (l+y-v2/a2) 


^oq  v 


6  I  cos  0  -a  I - ¥ — —  \  (sin  f6  x-cos 

1  Vl+Y-vVa*  ) 


and 


,(4)  _ 


—  xl  a~ 

"4  3  qZ 

a 


bf  COS 

exp  (I  girl)  V 


JL 


i+r-  ^ 


1  +  Y  -  *5 


(5.94) 


where  now  v  <  a.  This  effect  does  not  occur  when  a  uniformly  moving 
charged  particle  traverses  an  isotropic  medium.  The  accompanying 
Poynting  vector  is 
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where  again  one  of  the  6  -functions  has  been  associated  with  the  time 
T  and 
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gives  the  radiation  over  all  frequencies  in  the  ELF  range  for  which 

v  <  a(v).  Note  that  for  all  extensive  purposes,  this  expression  is 

proportional  to  1/y  and  thus  is  a  very  small  quantity  compared  to 

8  —1 

the  ordinary  Cerenkov  result,  (y  =  10  for  =  10  sec  )  Thus  the 

radiation  for  all  values  of  v  <  a  is  contained  in  the  surface  of  a 
i  -4 

cone  of  angle  9c«10  radians  and  is  of  very  low  intensity.  For 
v/a^l,  the  radiation  from  a  single  charged  particle  for  all  ELF 
frequencies  up  to  100  cps  is  of  the  order  of  1CT^  watts.  Even  a 
beam  of  particles,  each  contributing  collectively,  would  not  raise 
this  contribution  to  a  value  that  could  be  detected  above  the  back¬ 
ground  noise.  In  fact,  it  is  highly  doubtful  whether  either  the 
radiation  or  the  associated  magnetic  field  could  be  detected  by  any 
conventional  means,  assuming  such  an  attempt  was  made. 

SUMMARY  OF  RESULTS 

A  uniformly  moving  charged  particle,  traversing  the  ionosphere 
in  the  direction  of  the  lines  of  force  of  the  earth's  geomagnetic 
field,  produces  low  intensity  radiation  at  ELF  frequencies.  If 
the  particle  velocity  is  greater  than  the  phase  velocity,  the 
radiation  is  contained  in  the  surface  of  a  cone  of  angle  ©c=  cos“^(a/v). 
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For  large  values  of  v/a ,  the  radiation  and  the  associated  fields  are 
nearly  identical  with  those  produced  by  the  normal  Cerenkov  effect 
which  is  observed  in  isotropic  media.  As  v  approaches  a  from  above, 
the  magnetic  field  builds  up  to  very  large  values  and  reaches  a  max¬ 
imum  for  some  value  of  —  <  ,  I  y  .  Because  of  improper  treatment  of 

v  V  Y+l 

the  problem,  the  expression  derived  for  the  magnetic  field  becomes 
arbitrarily  large  as  v  — >  a  and  the  radiation  remains  finite  and  non¬ 
zero,  whereas  it  should  vanish  because  of  the  symmetry  of  the  problem. 
The  part  of  the  electric  field  that  contributes  to  the  total  power  lies 
wholly  in  the  magnetic  field  direction  and  the  direction  of  the  Poynting 
vector  is  such  that  it  cannot  represent  the  radiation  from  a  transverse 
wave.  In  the  normal  Cerenkov  effect,  the  wave  is  transverse  in  char¬ 
acter  and  the  associated  electric  field  is  not  restricted  to  the  direc¬ 
tion  of  the  external  magnetic  field. 

Because  the  ionospheric  plasma  behaves  much  like  a  conical  wave¬ 
guide,  radiation  of  much  lower  intensity  is  produced  even  when  the 
particle  velocity  is  less  than  the  phase  velocity.  This  latter  effect 
does  not  arise  when  the  particle  traverses  an  isotropic  medium.  The 
radiation  is  contained  within  the  surface  of  a  cone  of  very  small  angle 
which  under  no  circumstances  lies  outside  of  the  Alfven  cone  of  angle 
0^  =  tan”'1'  —  Again  the  expressions  for  the  magnetic  field  and  the 

W 

total  power  are  incorrect  as  v  approaches  a,  this  time  from  below.  The 
intensity  of  this  effect  is  so  low  as  to  render  it  undetectable  by 
conventional  means. 

The  approximations  used  in  arriving  at  these  results  were; 

a)  collisions  were  neglected 

b)  plasma  considered  at  zero-temperature 

c)  ELF  approximation — wave  frequencies  less  than  the  ion  cyclotron 
frequency 

d)  source  localized  within  the  Alfven  cone. 
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CHAPTER  VI 
FINITE  TEMPERATURE 

The  inclusion  of  finite  temperature  leads  to  further  modification 

of  the  electric  susceptibility  tensor  which  in  turn  leads  to  temperature 

dependent  Landau  or  "thermal”  damping  of  the  ELF  modes  previously  discuss* 

(3) 

It  has  been  shown  that  the  purely  thermal  effects  of  collisions  are 
important  in  the  F-layer  where  the  temperature  of  the  electron  gas , 

T(°K)  is  equal  to  or  greater  than  2  x  lO-1^  a^(c.g.s)  ,  where  a  is  the 
Alfven  or  phase  velocity  for  MHD  waves  propagating  in  the  medium.  For 
the  zero-temperature  case,  the  various  diffusion  coefficients  (one  for 
each  species)  appearing  in  the  Fokker-Planck  operator  were  set  equal 
to  zero.  With  the  existence  of  finite  temperatures,  the  ratio  ~  appears 
both  in  the  kinetic  Green’s  function  and  the  ensuing  susceptibilities. 


THE  ELECTRIC  SUSCEPTIBILITY  TENSOR 


(3) 


In  the  report  of  Cantor,  Keilson  and  Schneider,  the  temperature- 
dependent  conductivity  tensor  was  derived  using  the  Fokker-Planck 
Equation  as  a  starting  point.  We  refer  to  that  report  for  details. 

The  result  of  interest  is 
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The  conductivity  tensor  is  related  to  the  susceptibility  tensor 
defined  by 


A 


X  = 


i 

€  CO 

o 


a 


(6-3) 


which  appears  explicitly  in  the  matrix  equation  relating  the  electric 
field  to  the  external  current (Equation  (2-27)). 
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and  A*kk*A,  which  appear  in  the  integrand,  are  both  matrices. 


cos  &T  sin  &T 

dA  -  I  -sin  Qt  cos  Or 

dT 


(6-5) 


A'kk*A  =  Q 

Ol  T 

Q 


sin2QT 


xx 


yy 


sin2&T 


k2  - 
x 


k‘  - 


(1-cos  Qt)  ,2 


Q4 


(1-cos  Qt) 
Q2 


y 

2 


(6-6) 

(6-7) 

(6-8) 


.  _2.  2 
Qzz  =  T  kz 


2  sin^T I 

Q: 


Qxy  \  o2 


■U-cosQT,).  +  !jLJZrsin2QT  +  (1.  :osQT)j^ 


kf,  sinQT(l-cosQT)  k  k 


yx 


xz 


Qd 


ri  Jv  r 

Q2  L 


sin  ^T  +  (1-cosQtM  (6-11) 


k  T 
z 

“ST 


k  sin&T  +  k  (l-cos&r) 
x  y 


(6-9) 

(6-10) 


(6-12) 


88 


S-2023-1 


®zx 

k  T 
z 

=  ir 

k 

X 

sinftT  -  ky(l-cos^T) 

(6-13) 

°yz 

k  T 
z 

=  TT 

k 

y 

sinQT  -  kx(l-cos&T) 

(6-14) 

°zy 

k  T 
z 

=  TT 

k 

y 

sinQT  +  kx(l-cosQT) 

(6-15) 

First  let  us  consider  the  case  of  zero  external  magnetic  field* 

In  the  absence  of  both  t]  and  Q,  the  matrices  B  and  A  are  diagonal  and 
isotropic 


B  = 


:  A  = 


T  1 


(6-16) 


Then 


oo 

J  dT  exp  (-iooT)  exp(-  ^  k2T2) 
o  z 


1  -  lit2  k  k 

/V 

(6-17) 


The  next  step  in  evaluating  this  integral  is  to  go  to  the  "low 
temperature"  limit.  The  low-temperature  limit  is  dictated  by  the 
electron  and  ion  temperatures  prevalent  in  the  ionosphere (>  1000°K) 
and  is  consistent  with  the  low  frequency  limit  under  which  we  are 
laboring.  In  order  to  expand  exp  ("^  k2T2)  in  a  power  series  and 
neglect  all  terms  after  the  second,  the  inequality 

|ik2T2  «  1  (6-18) 


must  hold.  But  T  is  a  proper  time  and  is  inversely  proportional  to 
0),  the  only  frequency  occurring  in  the  integrand.  Further,  fx  =  — 
is  just  the  square  of  the  thermal  velocity, 

Therefore,  each  species  must  satisfy 


^  > 


(6-19) 
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where  \  -  is  the  wavelength  of  the  ELF  oscillation  having  angular 
frequency  to. 

+  5 

The  thermal  velocities  for  02  ions  and  electrons  are  1*3  x  10 

and  3  x  107  cm/sec ,  respectively  assuming  a  temperature  of  1000°K  in 

the  F-region.  Even  m  the  exosphere ,  with  a  temperature  of  possibly 

3000°K.  the  thermal  velocity  of  electrons  does  not  exceed  5  x  107  cm/sec. 

Therefore ,  in  order  for  the  inequality  to  be  satisfied,  the  wavelengths 

7  -1 

must  be  greater  than  10  cm,  for  a  wave  with  frequency  go  =  30  sec 

However,  ELF  magnetohydrodynamic  waves  of  this  frequency,  traversing 

Q 

the  exosphere  at  the  Alfven  velocity  of  5  x  10  cm/sec,  do  indeed  have 
wavelengths  of  hundreds  of  kilometers  and  hence  the  two  approximations 
of  low- temperature  and  low  frequency  are  consistent.  For  the  low 
temperature  limit  actually  places  a  restriction  on  the  wave  number, 
k?  which  in  turn  leads  to  long  wavelengths  which  are  consistent  with 
the  additional  assumption  of  extra-low  frequencies.  Let  us  speak  of 
frequencies  below  30  sec  ^  as  lying  in  the  ultra-low  frequency  or 
ULF  band-  Therefore,  from  Equation  (3-38),  after  expanding  the  ex¬ 
ponential  and  performing  the  T-integrals: 

2  2 

A  -GO  GO 

X  =  — #  1  -  -f  n  (1  k2  +  2  k  k)  (6-20) 

CO*  *  GO*  ^  ~  ^ 


As  was  mentioned  previously,  possible  modes  of  radiation  are 
closely  related  to  the  normal  modes  and  an  investigation  of  the 
dispersion  relation  governing  the  latter  can  yield  valuable  infor¬ 
mation  concerning  the  former.  The  radiation-like  denominators 
appearing  in  the  various  Greenns  Functions  already  discussed  are, 
when  set  equal  to  zero,  nothing  more  than  the  dispersion  relations 
governing  the  normal  modes. 


The  matrix  equation  involving  the  electric  field  is  again 
2 


(k2  - 

c  ®  ~  ~  c 


Define  two  scalar  quantities  by 


S  =  K2  -  K 


co2 

-i+  _ 

c  c  CO 


2  2 


0  for  the  normal  modes 

(6-21) 

(6-22) 
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U  =  1 


2oo2 

~2^2  ^ 

c  oo 


(6-23) 


Then,  Equation  (6-21)  can  be  rewritten  in  terms  of  these  two  scalars 

(SI  -  U  k  k)  E  =  0  (6-24) 

cz 

The  necessary  and  sufficient  condition  for  the  existence  of  nontrivial 
solutions  is  that  the  determinant  of  coefficients  vanish.  Therefore 

det  (SI  -  U  k  k  )  =  0  (6-25) 


and  some  simple  algebra  leads  to 

det  (SI  -Uk  k  )  =  S2(S-Uk2) 


(6-26) 


Hence 


S 


2 


0 


(6-27) 


and 


S 


-  Uk2  =  0 


(6-28) 


The  double  mode  is  given  by 


S  =  k  + 


2  2 

or  a/  2 

-JB  +  _B  u 
2  +  2^2 
C  C  00 


00 


=  0 


(6-29) 


or  solving  for  k  in  terms  of  ao  in  the  limit  of  the  ULF  approximation, 

00  <V  CO 


,2  co- 

k  +  —  =  0 

V- 


(6-30) 
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which  is  analogous  to  the  damped  mode  described  by  the  Yukawa  potential 
which  was  previously  discussed.  The  third  mode  is  given  by 


00 


S  -  Uk  =  3 


2,2 

c  co 


|i.k 


O)2 

_J2 

2 

c 


00 

-~2=0 

c 


(6-31) 


and  the  resultant  dispersion  relation  is 

*2  +  3  F  =  0 


(6-32) 


which  again  represents  an  isotropic  damped  wave.  Therefore,  the  case 

2  2 

of  low  temperature,  in  the  limit  p,k  <<  co  ,  with  the  absence  of  an 
external  magnetic  field  results  in  two  independent  modes  both  of  which 
are  isotropically  damped.  The  damping  is  frequency  dependent  and 
depends  on  the  thermal  velocity  v^  =  10^  cm/sec. 

Returning  to  the  excitation  problem,  the  two  Green's  Functions 
for  the  low-temperature,  zero-magnetic  field  case  are 

/>  jj,.  _ik*r 

H^rco)  =  J 


d3k _ 

(2tt)  3  CO2  co2 

k2  (1+  -%■  J^)+ 

CO  C  c 


(6-33) 


1 

4Tr  |  r  | 


CO 

1  +  -J2_  J±_ 
,,2  2 

CO  c 


exp  i 


co  [  r  |  /c 

.  *r  — 


00 


1  + 


co 


1/2 


and 


Ho  = 


+  3 


exp 


-“Dlr f/c 

2  \  1/2 


co 

1+  3  Ji ±_ 
l  +  J  2  2 

00  c 


(6-34) 
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But 


CO2 

_E 

co2 


th 


D 

10  for  co  =  10  sec 


-1 


Therefore 


H, 


exp  - 


2LL£  I 

vth 


co2 

_e  jl 

,2  2 

CO  c 


47T  Ir  I 


and  H. 


is 

co2 


exp  - 


4*[r  [' 


00  lr  I 

Vth 


JL 

2 

c 


(6-35) 


as  constrasted  to  the  zero-temperature  case 

o^pfrf'j 

c  /  (6-36) 

71 

The  damping  in  the  latter  case  is  much  greater  than  the  thermal 

damping  occurring  when  finite  temperatures  are  included.  However,  the 

—8 

new  Green's  Functions  have  a  constant  factor  of  10~  multiplying  them, 
which  is  absent  from  G2*  Thus,  the  inclusion  of  low  temperature  results 
in  a  greater  damping  distance  but  the  amplitude  of  the  initial  oscillatic 
is  greatly  depressed,  and  the  wave  velocity  is  greatly  reduced. 


LOW  TEMPERATURE , FINITE  MAGNETIC  FIELD 

We  now  consider  the  case  of  a  finite,  external  magnetic  field, 

Bq,  whose  magnitude  is  contained  in  the  ion  and  electron  gyrof requencies 


Most  ELF  magnetohydrodynamic  waves  that  have  been  detected  at  the 
earth  have  been  observed  to  propagate  along  the  lines  of  the  earth's 
general  magnetic  field.  The  magnetic  field  direction,  has  a  low- 
frequency  conductivity  much  greater  than  that  in  any  other  direction 
and  serves  almost  as  a  waveguide.  For  every  possible  ELF  wave  here¬ 
tofore  discussed,  excluding  the  isotropic  mode,  the  propagation  is 
in  the  magnetic  field  direction.  Since  one-dimensional  propagation 
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along  the  field  lines  seems  to  be  the  recurrent  and  dominant  physical 
phenomenon,  we  shall  limit  the  ensuing  discussion  to  this  particular 
case  of  interest  by  imposing  the  condition 

k  -  k„  z  (6-37) 

on  the  propagation  vector. 


This  greatly  simplifies  the  calculations  for  the  matrix  Q  =  A.kkA. , 
now  has  only  one  nonvanishing  component, 


Q 


zz 


(6-38) 


as  can  be  seen  directly  from  Equations  (6-5) - (6-15) .  The  nonvanishing 
components  of  the  special  susceptibility  tensor,  X,  are  then  given  by 


zz 


O)2 

_E 

60 


OO 


V 


dTe 


iooT 


T 

2 


2  2 

(1-fik^O 


(6-39) 


A 

X 


xx 


OO 


r  icoT 

dTe 

% 


cos  £2  T 


(6-40) 


oo 


* 

V 


dTe 


icoT 


sin  £2  T 


(6-41) 


It  will  be  shown  shortly  that  each  of  the  above  integrals  can 
be  explicitly  written  in  terms  of  its  real  and  imaginary  parts,  and 
that  in  each  case,  one  of  these  parts  can  be  evaluated  exactly  whereas 
the  low-temerature  and/or  low  frequency  approximations  discussed  above 
must  be  invoked  to  effect  the  other  solution. 


For  example,  the  z-z  component  may  be  rewritten! 


2  oo 

A  00 

x  =  ■-£ 

zz  CO 

o 


2  t 


dTe 


2  2  2 
(-sin  cot  +  i  cos  cot)  (l-(xk  T  ) 
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=  M1+  ^t?)  / 


oo  ,  2  T 

-VK  2 

dTe  (-sin  oot  +  i  cos  cot) 

(6-42) 


The  imaginary  part  of  this  expression  is  easily  evaluated 

.2 

— 

4 


oo 


,  -  2  T‘ 

r  2  i 

J  cos  coTe  dT  =  t 


oo  2 

f  ^rr  T  2//2/  i^T  -iO)Tx 

J  dTe  (e  +  e  ) 


oo 


=  exp  - 


co 2  V  277  1/2 


2Hk, 


\^K 


(6-43) 


Therefore  from  Equation  (6-42) 


Xzz  =  -<1+ 


/  .  2, 2  \  oo 

L  ^kz  9  2 


^  000' 

2, -3/2 


J  sil 


2  /  ^p  J  sin  “T  e 

o 
2 


-[ikz  T2/2 


,2, 


nr 


dT  +  it o-]J  f 


(|xk  )  exp  - 


00 


2p-k 


(6-44) 


To  evaluate  the  real  part ,  we  go  to  the  low-temperature  approximation 

^*z2t2  «  1 


which  puts  a  restriction  on  the  wavelength  X  > 


2T rv 


th  i 


^  -(for  each  species)  , 
as  discussed  previously.  The  exponential  can  then  be  expanded  in  a 


power  series 


\  00  /  .  2  2n 
Re  Xzz  =  “p(X  +  ^kz2  z~2j  J  dT  sin  “V" 


000 


o2 

p 


2' 

a2i  ii 


oo2 

«;|i*  l£  + 

9oo  /  V  00/  oo 


(6-45) 
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This  is  precisely  the  zero  temperature  result  for  Xzz.  Thus,  the 

presence  of  a  finite  temperature  serves  to  make  X__  a  complex  quantity 

zz  A 

whose  imaginary  part  is  wave-number  dependent.  The  evaluation  of  X 

A  “ 

and  XW(B  follow  similar  lines 
zz 


i. 


dTe 


v  2  2 

-fxk  T 


sin  (Q-oo)  T-sin  co)  T 


+  i  (cos  (£2-oo)  T  +  cos  (Q+oo)t) 


(6-46) 


Again  the  imaginary  part  can  be  determined  exactly.  Making  use  of  our 
previous  result 


oo 


J  cos  aTe 
o 

2 


-  V 
2 1 


1/2 


A  00 

ImXxx  =  ^ 


v  .  ,  2,-/2 
8  (^kz  > 


exp  l-  2^ 
2 


(6-47) 


-pi- 

2PV 


(Q+oo)  “ 

+  eXP("  2(ik  2 


(6-48) 


In  the  limit  of  extra  low  frequencies,  co  «  Q  <  oo^.  This 
expression  reduces  to 


ImXxx(ELF) 


60  I 

-  V" 


CO 


(tikz2)"1//2  exp 


£2" 


2(1  k 


z  J 


(6-49) 


To  evaluate  the  real  part,  we  make  use  of  the  low  frequency  approximation 
to  expand  sin  (Q-oo)t  and  sin  (£2+oo)t  in  powers  of  cot  about  co  =  0. 

sin(&+oo)T  =  sinQT+ooT  cos  Q  T  +  ...  (6-50) 


sin(Q-co)T  =  sin&T-ooT  cos  !)  T  +  ... 


(6-51) 


Then,  from  Equation  (6-46) 
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Real  X. 


xx 


oo 

-  a.'2  f  dTe 
P  oJ 


-v-K  T 


2„2/2 


oo 


OO 


=  + 


T  cos  £2  T 

2  2/2\ 


J  dT  ^  U  *  2  J  cos  a  T 


(6-52) 


and  integrating  by  parts 

2 


CO 


°r  -)j.kz2T2/2 

Q  J  dT  sin^Te  -1 


(6-53) 


2  2 

Again,  going  to  the  low  temperature  limit,  we  assume  (ikz  «  2 
which  is  again  compatible  with  our  low-frequency  assumption  in^that  it 


th 


puts  a  less  stringent  condition  on  the  wavelength,  namely  X  >  — ^ 
or  \  >  10  cm  (for  ions)  which  is  certainly  satisfied  by  any  extra-low 
frequency  magnetohydrodynamic  disturbance  propagating  at  the  Alfven 
velocity.  Then 


co 


Real  X. 


xx 


V>K 


oo 


Q  *  dT  sinOT(l-uk„ 2T2/'2)  -1 

V  Z§ 


oo 


a* 


(6-54) 


The  real  part  of  X  can  be  determined  exactly.  From  Equation  (6-41) 

xy 


Re  L 


,s2  °o  .  2  2/2 

k  J  -  ^ 


£cos  (Q-oo)  T  -  cos  (£2+co)  tJ 


oo2 

_E_ 

0) 


/§  dkz2r1/2(exp 


(Qko)  ‘ 


2^K 


-  exp 


(Q-oo) ' 
2^kz'  J 
(6-55) 
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and  in  the  ELF  approximation 


go2& 


Re  X  -  +  - 5 — o 

^ (ELF)  (|ik  *)  ■ 


72# 


exp  - 


q* 


2^k, 


(6-56) 


2  oo 


A 

Im  X.. 


xy 


-  2co  J  d T  exp(-jakz2T2//2)  sin  (&4-oc>)  T+  sin(Q-oo)T 


(6-57) 


Again  we  expand  sin(£2-^>)T  and  sin(ft+co)T 

2  oo 


Im 


Xxy  =  J  dT  exp(-|ikz2T2//2)  sinOT 


In  the  limit  of  low  temperature  and  low  frequency  the  inequality 
2  2 

^k  «  Q  is  satisfied  and  hence 


2  oo 


Im  X 


xy 


=  -jE  J  dT  sinQT(l-(ikz2T2//2) 


2  2  2 
GO*  Llk  *G0* 

"  +  Li _ £ 

Q  GO 


(6-58) 


The  first  term  vanishes  when  summed  over  the  two  species  because  of 
the  condition  of  charge  neutrality  in  the  plasma.  Therefore 


Im  X 


,  2  2 
r  z  p 


*y  a3 


(6-59) 


oo 


Thus  ,  in  the  low- temperature ,  low-frequency  limit,  the  electric 
susceptibility  tensor  is  given  by 


zz 


CO2 

~2  +  i 
<&  Oik  *) 


go2go 

_E _ _ 

2,3/2 


exp  - 


GO 


2jik, 


(6-60) 
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A 

X. 


xx 


A 

X 


yy 


a,2 

- E 

"  G2 


+  1 


CO2 

_E 

60 


<Ji  <»*. 2»-1/2 


exp  - 


QJ 


2^K 


(6-61) 


A 

X 


xy 


A 

X 


oo2Q 


yx 


<HV > 


273/2 


fi 


exp  - 


Q“ 


2H-K 


+  1 


_  2  2 
1  z  p 

£2^0 


(6-62) 


It  should  be  noted  that  the  exponential  term  in  each  of  these 
three  expressions  is  a  very  small  quantity  because  of  the  low-temperature 
approximation 

2  2  2 

1-ik  «  a)  ,  Q  . 


In  the  zero- temperature  case,  in  the  absence  of  collisions,  X 

A  A 

and  X„  were  both  real  and  wave-number  independent,  and  X  vanished, 
z  z  7  xy 

From  the  structure  of  the  above  expressions,  one  can  easily  see  that 

A 

the  dependence  of  X(kw)  on  k  is  entirely  contained  in  the  temperature 
dependence  as  was  previously  stated. 

Again,  we  shall  first  obtain  the  dispersion  relations  from  the 
normal  modes  and  then  construct  the  temperature-dependent  Green’s 
Functions. 

The  normal  mode  equation  is  again 


E  =  0 


J\ 

k  -  — o 

i  .  n 

CM  1( 

31 

i 

A 

X 

V  ° 

1  ~  z  z 

z 

c 

i 


The  vanishing  of  the  determinant  of  the  coefficients  leads  to 


V1+Xzz> 

c 


f  2  2  \  2  /,  2  \  2 


=  0 


(6-63) 

The  resultant  modes  are  then  given  by  the  dispersion  relations 
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,  A 

<1+xzz> 


(6-64) 


m  A 

-  %  (l+X  + 
2  xx- 
c 


X  )  = 
xy 


(6-65) 


The  first  of  these  equations  represents  a  longitudinal  oscillation 
similar  to  the  plasma  oscillation  occurring  in  the  zero-temperature 
plasma. 

We  shall  restrict  our  attention  to  the  two  modes  given  by  Equation 
(6-65)  and  attempt  to  determine  the  criteria  for  thermal  damping 


In  the  limit  of  low  temperatures  and  low  frequencies 

V  2  O2 

“  «  2 


Therefore 
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X 


xx 
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A 

X 


xy 


c2 
- B 

"  a2 


|j.k  V 


(6-66) 


This  expression  is  actually  a  sum  over  the  two  species,  but 
the  values  offi,(/0p  and  a  m  both  the  P-region  and  the  exosphere  are 
such  that  the  ion  terms  predominate.  Typical  values  of  these  para¬ 
meters  in  both  regions  are  listed  below  in  Table  1. 
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co 


P- 

co 

P+- 

Q- 

Qf 

ji- 


F-Reqion 

i  -  ,n8  -1 

1 . 4.  x  10  sec 


7.0  x  105  sec"1 

7  -1 

1.0  x  10  sec 

3.0  x  102  sec"1 

„  _  ir>15  2.  2 

2. 5  x  10  cm  /sec 

1.0  x  1011  cm^/sec^ 


a  (Alfven  1.0  x  10  cm/sec 

velocity 

T (kinetic  temp) 1000°K 


Exosphere 
1.0  x  105  sec'1 


7.0  x  10^  sec'1 
1.0  x  10^  sec'1 


3.0 


sec 


-1 


r  ^  .  -15  2  ,  2 

5.0  x  10  cm  /sec 
1.0  x  1011  cm^/sec^ 

g 

1.0  x  10  cm/sec 


3000°K 


Henceforth,  the  quantities  oo^,  and  ji  appearing  in  the  dispersion 
relation  shall  be  the  ionic  contributions  only. 
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r  z  p 
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(6-67) 


2  2 

00  00  0 

p  u 

2,  2,3/2 


exp  - 


If  we  neglect  the  imaginary  part  of  Equation  (6-68)  then  to  a  first 
approximation 


2  .  .2  y  2 
_ _ 2LZa 

+  JL  “ 
1  -  2  Q 

a 


(6-68) 


Since  in  the  ELF  limit  oo  <<  Q,  the  denominator  can  be  expanded 
to  obtain 


,2  oo  /  ,  —  u.  co 

kz  **  ~2  (  1  +  2 

a  \  a 


(6-69) 
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As  ^2  is  of  or<^er  °f  kz  is  nearly  equal  to  — 2  in 

a  a 

the  first  approximation.  One  can  therefore  obtain  the  next  higher 

2 

approximation  by  replacing  k  everywhere  in  the  imaginary  part  of 

2  Z 

Equation  (6-70)  by  co_  .  Therefore 
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2 
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2 
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(6-71 ) 


and 


2 

a  ^  p 

-2-372  ~f-  exp 

C  (J,  0) 


oVV 

2pLC02/j 


(6-72) 


4 

It  can  easily  be  seen  from  this  last  equation  that  for  go  <  10  sec 
the  exponential  dominates  the  damping  term  and  tends  toward  zero.  Thus 
for  any  frequencies  in  the  ultra-low  frequency  (ULF)  range  for  which 
the  low  temperature  and  low  frequency  approximations  are  satisfied , 
the  damping  distance  becomes  very  great,  and  the  low-temperature 

4 

thermal  damping  is  negligible.  Only  for  angular  frequencies  of  10 
and  greater  does  the  imaginary  part  of  the  wave  vector,  as  given  by 
Equation  (6-72)  become  the  order  of  the  real  part.  Hence,  the  effects 
of  low-temperature  thermal  damping  may  be  neglected  for  ultra-low 
frequency  one-dimensional  propagation  along  the  magnetic  field  lines. 

The  other  mode  contained  in  Equation  (6-63)  is  analogous  to  the 

longitudinal  plasma  oscillation  resulting  in  the  zero  temperature 
w  2 

case  1  -  — § —  =  0. 

ay 

A 

Now,  in  the  low-temperature  limit,  Xzz  is  wave-number  dependent, 
and  the  plasma  oscillation  is  replaced  by 
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i 


0 


oo 

-  — 2  (1-1 

O) 


0 


(6-73) 


(6-74) 


with 


Then , 


Assume  a  solution  of  the  form  x  -  Ai 


0  _  2/2  ■ 2 

,.3,  +  A  =  -  l  A /  — 

-A  l  e  '  y  tt 


(6-75) 


(6-76) 


Therefore,  solving  for  A  by  trial  and  error,  we  find 


A  =  .95 


Thus 


CO 


v.,  k 
th  z 


=  .95  i 


(6-77) 


This  longitudinal  mode  represents  a  damped  oscillation  with  a 
damping  distance  of  100  meters,  travelling  at  the  ion  thermal  velocity 
Thus  for  the  case  of  one-dimensional  propagation  along  the  field  lines 
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the  zero- temperature  Alfven  mode  k  =  ^  is  unchanged  by  the  inclusion 

z  d 

of  temperature  in  the  ULF  limit,  whereas  the  longitudinal  plasma 
oscillation  transcends  to  a  damped  oscillation.  Again,  the  most  im¬ 
portant  feature  of  the  above  discussion  is  that  the  one-dimensional 
magnetohydrodynamic  wave  propagating  at  the  Alfven  velocity  is  unaffected, 
in  the  ULF  limit,  by  f*e  presence  of  a  finite  temperature. 
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CHAPTER  VII 

A  MANY-BODY  DERIVATION  OF  THE  BOLTZMANN  EQUATION 
FOR  HARD- SPHERE  MOLECULES 

12  3 

The  distribution  functions  f  ,  f  ,  f  • • • ,  are  connected  by  a 
chain  of  equations  which  describe  in  an  obvious  way  the  motion  of 
one  or  two  or  three  "typical  particles'1  in  the  system  relative  to 
the  rest  of  the  system.  The  first  two  equations  in  this  chain  have 
the  form 

(^T  +  v  •  V  +  a  •  V  )  f  (rvt)  +  —  f  [  F(r-r’)  •  V  f^frvr’v’t)  =  0 

ot  v  m  J  «J  v 

r 1  v* 

(7-1) 

and 

+  v  .  7  ♦  V  •  7.  ♦  a  •  7V  ♦  a'  •  V  ♦  1  p(r-r')  •  *  1  P(r'-r) 


x  f^(rvr’v't)  +  —  f  I  (Ftr-r11)  •  V  +  F(r!-r!f)  •  V  .) 

m  J  J  v  v 

r  I  1  yll 

x  f 3  (rvr'v  'r 1  !v»  ’t)  =  0  (7-2) 

The  operator  on  f  in  Equation  (7-1)  describes  the  motion  of  a  single 
particle  under  the  influence  only  of  an  external  force.  The  opera- 
tor  on  f  in  Equation  (7-2)  describes  the  motion  of  two  particles 
moving  under  the  combined  influence  of  the  external  force  and  the 
interparticle  force.  The  integrated  force  terms  are  due  to  the  com¬ 
bined  average  influence  of  more  or  less  "near”  particles,  and  is 
less  significant  at  low  densities  than  high.  The  reason  that  a 
knowledge  of  a  high  degree  of  correlation  is  needed  in  these  terms 
is  that  the  distribution  of  "near"  particles  is  not  unrelated  to 
the  same  force  which  is  being  taken  account  of. 
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Our  intuition  tells  us  that  a  large  class  of  problems  ought 

to  be  solvable  using  only  Equation  (7-1)  and  some  simple  assumption 
2 

about  f  .  Likewise  another  large  class  of  problems,  very  likely 
encompassing  the  first  class,  ought  to  require  no  more  than 
Equations  (7-1)  and  (7-2)  and  some  simple  assumption  about  f  . 
Indeed  the  ’’self-consistent  force”  approximation,  and  the  ’’Boltz¬ 
mann  collision  hypothesis”  are  the  direct  expression  of  this  in¬ 
tuition.  We  will  proceed  on  the  basis  of  a  slightly  different 
philosophy  which  has,  in  practice,  no  more  or  less  valid  basis 
than  the  aforementioned  ideas,  but  has  certain  obvious  mathematical 
extensions.  This  is  the  idea  of  successive  noncorrelations. 

When  no  internal  forces  act  (we  use  a  subscript  zero  to  indi¬ 
cate  this  case)  ,  there  can  be  no  correlations  among  the  particle 
motions.  Either  by  inspection  or  intuitively  it  follows,  for  ex¬ 
ample,  that 

f£(rv, r'v't)  =  fQ (rvt)  x  f  (r'v't)  (7-3) 

What  information  does  the  nontrivial  function 

fj  c  (rv, r'v't)  =  f  (rvt)  f  (r'v't)  (7-4) 

imply?  Let  us  find  its  equation,  and  compare  it  with  Equation  (7-2) 


(X  +  v  .  V  +  V  •  V-  ♦  .  .  Vv  +  .  Vv1>4 


+  m  J  J  [F(r-r''>  •  Vv  f  (r'v't)  f2(rv,r'  'v'  ' t) 

r 1  '  v  ’  ’ 


+  F(r'-r")  •  V  ,  f  (rvt)  f2(r'v'r'  'v'  ’t)]  =  0  (7-5) 
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2  2 

We  see  that  fNC  is  approximately  equal  to  f  when  the  two  particles 

are  far  enough  apart  that  their  mutual  force  may  be  neglected,  and 

3 

when  the  system  is  such  that  in  the  f  term  only  the  particles  con¬ 
nected  by  a  force  need  be  treated  as  correlated?  that  is, 

f3(rv,r'v'  ,r'  '  v  "  t)  =»  f  2  (r '  v '  ,r'  'v'  't)  f  (rvt)  (7-6) 

in  a  term  in  which  r1  and  r1  f  are  connected  by  a  force.  Acting  on 

the  belief  that  in  a  multiparticle  system  with  two-body  forces  the 

2 

two-particle  function  f  will  be  of  greater  numerical  importance 

3 

than  f  ,  no  matter  how  strong  the  forces  are,  we  conclude  that  the 

approximation  expressed  by  (7-6)  will  be  of  extreme  generality,  and 

2 

especially  even  so  when  the  arguments  r,  r'  of  f  are  close  to  one 

2  2 
another.  But  fNC  is  then  a  bad  approximation  to  f  .  Consequently, 

Equations  (7-1)  and  (7-2)  together  with  assumption  (7-6)  in  the 
force  term  of  (7-2)  provide  a  reasonable  way  to  terminate  the  other¬ 
wise  endless  set  of  coupled  equations. 

We  can  now  derive  an  equation  for  the  quantity 

D(rv, r'v't)  =  f2(rv  r'v't)  -  f  (rvt)  f  (r'v't)  (7-7) 


Let  us  first  introduce  a  simplifying  notation.  Let  us  number  the 
particles  1,  2,  3,  etc.,  and  let  the  number  stand  for  the  full  set 
of  position-velocity  variables.  For  example,  1  =  (^v^) ,  2  =  (z^v^ 
3  =  (r^v^) • • • .  The  variable  t  is  the  same  in  all  functions  and  may 
be  ignored.  The  operator 


L  = 


at  +  vi 


vi  +  ai 


+  V, 


V2  +  a2 


(7-8) 


may  be  called  L(12«*0  .  It  describes  the  linear,  noncorrelated 
part  of  the  motion  of  the  "typical  particles"  1,  2,*-*.  The  oper¬ 
ator  M(12...)  will  describe  the  mutual  forces.  As  examples,  M(l)  =  0, 
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M (12)  =  Jj  F(r,-r,)  •  V  +  i  F(r,-r,)  •  V 
m  12  v^  m  2  1 


M(123)  =  J(F(r.-r,)  +  F(r,-r,))  •  V  +  £  (  )  •  V  +  J  (  )  •  V 

m  1  2  13  V1  m  v2  m  V3 


(7-9) 


with  the  variable  in  the  force  terms  appearing  in  cyclic  order. 

1  2 

The  superscripts  on  f ,  f  may  be  left  off  where  the  arguments 

of  the  functions  are  given,  for  then  the  "order”  of  the  distribution 
function  is  obvious.  Finally,  we  may  let 


1 

m 


P(rx-r2) 


A  (1-2) 


(7-10) 


so  that,  for  example, 

M(123)  =  A(12)  +  A  (13)  +  A(21)  +  A(23)  +  A(31)  +  A(32)  (7-11) 

In  this  notation,  the  two  equations  we  intend  to  solve  are 

L  (1)  f  ( 1)  +  j  A(12)  f  (12)  =  0  (7-12) 

2 


and 


[L(12)  +  M (12)]  f  (12)  +  J  (A(13)  f  (13)  f  (2)  +  A(23)  f  (23)  f  (1) ) 

3 


0 

(7-13) 


Let 


D  (12)  =  f  (12)  -  f  (1)  f  (2) 


(7-14) 
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Since 


L  (12)  (f  (1)  f  (2) )  +  J  (A (13)  f  (13)  f  (2)  +  A(23)  f  (23)  f  (1) )  =  0  (7-15) 

3 


we  find 


[l  (12)  +  M(12)]  D(12)  =  -M(12)  (f  (1)  f  (2))  (7-16) 

The  function  D(12)  may  be  introduced  into  the  equation  for  f(l)  in¬ 
stead  of  f (12) .  Before  we  do  so,  though,  we  can  increase  the  sym¬ 
metry  of  our  equations  by  replacing  A(12)  by  M(12)  in  Equation  (7-12) . 
Note  that 

M  (12)  =  A  (12)  +  A(  21)  (7-17) 

The  A (21)  term  gives  zero  under  the  integral  ,  for  this  contains 

T  2  y  2 

jF(r!-r)*  V  ,f  (w')  which  integrates  exactly f  leaving  f  (w1) 

evaluated  at  the  limits  of  the  v'  integrals.  In  any  physical  system 
there  will  be  no  particles  of  infinite  velocity,  so  that  the  inte¬ 
grated  term  is  zero.  Thus,  Equation  (7-12)  can  be  rewritten: 

L  (1)  £  (1)  +  J  M  (12)  f  (12)  =  0  (7-18) 

2 

and  the  substitution  of  (7-14)  gives 


L  (1)  +  J  M(12)  f  (2)  f(l)  +  J  M(12)D(12)  =  0  (7-19) 

2  -I  2 

The  operator  added  to  L(l)  arises  (of  course)  in  the  noncorrelation 
approximation.  Written  out  it  is 
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J  M  (12)  f  (2)  =  J  A(12)  f  (2)  =  Jj  J  J  F(r-r')  f  (r'v't)  •  V, ‘ 

2  2  r '  v ' 

It  is  often  referred  to  as  the  "average-f orce'1  term,  and  was  the 
basis  of  Vlasov’s  and  Landau’s  treatment  of  the  long-range  forces 
in  a  plasma.  As  we  mentioned  f^  is  (together  with  (7-6))  a  good 
approximation  to  when  the  particles  are  ’’far  enough  apart." 

Since  the  Coulomb  forces  do  not  have  a  finite  range,  it  is  diffi¬ 
cult  to  define  for  them  the  precise  region  of  validity  of  the  non- 
correlation  approximation.  It  proves  to  be  applicable  to  the  plasma, 
nevertheless,  because  of  the  plasma's  ability  to  shield  charges. 

The  effective  forces  in  the  plasma  are  of  finite  range,  so  that 
the  solutions  using  the  noncorrelative  approximation  are  consistent 
with  the  approximation.  But  when  the  forces  are  solely  or  primarily 
short-range,  this  term  has  contributions  only  from  one  or  two 
neighboring  particles?  in  a  distribution  theory  such  a  constribution 
is  of  negligable  weight.  The  last  term  of  Equation  (7-19)  then  is 
most  important. 

Our  main  problem  now  is  to  solve  (7-16)  for  D.  The  equation  is 
inhomogeneous,  and  may  be  solved  by  a  Green’s  function  for  the 
particular  solution.  The  homogeneous  solution  will  have  an  arbi¬ 
trary  constant  multiplying  it  which  must  be  adjusted  to  the  physi¬ 
cal  problem  ("the  boundary  conditions") .  In  fact,  there  cannot  be 
an  homogeneous  solution  --  the  constant  is  zero  —  because  only 
that  part  of  the  force  on  a  particle  [f(l)]  must  be  included  which 
is  specifically  due  to  the  self-consistently  positioned  neighboring 
particles  that  are  within  the  range  of  the  force.  An  homogeneous 
term  D°  would  make  no  reference  to  this  self-consistency,  but  the 
inhomogeneous  solution  of  Equation  (7-16)  would,  by  virtue  of 
the  source  term  on  the  right  hand  side. 

For  the  construction  of  the  Green's  function  we  shall  need  to 
indicate  the  time  variable  t  at  which  the  particle  variables  are 
specified.  Let  G(12t,  34tQ)  be  the  solution  of 
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(L  (12)  +  M(12))G(12t,34t0)  =  6  (t-tQ)  &  (1-3)  6  (2-4)  (7-20) 


In  the  usual  pictorial  manner  of  speaking,  G  is  said  to  take  particles 
3  and  4  at  time  t  into  particles  1  and  2  at  time  t.  That  is,  G  des¬ 
cribes  the  details  of  the  orbit  and  momentum  exchange  of  two  particles 
under  the  influence  of  the  inter-particle  force.  In  full  detail,  but 
leaving  out  the  external  acceleration  A,  we  have 


V1  +  V2 


1 

m 


F(r, 


"r2) 


1 

m 


F(r2-r1) 


V 

v 


) 

2 


x  G(rivir2v2t?  r3v3r4v4t0)  =  6  6  3  (rl-r3)  6  {vl_v35  6  lr2-r4) 


x  6  3 ^v2-v4^  (7-21) 

All  of  the  delta-functions,  with  the  exception  of  the  first,  are 
three-dimensional,  e.g. , 

&3(rl"r3)  =  Mxrx3)6  (yj-y3)& (zi"z3)  (7-22) 


We  write  the  solution  as 


D (12t)  -  -  J  J  J  G(12t,34t0)M(34)f(3t0)f(4t0))  (7-23) 

3  4  fco 

Application  of  L(12)  +  M(12)  to  this  equation  leads  to  Equation  (7-16) 
in  view  of  Equation  (7-20) . 

There  remains  one  element  of  uncertainty  in  G  —  it  may  be 
"retarded"  or  "advanced"  or  a  mixture  of  the  two.  Looking  at 
Equation  (7-23)  we  realize  that  D(12t)  should  depend  only  on  the 
earlier  particle  distributions,  since  all  of  the  dynamical  processes 
are  causal.  This  means  that  G(t-tQ)  should  be  "retarded",  i.e. , 

G  -  0  if  t  <  t  . 
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Our  interest  in  the  solution  of  the  two  coupled  equations  stems, 
in  the  present  instance,  from  a  desire  to  treat  collisions  more 
accurately  than  is  possible  with  the  Fokker-Planck  operator.  We  can, 
for  this  purpose,  try  to  solve  Equation  (7-21)  without  the  external 
acceleration  a,  which  is  almost  never  as  strong  as  the  internal  ac¬ 
celerations  jjj  F(r-r’)  ,  and  acts  mainly  to  provide  a  net  drift  super¬ 
imposed  on  the  "chaotic"  coll isional-dif fusion  processes.  Furthermore, 

the  quantity  D  gives  the  "closer-in"  effects  when  contrasted  with  the 

2 

"long-range"  effects  of  fNC*  Hence,  as  a  kind  of  orientation  in  the 
problem  of  collisions,  and  also  for  an  improved  model  for  collisions 
in  the  plasma,  we  treat  the  problem  of  hard  spheres  colliding.  This 
is  certainly  an  improvement  over  the  Fokker-Planck  procedure  which 
applies,  as  we  have  noted  before,  to  the  motion  of  large  molecules 
suffering  many  small  impacts  in  the  surrounding  medium  of  small 
molecules.  Whether  the  collisions  are  among  electrons,  among  ions, 
between  ions  and  electrons,  or  between  charged-particles  and  neutral 
molecules,  we  expect  some  improvement  in  description. 

Let  us  carry  out  the  analysis  as  far  as  possible  without  special¬ 
izing  the  force.  Hard- sphere  collisions  will  prove  to  simplify  the 
final  analysis.  Since  we  know  that  in  a  two-particle  collision  total 
momentum  is  concerned,  and  since  the  force  lies  along  the  line  between 
the  centers  of  the  (pt.)  molecules,  we  shall  introduce  appropriate 
variables.  It  is  enough,  for  the  present,  to  treat  one  kind  of  mole¬ 
cule  to  illustrate  our  ideas. 

The  center  of  mass  (R) ,  the  relative  distance  (r) ,  the  center 
of  momentum,  V,  and  the  relative  velocity  (v)  of  particles  1  and  2 
are 


R  =  2  (rl  +  r2> 


r  =  rl  -  r2 


V  = 


V1  +  V2 


v  =  (vx  -  v2) 


(7-24) 
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A  similar  set  of  variables ,  constructed  from  particles  3  and  4  will 
be  called  RQVorovo.  It  is  easy  to  verify  that  the  gradients  trans¬ 
form  as  follows 


n  _  ±  y  +  y 
Vr1  2  VR  +  vr 


V  =  “  V__  +  V 
V1  2  V  v 


V  =  —  v  _  v 
r2  2  VR  r 


7  =  ~  V  -  V 

v2  2  V  V 


(7-25) 


so  that 


v.  •  V  +  v0  •  V  =  v  •  V  +  V  •  V 
■L  r  2  r  2  r  R 


(7-26) 


and 


(7-27) 


Furthermore,  since  the  Jacobian  of  the  transformation  (7-24)  is 
unity,  the  delta-function  term  in  Equation  (7-21)  may  be  rewritten 
directly  in  terms  of  the  new  variables.  The  resulting  equation 
makes  use  of  Newton’s  third  law: 

F(r2  -  r1)  =  -F(r1  -  r2>  (7-28) 


and  becomes 

+  v  •  VR  ♦  v  '  *  m  F(r)  '  V  <3<RVrvt,R0V0r0V0t0) 

=  b  (t-tQ)  6(R-RQ)  6  (V-VQ)  6  (r-rQ)  6  (v-vQ)  (7-29) 

We  have  used  the  same  symbol  for  the  Green's  function  written  in  terms 
of  the  new  variables.  We  shall  speak  quite  generally  of  Equation  (7-29) 


113 


S-2023-1 


as  describing  a  collision,  although  it  will  describe  a  bound  system 
also  when  the  forces  can  support  one.  Inspection  of  Equation  (7-20) 
shows  directly  that  the  total  ’'momentum’’  V  is  unchanged  during  the 
collision.  Hence  G  contains  b  (V-VQ)  as  a  factor.  Further,  we  can 
also  remove  a  factor  expressing  the  fact  that  the  center  of  mass 
moves  in  a  straight  line: 

G  =  6 (V-V  )  6 (R-R  -  V  (t-t  )  ) G (rvt ,r  v  t  )  (7-30) 

o  o  o  o  7ooo 


where 


[wt  *  v  •  Vr  +  i  F<ri  •  Vv] 

=  b  (r-rQ)  6  (v-vq)  &  (t-tQ)  (7-31) 

For  a  force  which  can  be  derived  from  a  central  potential,  X: 

F  (r)  =  -  VX(  (r  1)  =  -  ~j-  X  •  V  |r|  =  -  X<  r  (7-32) 

in  which  r  is  the  unit  vector  along  the  radial  direction  (from  the 
origin  of  coordinates)  .  Consequently,  a  spherical  coordinate  system 
for  r  together  with  a  suitable  chosen  velocity  coordinate  system  at 
each  point  of  space  will  simplify  Equation  (7-31)  in  such  a  way  as 
to  make  clear  the  fact  that  the  force  between  the  two  particles  only 
causes  a  radial  acceleration.  The  transformations  for  the  coordi¬ 
nates  are 


x  =  |r|  sin  0  cos  $5 

y  =s  |r|  sin  0  sin  ^  (7-33) 

2  =  |r|  cos  0 

2  -1 

whose  Jacobian  is  (r  sin  0) 
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The  transformation  of  velocities  has  to  depend  on  that  for 
coordinates.  For  if  xyz  are  imagined  to  be  functions  of  time,  then 
so  are  r0g$,  and  the  velocities  vx  v^,  vz  =  x  y  z  would  be  related  to 
the  time  derivatives  r  6  jz!  by  means  of  the  following  equations  (ob¬ 
tained  by  differentiating  Equation  (7-33) ) : 


vx  =  r  sin  0  cos  f6  +  Or  cos  0  cos  &  -  jH  r  sin  0  sin  fi 

Vy  =  r  sin  0  sin  fi  +  6  r  cos  0  sin  ^  +  &  r  sin  0  cos  f6 

vz  =  r  cos  0  -  6  r  sin  Q  (7-34) 

2  -1 

whose  Jacobian  is  also  (r  sin  Q)  “  .  In  the  present  formulation  — 

which  might  be  called  the  hydrodynamical  treatment  of  the  two-body 

problem  —  the  quantities  r  6  0  have  to  be  considered  as  independent 

variables.  Call  them  v  v~  v,.  We  find 

r  0  p 


r  =  v  cos  sin  0  +  v  sin  &  sin  0  +  v  cos  0 

r  x  y  z 

rd  -  b  (v  cos  cos  ©  +  v  sin  f6  cos  0  -  v  sin  0) 
y  i  x  y  z 


v  —  - 

p  r  sin  0 


(-v^  sin  f6  +  Vy  cos  tf) 


(7-35) 


Now,  G(xyz,  vxvyv2)  becomes  a  new  function  G(r0^vrVgV^)  .  If  we 
recall  the  rules 


3x  ^(r(x)  ...vr(Xlvx)  -  ..)  =  |f  |?  + 


a  a  v  A 
0r  0G  .  _r  5G__ 

5x  9v 


(7-36a) 


and 


a  v. 


+  •  *  • 


(7- 36b) 
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we  shall  have  no  trouble  converting  Equation  (7-31)  to  its  new  forms 
[dt  +  vr  £  +  ve  QQ  +  h  +  (rv0  +  r  sin2  0  -  m  X,)  9^7 


2Vr-VQ  9 

+  (-  -  +  vt  sin  0  cos  0)  r — 

r  p  o  v 


0 


2v  V/ 

+  <- 


r  -  2  vov^  cot 


&  (t-tp)  6  (r-rQ)  &  (Q-Qq)  &  (gMQ)  6  (Vp-Vp0)  &  (vq-vq°)  5 

(r2  sin  0  ) 2 
o  o 


(7-37) 


From  angular  momentum  conservation  we  can  deduce  three  variables 
and  that  greatly  simplify  Equation  (7-37) .  Converting  r  x  v  =  L 
to  our  new  coordinates  gives 


2 

Lx  =  -r  (Vg  sin  0  +  v^  sin  0  cos  0  cos  0) 

2 

Ly  =  r  (Vg  cos  -  v^  sin  0  cos  0  sin  0) 

L  =  r^  v v  sin^  0  (7-38) 

Z  0 


We  can  use  L  L  L  to  eliminate  vn  v  ,  and  Q  not  0  because 

a  x  y  Q  P  L 

Lx=  Ly  and  Ly  =  -  5^  LXJ  *  Any  function  of  Lx  or  Ly  or  Lz  satis¬ 
fies  the  homogeneous  part  of  (7-37) .  We  therefore  choose  a  factor 
of  G  that  will  reduce  to  part  of  the  inhomogeneous  term  when  the 
remaining  variables  are  correctly  given.  In  fact,  we  choose 


G(rO t  vrv0v^)  =  6  (Lx-Lx°)  6  (Ly-Ly°)  6  <  VLz°> 


x  G(r jrfvr  LxLyLz) 


(7-39) 
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General  formulas  of  use  are 


L2  =  L2  +  L2  +  L 2  =  r4  (v2  +  v2  sin2  0)  ; 


-  Lx  sin  /  r  L  cos  *$  =  r  vQ  v  .  =  — - = — 

y  ^  sin  Q 


1?  -  (-L  sin  ^  +  L  cos  rf)  ^  sin^  0  =  - rr 

x  y  P 


sin  0 


From  these  we  conclude  that 


v*$  -  [Lz  +  <LX  cos  t  +  Ly  sin  *$)  2] 


(7-40) 


With  the  new  choice  of  variables, 
2 


JL  ,  /Lr  2X1\  JL  ,  JL 

3t  ^3  m  '  3v^  Vr  dr 


2  y  2 
Lz  +  ^Lx  COS  ^  +  Ly  sin  ^ 

r2L 


3?]  5(r*V 


=  r2  sin2  jz$q  |cos  20q|  |v^  j  6  ( t-tQ)  6  (r-rQ)  &  (vr-vr°)  6  (*$-*$0) 


(7-41) 

The  appearance  of  |v^°|  requires  one  comment.  We  see  from 

Lz  =  L»z°  that  the  sign  of  is  a  constant  —  i.e.,  the  particle 

always  continues  in  the  same  angular  direction  around  the  center 

of  force.  Th„  ©dependence  could  be  greatly  simplified  because,  as 

is  well  known,  the  motion  of  the  particle  is  in  a  plane  fixed  in 

space,  determined  only  by  the  angular  momenta.  In  fact,  if  Vg  =  0 

and  ©  =  tt/2  initially,  then  vQ  =  0  and  ©  =  tt/2  for  all  times  i.e., 

^  2 

L  =  L  =  0,  and  the  coefficient  of  3/9*$  is  simply  L  /r  ,  the 
x  y  z 

angular  velocity  in  the  plane. 
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Equation  (7-41)  can  be  simplified  in  its  r-dependence  by  use  of 
the  constantcy  of  energy.  Apart  from  a  factor  of  m/2,  the  reduced  mass 
the  energy  is 


E  = 


1  L_ 

2  r2 


2X(r) 

m 


(7-42) 


and  any  function  of  E  satisfies  the  inhomogeneous  part  of  (7-41)  , 
Using 


we  readily  verify  that 


6 


(7-43) 


(7-44) 


solves  (7-41)  .  If  the  actual  discontinuity  comes  from  a  factor  H(t-t0) 
that  insures  that  G  is  retarded  --  then  (7-44)  gives  6 (r-rQ)  exactly 
as  written.  Then 


& 


(7-45) 


in  turn  gives 

6  (vr-vr0) 

exactly.  Of  course  the  integral  in  Equation  (7-44)  must  be  taken 
over  the  complete  orbit  between  rQtQ  and  rt,  so  that  vr  in 
Equation  (7-43)  must  be  understood  to  be  the  positive  root  when  r 
is  increasing,  and  the  negative  root  when  r  is  decreasing.  Inter¬ 
estingly,  the  |v  °|  of  (7-44)  and  (7-45)  cancel.  All  that  is 
needed  to  complete  the  solution  is  a  6 -function  of  the  form 
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6 


*  dr' 

J  g(r') 


(7-46) 


Lro  *o  J 

which  satisfies  the  homogeneous  part  of  Equation  (7-41)  ,  and  reduces 
to  the  remaining  factors  of  the  inhomogeneous  part  namely 

ro  sin2  ^o  |cos  20ol  | V/l  6  ^  ’  *o} 

when  r  =  r  .  Clearly  (7-46)  becomes 
o 


|f  (*So)|  6  (>z5-^o)  (7-47) 

in  this  case.  Now  (7-41)  is  readily  solved  by  quadratures,  and  we 
find  that 


(again  with  the  appropriate  choice  of  sign)  while 
f(0)  =  |l2  +  (Lx  cos  f6  +  Ly  sin  f6)  2J 

Using  (7-38)  we  find  that 

f  W  = 

so  that  the  final  result  is 

G  =  H (t-t  )  & (L  -L  °)  6 (L  -L  °) 6 (L  -L  °) & (E-E  )  x  6 (t-t 

U  X  X  Jl  Jl  z  z  u  u 

x  6(  J  fifr  -  J  7~ff  ■>  sin2  lcos  20ol 

M  r  'r  '  r' 

o 


(7-48) 


(7-49) 


(7-50) 


(7-51) 


o 
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where 


H(x)  =  1 


=  0 


if  x  >  0 


if  x  <  0 


and 


H  (x)  =  b  (x) 

This  remarkable  formula  is  correct  for  any  two-body  central  force 
problem.  Let  us  use  it  to  derive  the  Boltzmann  Equation  for  a 
hard-sphere  gas. 

In  a  collision  between  "hard"  spheres  the  interaction  acts 
only  for  one  instant  of  time,  and  the  pair  of  spheres  does  not  move 
any  distance  during  this  instant.  Actually,  no  atoms  are  perfectly 
hard,  but  instead  may  be  visualized  as  having  a  core  of  increasing 
hardness.  The  only  effect  of  a  collision  of  hard  spheres  will  be 
to  reverse  the  relative  radial  velocity,  the  reversal  being  the 
more  sudden  the  "harder"  the  spheres.  When  the  centers  of  the 
spheres,  of  radius  a,  are  more  than  2a  apart,  the  potential  between 
them  in  zero.  At  r  =  2a  the  potential  is  repulsive  (positive)  and 
rises  steeply  to  a  large  value.  The  simplest  analytical  representa¬ 
tion  which  captures  the  physics  of  the  collision  process  is  as 
follows: 


X(r)  =  (2a-r)K  r  <  2a 

=  0  r  >  2a  (7-52) 

Hence,  the  radial  force  is  a  constant  (large I) : 
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F  =  -X'  (r)  =  K  r  <  2a 
=  0  r  >  2a 

=  K  H(2a-r)  (7-53) 

where 


K  >  0 


Let  us  write  out  Equation  (7-19)  with  Equation  (7-23)  substituted 
to  see  how  this  choice  for  F  simplifies  the  equations 

(^  +  v  •  V)1  f  (it) 


5  J  J  I  J  p(ri-r2>-(V  -  Vv  >  °ret<12t’  34to> 

m  .  1  2 


t  2  3  4 
o 


X  F(r3-r4)  •  (Vv  .  V..  )  2  (3*. )  f  (4tQ) 

-  m  J  P(rrr2)-(Vv  -7  )f(2t)f(lt) 


1  2 


4K 


m 


J  J  J  J  J  J  H(2-r>  H(2a-r0) 


2  R  V  v  r  t 
o  o  o  o  o 


(7-54) 


9v  ret  5v 
r  r 


ff  -  f  j  «<2a-r> 

2 


dv„ 


ff 


(7-55) 


Making  use  of  the  fact  that  K  is  very  large  (compared  to  other 
forces,  e.g.,  the  centrifugal  forces)  we  can  simplify  Equation  (7-51). 
In  particular,  in  6(t-**0  let  us  replace  dr  by  dvr*  Using  B*Eq, 
we  find  that 


12] 


S-2023-1 


dv 


m 

r  2K 


when  X  is  sufficiently  large.  Hence, 


Mt-to  -  J  2e,Mb-eo) 
r 


m 

2K 


Mt-to  -  2K(vr-vr°)Mr-rc 


_  JS,  (v2  _ 
4X  v  r 


(7-56) 

Because  of  the  first  factor  of  (7-56) ,  we  can  replace  H(t-tQ)  in 
(7-51)  by  H(vr-vr°)  .  Equation  (7-56)  shows  that  as  X-*oo  ,  t-*tQ, 
r-*rQ,  so  that  in  (7-51)  ,  Vg->Vg°,  v^-*v^°  and  the  net  con¬ 

stant  is  unity 1 

Upon  substituting  G  into  Equation  (7-55) ,  the  term  resulting 

from  9/3v  H(v  v  ®)  =  6 (v  -v_  )  exactly  cancels  the  2nd  term  on 
r  r  r  r  r 

the  right-hand  side  of  (7-46)?  The  remaining  terms  must  be  evaluated 
in  the  X-*oo  limit.  Several  of  the  integrals  may  be  performed  di¬ 
rectly,  leaving  the  structure 

2a  oo 

1  "  ?  J  J  J  dro  J  dvrQ  H(vr“vr) 

2  t  o  -oo  °  ° 

o 

*  ‘[‘-‘o  -  at  (Vvr0>]}  3^-  £t  [ro  vrQ 

o 

(7-57) 


NOW 


J  ar0  6  [r-ro  -  «  (v?  -  vr0>]  -  1 

O 
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if  r  -  *  <vr 

write 


is  less  than  2a,  =  0  otherwise. 


Hence,  we  can 


H 


r  + 


m 

4K 


for  it.  Performing  also  the  t  integral  leads  us  to 


I  = 


2  -oo 


x  H 


r  + 


ff 


J2-  (v2 
4K  KVr 


-  vj  )  ,Orf;  vr  vQv0,  V,R  +  V  ^(vr-vr  ),t^(vr-vr 


m 


(7-58) 


Now  in  this  equation  d/dvr  must  be  understood  to  act  only  on  the  vr 
dependence  that  existed  in  ff  before  we  did  the  t  integral?  but 
3/0v acts  on  the  vr  dependence  in  H  and  ff.  To  complete  the  deriva¬ 
tion,  we  have  to  note  that  the  K-*oo  limit  may  be  performed  in  ff 
before  3/Qvr  or  3/3vrQ  are  applied,  but  not  in  H.  The  reason  for 
this  is  that  ff  will  be  a  smoothly  varying  function  of  its  arguments, 
whereas  H's  properties  change  discontinuously  when  it  loses  its 
v  dependence!  So,  set  K  =  oo  in  ff ,  and  evaluate 


?  £ -  -y  -  isbr  • 


H 


4K  r 


-  vr  6  [2a-r+«  {vr  -  vr ■>] 


(7-59) 
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Since  r  <  2a ,  (7-59)  requires  the  vr  integral  to  be  restricted  by 
2  2  ° 

the  condition  vrQ  >  vr  *  Subject  to  this  restriction,  however,  (7-59) 
may  be  replaced  by  its  limiting  value 

vrb  (r-2a) 


Since  H(vr  -  vrQ)  implies  vf  >  vrQ,  the  vrQ  integral  extends  to  vr 

when  v^  is  negative,  and  -v  when  vr  is  positive.  Thus,  in  either 

case  of  the  v^  integral 
o 


<£  + 


V, 


V f  (rivit) 


=  J  J  vr  6  (r-2a)  ff  £rOjrf,  -  vrvQv^,  V,R,tJ  (7-60) 


r 2  V2 


where 


V  = 


Vl+V2 


v  -  vi  v2 


R  = 


rl+r2 


r  “  rl_r2 


and 


ff  =  f(R  +  h,  V  +  y)  f(R  -  i  r,  V  - 


— ) 
2  ' 


with 


V*  =  -  |vr|  vev^ 


(7-61) 


The  collision  term  has  a  clear  meaning.  The  vr  >  0  and  vr  <  0 
contributions  determine  respectively  the  "IN"  and  "OUT"  contributions. 
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The  rate  at  which  particles  are  scattered  out  of  the  r^v^  component 
of  f  involves  the  velocities  v^  and  before  collision;  but  for 
particles  scattered  into  the  r^v^  component  depends  on  the  velocities 
v^  and  v^  that  must  exist  before  collision  in  order  that  v^  and  V2 
occur  afterwards!  Now  R  +  1/2  r  =  r^  and  R  -  1/2  r  =  r2  ,  and 


v  +  T  -  vi 


v  -  v  =  v. 


when  v  <0 
r 


V  +  ~r~  —  v,  -  e  (e  *v) 


V  *  “2"  “  v2  +  e(e-v) 


when  vf  >  0 


where  e  is  the  unit  vector  joining  the  centers  of  the  colliding  hard 
spheres . 

One  can  easily  verify  from  these  equations  that  the  conservation 
of  mass  flow  still  holds.  The  methods  which  are  usually  applied  to 
the  Boltzmann  equation  are  thoroughly  described  in  Chapman  and 
Cowling’s  work  on  ’’The  Mathematical  Theory  of  Nonuniform  Gases” 
(Cambridge  Press,  second  edition,  1952) .  We  shall  consider  an  at¬ 
tempt  at  deriving  directly  from  the  linearized  equation  the  dis¬ 
persion  relation  describing  the  possible  modes  of  the  system.  At 
a  later  date  the  same  procedure  will  be  considered  for  the  more 
interesting  problem  of  the  plasma. 
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CHAPTER  VIII 

LINEARIZATION  OF  THE  BOLTZMANN  EQUATION  FOR  HARD  SPHERES 

Suppose  that  a  gas  of  hard- sphere- like  particles,  at  a  finite 
temperature,  is  sharply  disturbed  so  that  a  small  group  of  particles 
is  set  moving  in  a  definite  direction  with  high  speeds 9  In  fact, 
suppose  that  a  signal  is  communicated  to  the  gas  by  varying  the  amp¬ 
litude  of  the  disturbances  with  time.  The  particles  so  excited  repre¬ 
sent  a  certain  amount  of  momentum  and  energy  injected  into  the  gas 
which,  if  the  gas  is  isolated,  will  be  conserved  throughout  the  future 
history  of  the  whole  gas.  When  the  density  is  low,  so  that  the 
mean-f ree-path  is  considerable,  the  disturbance-signal  will  be  found 
across  the  gas  with  very  little  change.  When,  however,  collisions 
come  into  play  quickly,  the  energy  and  momentum  gets  shared  over 
larger  and  larger  groups  of  particles.  The  limiting  situation  may 
be  described  as  follows:  the  final  motion  of  the  gas  as  a  whole 
resembles  the  gas  motion  before  disturbance,  but  with  a  slightly 
higher  average  energy  per  particle  (a  higher  gas  temperature) -  The 
average  momentum  is  lost  to  the  vessel  containing  the  gas.  No  sig¬ 
nal  is  received  across  the  gas  if  this  complete  thermal ization  takes 
place  in  times  less  than  the  least  time  possible  for  a  remnant  of  the 
original  signal  to  get  through  in  any  recognizable  shape. 

Now,  except  for  the  very  low  density  case,  the  original  signal 
will  be  thermalized  locally;  that  is,  the  impulse  of  energy  and 
momentum  rapidly  will  be  distributed  to  a  larger  group  of  particles, 
though  not  to  the  whole  gas.  In  this  state  of  its  evolution,  the 
impulse  will  consist  of  a  small  region  of  the  gas  of  slightly 
higher  temperature  than  the  rest  of  the  gas,  and  of  small  nonzero 
average  momentum  compared  to  the  rest  of  the  gas.  Whether  any  part 
of  the  signal  will  be  detected  after  this  stage  has  been  reached 
can  be  studied  by  treating  solutions  of  the  gas  equation  which 
differ  only  slightly  from  the  equilibrium  solution.  This  approach 
permits  us  to  linearize  the  gas  equation  (in  this  case,  Boltzmann’s 
equation)  and  simplifies  the  ensuing  analysis. 
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That  such  a  treatment  must  —  at  a  finite  temperature  —  lead 
to  the  possibility  of  transmission  of  signals  is  known  from  direct 
experience  with  the  phenomena  of  sound,  and  also  from  the  possibility 
of  deriving  the  sound  velocity  from  the  linearized  form  of  Euler's 
Equations,  together  with  the  thermodynamic  equation  of  state  for 
an  ideal  gas,  By  an  "Ideal"  gas  we  mean  one  without  any  important 
interaction  except  that  its  particles  must  be  capable  of  interchang¬ 
ing  energy  and  momentum  to  bring  about  an  equilibrium  velocity  dis¬ 
tribution.  (Our  hard-sphere  gas  approximates  the  ideal  gas  in  a 
reasonable  way.)  The  actual  calculation  is  simple.  The  time  deriva¬ 
tive  of  the  mass  continuity  equation  reads 


(8-1) 


The  linearized  Euler  Eq.  relates  3/5 t  j  to  the  pressure,  p. 


(8-2) 


The  equilibrium  equation  of  state,  for  an  ideal  mono-atomic  gas,  is 


P  = 


(8-3) 


But  adiabatic  variations  of  the  pressure  follow  the  familiar  law 


&E  .  JB.  Ml/v)  _  5  fle. 
p  cy  (l/v)  3  p 


Thus,  the  density  fluctuations  are  given  by 


5  KT 
3  m 


^28 


(8-4) 
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which  has  the  dispersion  relation 


co 


2 


s 


k 


2 


(8-5) 


with 


-VF 


KT 

m 


(8-6) 


v  is  also  one  definition  of  the  average  thermal  velocity.  The  prob- 
lem  we  set  for  ourselves  is  to  derive  the  dispersion  relations  which 
follow  from  the  hard-sphere  equation?  to  investigate  the  way  parti¬ 
cles  of  finite  size  contribute  damping  terms  to  equations  like  Eq,  (8-5)  . 
This  means,  in  effect,  deriving  the  thermodynamic  equations  of  state 
from  first  principles.  That  this  can  be  done  to  some  extent  was 

shown  by  an  analysis  of  Wang,  Chang,  and  Uhlenbeck,  an  outline  of 

( 8) 

the  method  being  given  in  lecture  notes  of  Uhlenbeck.  We  shall 
take  that  method  as  a  guideline  for  our  researches  into  the  micro¬ 
scopic  derivation  of  macroscopic  wave-motions.  In  the  present  Chap¬ 
ter  we  shall  give  a  brief  treatment  of  it. 

We  have  derived  the  Boltzmann  Eq.  in  the  preceding  section. 

We  transcribe  it  here  in  a  new  form  by  letting 


and 


e  = 


rl  I  r2 

lrl  -  r2l 


v  = 


(rl  -  r2> 

lrl  -  r2 I 


(V,  - 


v2> 


=  e 


(V1  -  V2} 


v'  = 


V1  ‘  vr  e 


vi  = 


v2  +  vr  e 
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The  result  is 


+  v1  •  V1)f(r1v1t) 

=  J  J  vr  6(lrl  ~  r2 1  "  2a) 

r2  v2 

f  (r1v1t)  f  (r2v2t) 

v  <0 
r 

f(r1v1't)f(r2v2't) 

vr  >  0 

(8-7) 

We  shall  take  the  step  of  making  the  hard-sphere  gas  into  an 
ideal  gas.  This  can  be  done  by  imagining  the  radius  a  of  the  spheres 
to  be  so  small  that  we  can  replace  r2  by  r^  in  f  because  of  &(  |ri”r2 
The  direction  of  the  vector  -  r2,  or  of  is  still  quite  important, 
however . 
of  r„  and  r. 


After  setting  in  f(r2)  r2  ^  we  can  then  introduce  instead 


ot  r^  and  r2  the  variables  r  =  and  r  =  r^  -  r 
r  gives  (2a) 2  and  an  angular  integral  JdQe 


2«  The  integral  over 
over  the  directions  of  e. 


This  leads  to 


V)  f  (rv1t) 


(2a)  2 


d&  (dvjv 

0  C  J  4  I 


f  (rv-^t)  f  (rv2t) 
f (rv^t) f (rv^t) 


vr<0 

vr>0 


(8-8) 


The  cumbersome  restriction  on  vr  can  be  removed  now  by  replac¬ 
ing  vr  by  ±.  where  indicated,  and  extending  the  integrals  over 
the  whole  range  with  an  extra  factor  of  1/2 . 


(^  +  v1  •  V)  f(rv1t)  =  1/2  (2a)  2  J  dQ  J  dv2  |  e  •  (v^-v^l  jf  (rv£t)  f  (rv^t) 


f  (rv^^t)  f  (rv2t) 


} 


(8-9) 
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We  linearize  Eq.  (8-9)  by  setting 


f(rvt)  =  n  f  (v)  (1  +  g(rvt)) 
o  o 


where  nQ  is  the  density  of  particles  N/V,  and  fQ(v)  is  the 
Maxwell-Boltzmann  function: 


fo(v) 


e 


-av 


7 


a  =  m/2KT 


J(d3v)fQ(v)  =  1 


We  neglect  terms  involving  products  of  two  gfs.  The  terms 
g  are 

jdsj  dv2  e.(vrv2)  {f0(vi)f0(v')  ~  WW}  » 
e  v2 

which  vanishes  identically.  For 

(vj)  2  =  -  2  vr(v1*e)  + 

(vp2  =  v2  +  2  vr(v2*e)  + 


and 

(vp  2  +  (v p  2  =  v2  +  v2 

since 


(vl  -  v2> 


e  =  v. 


(8-10) 

normalized 


(8-11) 

without 

(8-12 


(8-13) 
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Hence , 


£o(vi)£o(v2>  -  “WW  • 


(8-14) 


Therefore , 


^3t  +  V1  ’  ^9^1^  ~  ~2  (2a)2  J  dQ  J  dv2  |e*  (v1-v2)|  fQ(v2) 


e  v. 


x  jg(rv^t)  +  g(rv^t)  -  gtrv.^)  -  g(rv2t)j  (8-15) 

in  which  we  have  used  (8-14)  again.  The  r  and  t  transforms  of 
g(rtv)  ,  which  we  call  G^(v)  ,  satisfies: 

(co  -  v1>k)G(v1)  =  i  (2a)2  J  dQ  J  dv2  |e •  (v1~v2)  |  fQ(v2) 

e  v„ 


|g(v|)  +  G(v^)  -  G(v1)  -  G(v2)j 


(8-16) 


or 


(oo  -  v^*k)G(v)  =  i  J  G(v)  (8-17) 

in  which  J  is  the  linear  integral  operator  acting  on  G. 

There  is  a  formal  approach  to  (8-17)  that  is  very  suggestive. 
Let  us  suppose  that  there  exist  orthonormal  eigenfunctions  of  J, 
and  let  them  be  \J r  with  X^  the  eigenvalues.  Then 

Ctyn(v)  =  Kntn(v) 


where 


(dv)f  (v)\|f  (v)f  (v)  =  f> 

«  T  n  T  m  o  mn 


(8-18) 
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There  are  five  functions  of  velocity  that  automatically  make 
the  collision  term  vanish.  They  are,  apart  from  constants  and 
orthogonalization , 

g  =  1 ,  a  constant 

g  =  v,  the  three  components  of  velocity 
2 

g  =  v  ,  the  kinetic  energy. 

Hence,  there  are  five  zero-eigenvalue  eigenfunctions  of  J.  The 
remaining  eigenvalues  are  negative.  For 

,  J(dv1)fo(v1)|i(v1)  J  [(ti(v1)] 

J  (dVl) fo^i 


and  the  numerator  can  be  written 

J^vif^J^)  =  -1/4  ( 2a)  J  dQ  J  dVj^  J  dv2  |e- (v^v^  |  f^v^  fQ(v2) 

e 

x  +  ^(v^)  -  |i(v1)  -  ti(v2)} 

We  get  this  result  by  first  interchanging  v^  and  v2*  then  interchang¬ 
ing  v^v2  and  vj^  ,  and  taking  the  average  of  the  four  results. 

That  most  of  the  eigenvalues  are  negative  means  only  that  the 
collision  term  causes  damping  of  most  disturbances. 

The  only  disturbances  that  the  collision  term  cannot  influence 

2 

are  changes  in  density  (1) ,  momentum  (v)  and  energy  (v  )  ,  i.e., 
temperature,  of  the  gas. 

If  we  knew  the  orthonormal  eigenfunctions,  and  if  we  knew 
that  they  formed  a  complete  set,  we  could  expand  an  arbitrary  G(v) 
in  a  series: 
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G(v)  =  >  c(ku)  t_(v) 
Li  m  T  m 

m 


(8-19) 


Then  (8-17)  becomes 

Y  c  (oo-v  *k  -  i U|m(v)  =  0  (8-20) 

l-i  m  m  7  m 

m 


or 


(oo  + 


i\  )  c 
m  m 


-  k 


i 


n 


v  c 
mn  n 


=  0 


(8-21) 


where 


v  =  I  (d3v)\J/  (v)v^  (v)f  (v)  (8-22) 

mn  J  7  m  7  n  o 

is  a  matrix  element  of  v. 

The  dispersion  relation  for  wave  motiones  in  the  gas  are  given 
by  the  condition  that  (8-21)  have  a  solution: 

det(o)  +  i\  )  6  -  k  •  v  =  0  (8-23) 

m  mn  mn 

The  problem  of  the  completeness  of  the  f  has  not  been  solved , 
neither  in  general  nor  for  hard-sphere  collision  terms.  Nor  is  it 
even  known  that  the  spectrum  of  kfs  is  discrete,  although  both 
completeness  and  discreteness  are  entirely  plausible.  Still, 

Eq.  (8-23)  is  very  revealing.  Let  us  restrict  ourselves  only  to 
the  five  eigenfunctions  for  which  \  =  0.  If  we  solve  (8-23)  for 
this  very  restricted  subspace  of  the  full  function  space  of  the 
>]/m,  we  find  —  with  Wang,  Chang  and  Uhlenbeck  —  the  dispersion 
relation  for  sound  waves. 
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Since  fQ(v)  is  normalized,  ^  =  1  is  also  normalized  (relative 
to  fQ(v))  .  ^3  and  are  porportional  to  vx  v  vz,  respectively. 

They  are  therefore  orthogonal  to  \|/^.  For  normalization  we  find 


tv  =  \[2a  v 

2  2  3 

Now  v  is  not  orthogonal  to  ^ ,  but  v  ~  Ja  is#  T^ie  normal^zed  func¬ 
tion  is 

f5  =  a  ^2/3  (v2  - 


The  matrix  elements  of  vx  are 

1  2  3  4  5 
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Or 


3,  2 

00  (go 


5k 


6a 


0 


But  a  -  m/2KT  ,  so  that 


GO 


2 


5  KT  k2 
3  m  x 


By  symmetry, 

ca2  =  f  SI  k2  (8-25) 

o  m 

Thus,  there  is  every  hope  that  a  purely  microscopic  approach 
to  the  problem  of  wave-motions  in  a  multi-particle  system  will  suc¬ 
ceed,  and  without  undue  analytical  complications.  The  alternate 
approach,  the  Macroscopic  methods  are  best  suited  to  the  formulation 
of  boundary  conditions  when  such  conditions  are  easier  to  obtain 
empirically  than  theoretically. 

To  continue  with  the  solution  of  the  linearized  Boltzmann 
equation,  let  us  note  that  the  nonzero  collision  terms  are  numeri¬ 
cally  very  important.  Returning  to  Eq.  (8-15) ,  we  can  estimate 
the  magnitude  of  a  typical  term,  the  third  on  the  right-hand-side. 
Introducing  transforms,  bringing  it  to  the  left  side,  and  replacing 

vi~v2  ky  v'  and  vi  ky  v>  we  9et 

go  -  v  •  k  +  nQ(2a)2  J  dvf  J  dQ  |e  *v '  |  f  (v-v 1 )  G  (v)  =  ••• 

v T  e 


Now 


J  dQ^  |e*vf|  =  27r  fv f  | 

as  can  be  seen  by  letting  vT  be  along  the  polar  axis  of  the  e- 
coordinates  in  a  spherical  system.  At  zero  temperature 
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( 3) 

f  (v-v')  =  6  (v-v') ,  and  the  complex  factor  contains  the  fre¬ 

quency 

v  =  Trno(2a)  2  (v|  (8-26) 

Now,  we  do  not  expect  a  sound-wave  at  zero  temperature,  but  the 

quantity  of  (8-26)  will  be  typically  important  as  a  measure  of  the 

collision  rates.  To  get  an  idea  of  the  difficulties  that  are  in 

store  in  any  microscopic  treatment  of  gases ,  consider  the  magnitude 

of  v  for  air  at  room  temperature.  To  sufficient  approximation  we  can 

—8 

treat  a  nitrogen  (N0)  gas.  let  2a  v  3,5  x  10  cm  be  a  typical  atomic 

-24  3 

diamter,  m  ^  43  x  10  gms  per  molecules;  a  density  of  0.00129  gm/cm 

leads  to  n  ~  3  x  10^  molecules/cm^;  a  temperature  of  300°K  leads 
°  4 

to  a  thermal  velocity  of  4.6  x  10  cm/sec;  so  that  a  thermal  mole¬ 
cule  experiences  a  collision  frequency 

v  4  x  10^/sec  (8-27) 

3 

With  such  a  value,  how  can  a  sound  wave  of  10  cps  frequency  travel 
more  than  about 

vth  -5 

d  ~  ^  10  cm  (8-28) 

For  v  in  G(v)  would  completely  dominate  go,  in  fact  wiping  out  all 
reference  to  col  For  this  not  to  happen  in  practice,  the  remaining 
terms  in  (8-13)  must  cancel  off  most  of  this  damping  effect.  What 
happens  in  practice  is  this:  a  disturbance  of  the  air,  g(rvt)  ,  can 
be  thought  of  as  the  superposition  of  characteristic  disturbances, 
or  ’’eigen-disturbances”.  These  include  changes  that  leave  a  perman¬ 
ent  alteration  in  the  density,  momentum  and  energy  of  the  air 
(X  =  0)  *  together  with  all  other  possible  changes  in  the  local 
velocity  distribution,  which  are  transient  (X  ^  0)  .  The  latter 
are,  in  fact,  rapidly  damped,  but  the  former  survive  as  long  as 
the  conservation  laws  permit.  A  locally  produced  sound  dies  away 
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with  distance,  since  the  "density-momentum-energy"  is  shared  over  larger 
and  larger  surfaces  with  time.  But  an  ideal  plane-wave  remains  intact 
because  each  element  of  the  advancing  wave  front  receives 
equal  contributions  of  "density-momentum-energy"  from  all  sides. 

Sound  waves  will  exist  in  a  plasma  also,  as  can  be  understood 
from  the  great  generality  of  the  processes  involved.  But  more  often 
than  not  the  excitations  of  the  plasma  that  are  of  most  interest  are 
those  for  which  \  /  0 ,  and  for  which  a  knowledge  of  (and  ty/)  are 
of  the  greatest  importance.  In  this  report  we  have  dealt  with  hydro- 
magnetic  non-damped  wave-motions  derived  from  a  model  of  the  plasma 
which  ignored  collisions  entirely.  That  we  were  justified  can  be 
argued  only  by  showing  that  indeed  collisions  are  unimportant  (at 
zero  temperature) .  The  sequel  to  this  report  will  consider  this 
question.  We  can,  however,  note  briefly  that  unlike  the  hard-sphere 
case,  for  which  a  collision  frequency  exists  for  a  particle  of  non¬ 
zero  velocity  in  a  gas  at  zero  temperature,  there  is  no  intrinsic 
radius  of  repulsion  for  electrons  and  ions.  In  fact,  the  long- 
range  forces,  and  the  mixture  of  repulsion  and  attraction  in  the 
plasma,  produces  a  variety  of  collective  phenomena  which  the  hard- 
sphere  model  can  in  no  way  duplicate.  The  search  for  solutions  to 
the  problem  of  the  wave-motions  in  a  plasma  promises  great  rewards. 
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